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ABSTRACT 


Atmospheric winds at heights between 25 and 120 km have been retrieved 
with precisions of 5/ms - 1 from the Doppler shifts of atmospheric absorption 
lines measured from a satellite-borne instrument. Lines of the 1/3 CO2 and 
1/2 H2O rotation- vibration bands caused by gases in the instrument allowed 
the instrumental frequency scale to be absolutely calibrated so that accurate 
relative speeds could be obtained. By comparing the positions of both sets of 
instrumental lines the calibration of the frequency scale was determined to be 
stable to a precision of < 2 x ldf 5 cm -1 during the course of each occultation. 
It was found that the instrumental resolution of 0.015/cm _1 after apodization, 
the signal to noise ratio of about 100 and stable calibration allowed relative 
speeds to be determined to a precision of 5 /ms" 1 or better by using small 
numbers of absorption lines between 1600 and 3200^/cm -1 . Absolute absorption 
line positions were simultaneously recovered to precisions of 5 x 10 -5 cm -1 or 
better. 


The wind speed profiles determined from four sunset occupations and 
one sunrise occultation show remarkable similarities in the magnitudes and 
directions of the zonal wind velocities as functions of height. These wind profiles 
appear to be manifestations of atmospheric tides. The periodic heating of the 
atmosphere primarily by ozone absorption and water insolation can force global 


circulation of the upper atmosphere as described by Lindzen and Chapman 
(1970). The diurnal and semidiurnal Fourier components of the heating rates 
as predicted by Groves (1983) were shown to be significant in driving the tidal 
circulation. The linear or “classical” tidal theory [Lindzen, Chapman, 1970], 
however, predicts wind speeds which grow exponentially with height due to 
the exponential decrease in the hydrostatic background pressure. The growth 
was found to lead to wind speeds much larger in amplitude than revealed by 
the spectroscopic measurements. The classical theory must break down in a 
nonlinear way and the consequences of nonlinear and adiabatic instabilities 
present in such growth in amplitude is discussed. 
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INTRODUCTION 
Mission Objectives 

The Atmospheric Trace Molecule Spectroscopy experiment (ATMOS) was 
designed to obtain high resolution infrared absorption spectra from an orbit- 
ing platform in space. The instrument was managed for the National Aero- 
nautics and Space Administration (NASA) by the Jet Propulsion Laboratory 
(JPL) under the direction of C. B. Farmer and was first flown in the spring 
of 1985 as a payload of the shuttle Spacelab-3 mission. This project has also 
become a part of NASA’s Upper Atmospheric Research program (UAR). 
The major goal of the UAR is to obtain quantitative data pertaining to the 
state of the upper atmosphere over a long period of time to determine the ef- 
fects of natural and perhaps man-made changes of the atmosphere’s physical 
state. Information has been obtained from the Spacelab-3 experiment by the 
ATMOS experimental team concerning the volume mixing ratios (VMRs) 
and concentrations of a number of minor and trace gases in the atmosphere. 
Temperature and pressure variations with height were also obtained by ana- 
lyzing the strengths and equivalent widths of the absorption lines and these 
measurements are related to upper atmosphere dynamics. The purpose of 
the research described here was to determine if accurate wind speeds along 
the line of sight can be obtained to high precision. This was achieved. In 
addition the influence of the diurnal and semi-diurnal Fourier components of 
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the heating of the lower atmosphere on the upper atmospheric circulation is 
discussed in this report. 

The ATMOS instrument is a fast Fourier-transform (FFT) interferometer 
which covers a wide spectral range from 600 to 5000 cm -1 . This spectral 
range is further divided into smaller regions by the use of specialized filters 
to increase the signal to noise ratio of the processed spectra. With a scan 
time of about 2 s the spectra have an unapodized resolution of about 0.01 
cm -1 , a digital frequency spacing of 0.0075 cm -1 and a vertical resolution of 
about 4 km. A variety of programs is available on a Prime 955 computer at 
the ATMOS data analysis center to study the spectra. 

The mission’s strategy was that, by using the sun as the source of ra- 
diation, absorption spectra of about 40 molecular species in the earth’s at- 
mosphere could be observed in successive spectra during sunset and sunrise 
solar “occultations”. An occultation, or more precisely a limb observation, 
is defined as when the earth’s atmosphere intersects a line between the or- 
biter and the sun during an orbit. From these spectra the properties of the 
atmospheric constituents and the dynamical state of the atmosphere in re- 
gions near the tangent points of observation were studied. A tangent point 
is defined as the point above the earth’s surface where a light ray from the 
sun to the orbiter is closest to the earth’s surface. This distance above the 
earth’s surface is referred to as the tangent height. A description of the opti- 
cal filters used to define several spectral regions and a list of some molecules 
whose infrared transitions occur in these regions are given in Table 1 [Farmer, 
1987]. 

During the course of the mission two latitude bands were sampled. Mid- 
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latitude sunsets occurred near 30° N and sunrises near 50° S. Since the orbital 
period of the platform was about 90 minutes a wide range of longitudes was 
available for study in the two latitude bands. The Spacelab-3 Mission (SL-3) 
was launched on board the shuttle Challenger on April 29, 1985 at 21:00:30 
EST and Table 2 lists the geographical locations, times, and filter region 
for each occupation. For each occupation a set of spectra was obtained 
in a period of a few minutes and hence a cross section of the atmospheric 
state could be sampled as a function of altitude for the latitude and longitude 
about the tangent points. A total of 19 such occupations accumulated almost 
1200 atmospheric spectra. In addition, more than 1000 solar spectra were 
obtained at altitudes far above those where infrared absorption by telluric 
gases is significant. 


Overview of Thesis Research 

After the completion of the mission in May 1985 the raw data were pro- 
cessed and soon afterwards the first transformed spectra were available for 
analysis. It was decided to study the possibility of measuring winds from 
the spectra and to estimate the accuracy which can be expected. Due to the 
relative motion of the shuttle orbiter with respect to the atmosphere as well 
as possible motions caused by atmospheric circulations the absorption line 
positions were expected to change significantly from spectrum to spectrum. 
It was determined that the relative speeds between the observer and the 
atmosphere surrounding the tangent points could be determined by measur- 
ing the Doppler shifts of groups of absorption lines and the relative speeds 
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associated with these shifts. 

To measure the Doppler shift of a moving absorber, the frequencies of 
two lines must be compared, the rest frame position and the position in the 
moving frame. If c is the speed of light, v is the relative speed of the absorber 
to the orbiter, and c is the wavenumber (in cm -1 ) of the line in the rest frame, 
then the new wavenumber in the moving frame is shifted from the rest frame 
position by the amount 

Ac=-?. (1) 

c 

The relative speed v can be determined by this shifts in frequency by 

A? 

v = — c. 

C 

The accuracy of v is limited by how well known the rest frame position t is 
relative to the instrument’s “true” rest frame position for that transition. It 
was determined that even if the instrumental calibration is not “absolute”, 
and hence line positions determined using it were not absolute, accurate wind 
speeds could, nonetheless, still be determined. This is true since v depends 
not on the Doppler shift Ac alone, but on the ratio Ac/c. Average calibration 
offsets will cancel in the determination of v and this is discussed in §11.4. This 
technique was found to require that for at least one telluric absorption band, 
a set of corresponding lines due to gases inside the instrument needed to be 
present in the spectra. The accuracy was then transferred to other spectral 
regions by a “bootstrapping” technique which utilized the wind speeds de- 
termined from bands in which instrumental lines were available to determine 
the positions of new absorption lines relative to the instrument’s calibration. 

It was determined that a sufficient number of spectral lines was available 
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in the spectra (~ 20-40 lines/spectrum) to obtain precisions of 5 ms -1 or 
better by measuring the Doppler shifts of telluric absorption lines [Shaw, 
1985] using realistic estimates of the noise level. It was also determined that 
upper limits of fluctuations in the calibration stability of the instrument were 
needed since the precisions of the wind speeds would be ultimately limited 
by these. Such uncertainties were determined to be extremely small and did 
not contribute more than 1-2 ms -1 errors. 

Accurate relative speed profiles were estimated from spectra with tangent 
heights between 25 and 120 km and were found to have precisions from 3 to 
5 ms -1 throughout this range. In order to obtain wind speeds relative to the 
rotating frame of the earth the orbital motion of the shuttle and the rotational 
motion of the tangent points were estimated from telemetry data provided by 
Goddard Space Flight Center. This motion was subtracted from the relative 
speeds between the orbiter and the atmosphere estimated from the Doppler 
shifts to obtain wind speeds relative to the earth’s surface. Uncertainties due 
to errors in estimating the line of sight velocity components were found to 
be less than 2 or 3 ms -1 . The residual wind speed profiles were found to 
indicate essentially the zonal component of the atmospheric flow. 

Similarities of the wind speeds with height were found in all four sunset 
occultations analyzed. These profiles were found to have wavelike character- 
istics and further investigation of this phenomenon was done. (Since only 
one sunrise occultation was available in the filter used for wind speeds such 
comparisons with other sunrise measurements are unavailable at the current 
time.) Such variations in altitude can be described by the theory of atmo- 
spheric tides [Siebert, 1961; Chapman, Lindzen, 1970]. The source of the 
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excitations and their “wavelike” nature can be explained by the periodic ab- 
sorption of solar radiation by the atmosphere, particularly by ozone between 
20 and 90 km and water vapor below 10 km altitude. Since this happens on a 

global scale the response of the atmosphere is also on this scale. The notion 
* 

of tidal eigenfunctions and accompanying vertical structure equations were 
derived and solutions to these mathematical models were computed. There 
were similarities of the phases predicted between the theory and the ATMOS 
results. The “classical” theory of tides, however, was found to be inadequate 
in describing the magnitudes of the waves at heights above about 85 km since 
it ignores the contributions due to non-linearities and adiabatic instabilities. 
The wind speeds measured do contain information of the geometric distribu- 
tion of global heating rates regarding the absorption of solar radiation in the 
atmosphere. This can be attributed to the global scale of the atmospheric 
response to tidal forcing. Further study into the non-linearities in the fluid 
dynamics as well as future measurements of winds in a manner similar to the 
ATMOS project may yield important data about the geographical distribu- 
tion of the gaseous absorbers. 
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Table 1 Frequency transition ranges of the ATMOS optical filters together 
with the principal molecular species that can be analyzed in these regions 
[Farmer, 1987]. 


Filter 1 
600-1200 cm " 1 

Filter 2 

1100-2000 cm " 1 

Filter 3 

1580-3400 cm " 1 

Filter 4 

3100-4700 cm -1 

C0 2 

C0 2 

COj 

C0 2 

HjO 

HjO 

H a O 

HjO 

o 3 


0 3 



ch 4 

ch 4 




CO 




NjO 


nh 3 

H 2 Oj 

NO 

HF 

HNOa 

HOj 

(NOa) 

(HCN) 

(HN0 3 ) 




hno 4 




cci 3 f 

(NO) 

HDO 


CCLjFj 

NOj 

(HaOa) 


CHCC1F2 

N 3 O 5 

HjCO 


CH 3 C1 

HN0 3 



CC1 4 


HC1 


cof 2 

HOC1 

(CH 3 C1) 


COC1F 

cf 4 

OCS 


CIO 

S0 2 

HCN 


ClONOj 
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Table 2 Summary of occultation spectra obtained during the Spacelab-3 mis- 
sion. Latitudes and Longitudes refer to both high (~ 120 km) and low (~ 
0 km) tangent heights. MET - Mission Elapsed Time, GMT - Greenwich 
Mean Time. 


Occ 

Type 

Opt 

Fil 

MET 

Start 

GMT 
60 km 

Lat 

Low/High 

Long (E) 
Low/High 

SR01 

3 

0/07:13 

119/23:17 

49.9/48.5 S 

118.0/117.0 

SS01 

1 

0/08:08 

120/00:12 

34.5/35.9 N 

276.3/277.1 

SS02 

1 

0/18:52 

120/10:54 

32.6/34.1 N 

115.3/116.1 

SR02 

3 

0/19:27 

120/11:20 

49.1/47.5 S 

294.6/293.3 

SS03 

2 

0/21:55 

120/13:57 

33.6/31.9 S 

70.2/69.3 

SS04 

4 

1/00:55 

120/17:00 

33.0/31.3 N 

24.2/23.3 

SR03 

6 

1/01:34 

120/17:36 

48.6/46.9 S 

202.7/201.4 

SS05 

1 

1/02:28 

120/18:32 

32.8/31.1 N 

1.2/0.3 

SS06 

3 

1/07:01 

120/23:07 

31.9/30.1 N 

292.1/291.2 

SS07 

2 

1/08:33 

121/00:38 

31.7/36.0 N 

269.2/268.2 

SS08 

1 

1/10:05 

121/02:10 

31.4/29.6 N 

246.2/245.2 

SR04 

4 

1/15:18 

121/07:21 

47.4/45.6 S 

356.0/354.7 

SR05 

2 

1/16:49 

121/08:53 

47.3/45.5 S 

332.6/331.8 

SR06 

1 

1/18:21 

121/10:24 

47.2/45.4 S 

310.1/308.8 

SR07 

5 

1/19:53 

121/11:56 

47.1/45.1 S 

287.2/285.8 

SS09 

3 

2/01:22 

121/17:26 

28.5/26.6 N 

16.1/15.2 

SS10 

4 

2/02:53 

121/18:58 

28.1/26.2 N 

353.1/352.2 

SS11 

3 

2/04:25 

121/20:30 

27.8/25.9 N 

330.1/329.2 

SS12 

2 

2/05:51 

121/22:01 

27.3/25.3 N 

284.1/283.2 

SS13 

3 

2/07:27 

121/23:33 

27.3/25.3 N 

284.1/283.2 
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CHAPTER I 

SPECTROSCOPIC THEORY 


1.1 The Michelson Interferometer and Fourier 
Transform Spectroscopy 


In the classical Michelson interferometer Figure 1.1.1 light from the source 
oo is collimated at C into a beam which is incident on a half-silvered beam 
splitter B. The reflected ray travels to a stationary mirror S and it is reflected 
back to the beam splitter. The partially transmitted ray travels to a movable 
mirror M and is reflected back to the beam splitter. The two split beams 
recombine at O and are sent to the detector D. Since M is movable, the two 
beams interfere both constructively and destructively as the path difference 
x changes. 

For simplicity we consider an incident beam with frequency /. The 
wavenumber ? is related to the frequency / and the speed of light c by 





( 1 . 1 ) 


If we assume the amplitude A of this monochromatic incident beam is split 
evenly between both arms M and 5 with path difference x and if we neglect 
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losses in intensity in the beam due to dispersion then the complex amplitude 
of the beam at the detector is 

A d = ix(l + e“), (1.2) 

where the phase difference is 

6 = 27r?x. (1.3) 

The intensity at the detector D is 

/(s,x) = | Ad\ 2 = ^(1 + cos£)A 2 = ^(1 + cos£)J 0 (c). (1.4) 

Suppose that 1(c) represents the spectral radiance (power/cm -1 ) then 
total signal J 0 (x) received at the detector D is the sum of all contributions 
at all frequencies 

M x ) = l dq i(l + cos £) /(c), 

Jo & 

= i J a dfl a (f) + \J 0 * cos 6 1(f), 

= \ [/ + /(*)] • ( 1 - 5 ) 

The first term / is the total incident intensity while the second term /(*) is 
the modulated part and is given by 

J(x) = 2/ 0 (x) - 7. (1.6) 

In Fourier spectroscopy the interference pattern J(x), also known as an *n- 
terferogram , is measured and from it the the spectral radiance 1(c) is deter- 
mined. From the theory of Fourier transforms we have the relationship 

/ oo 

dc ; COS 27TCX /(c0> 

-OO 


(1.7) 
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where it is assumed that /(?) = 0 for ? < 0. This can be inverted [Mathews, 
Walker, 1970] to give 

/ oo 

dx cos 27r?x J(x) (1.8) 

■OO 

and the pair of equations (1.7-8) are a Fourier-cosine pair. This technique of 
spectroscopy has the multiplex advantage which improves the signal to noise 
since all spectral features are simultaneously recorded in the interferograms 
[Thorne, 1988]. 

In practice the maximum path difference is not infinite but is usually 
confined to a limited region say, —L < x < L. Then J(x) is zero outside this 
range and the spectrum (1.8) is 

r L 

/'(?) oc / dx cos 2ffc® J(x). (1.9) 

J —L 

Some spectral information is lost due to the finite path length with a subse- 
quent degradation of the instrumental resolution. If we consider the case of 
a monochromatic input signal Jo at wavenumber 

/(<) = /(>%-<')> ( 1 . 10 ) 

then the output interferogram is a cosine function 

J(x) = J 0 cos 2tt<;'x. (1.11) 

Since data are only available for — L < x < L the computed spectrum I\<;) 
is 

f L 

/'(?) OC Jo / dx COS 27T?X COS 27lVx 

J — h 

oc I 0 L [sine 2ir(?' + c )L + sine 2 tt(c > — s)L) (1.12) 
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where the sine function is defined as 


smey = 


siny 


(1.13) 


Since the first term on the right hand side of (1.12) in negligible compared 
to the second the output interferogram for the monochromatic wave is 


I'(c) <x I 0 L sinc27r(?' — c)L. 


(1.14) 


Figure 1.1.2 shows the first two zeros of the sine function ±1/2L. The reso- 
lution of the spectra is approximately 1/2 L and corresponds to the region of 
maximum instrumental response to a spectral feature. Outside this region 
the sinusoidal variation of the instrumental response is due to the finite path 
difference. As L ' — ► oo the sine function approaches a ^-function, and we 
retrieve all spectral information. 

Now consider the case where the radiation is described by a spectral 
density 1(c) 

i(c) = r dc' %-cW)- (i-i5) 

J — oo 

Then, from (1.7), the interferogram is 

I(x) <x f dq cos 2ircx f dc' I(c') % - ?')> 

J —oo J —OO 

= f dc' cos 2ttc'x /(?'). (1-16) 

J — oo 

Since data are only gathered for — L < x < L the resolved spectrum is, from 
(1.9) 

r L 

I'(c) oc J ^dx cos27rcx I(x), 

= J dc' J dx COS27TC* cos 2 tc'x I(c'). 


(1.17) 
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We see that (1.17) can be written 

(i-w) 

J —oo 

if we define the normalized instrument function 

r L 

i(q, c*) = I dx cos 2‘kc'x cos 2ir<;x, 

J —L 

= i [sine 27r(c + q')L + sine 2ir{q — q')L ] . (1.19) 


It can be shown [Thorne, 1988] that this form is normalized 

/ dct'(?,?') = l* (1.20) 

J — OO 

We see that (1.19) represents a response centered around q = ±?'. This is 
true since Fourier transforms of functions in the x domain always have mirror 
images about zero in the q domain. If we define 


*±(c,0 = L sine 2irL(q ± q'), 

(1.21) 

then the instrument function (1.19) is 


*(?,0 = *+(*»0 + *-(*»«')• 

(1.22) 


Usually we refer to the region q > 0 so that the contribution of t + is usu- 
ally negligible for q 1/2 L. We can then say the instrument function is 

approximately 

*(?,?') = »(? - ?') = 2t_(c,?') = 2 L sine 2irL{q - ?'). (1.23) 

The convolution between two functions -A(c) and B(q) is defined as 
-A(s) * £(?) = f dq' A(? -.?')£(«) 

J —OO 


(1.24) 
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so then. (1.18) can be expressed 

/'(?) = (1*25) 

The resolved spectrum /'(?) is the convolution of the true spectrum /(?) with 
the instrumental function i. 

The resolution 6<; = 1/2 L caused by the finite path length allows infor- 
mation from nearby wavenumbers to contribute to the partially resolved 
spectrum at s which is given by (1.25). The sine function has a “ringing” 
at wavenumbers not centered near the maximum of the instrument function 
which can contribute significant interferences between different wavenum- 
bers in the spectra. Figure 1.1.3 shows the ringing caused by this instrument 
function about an absorption line taken from the ATMOS spectra. The deep 
central spike is the instrumental response at the spectral line center while 
the side lobes are artifacts of the finite path difference and not real spectral 
features. This ringing can be mathematically filtered out through a process 
called apodization - literally meaning “cutting of the feet” [Thorne, 1988] - 
but with a complementary loss of spectral resolution . From the property 
of Fourier transforms it can be shown that the transform of a product of 
functions is the convolution of their transforms 

[A£](0 = A(c)*£(c). (1.26) 

We can now imagine multiplying the interferograms I{x) by an “apodizing 
function” a(x). An apodized spectrum can be retrieved from this apodized 
interferogram a(x) I{x) using (1.25) 


/'(c) = •( c) * [«/](?) = "(c) * Me) * /(c)]. 


(1.27) 
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If we define the apodized instrument function 

».(?, o = r *(? - o a (?" - o (i-28) 

then by using (1.25) the apodized spectrum (1.27) can be shown to be 

/'(c) = 0/(0- (1.29) 

The apodizing function a(x) can be freely chosen to achieve whatever mathe- 
matical filtering is desired but always at the cost of spectral resolution. There 
are several apodizing functions available to the user when computing the AT- 
MOS spectra and they are described in detail by Norton and Beer (1976). 
Figure 1.1.4 shows the absorption line used in Figure 1.1.3 after convolution 
with three different apodization functions available on the ATMOS computer 
and described by Norton and Beer [1976] were used. These apodizing func- 
tions give more smoothing of the spectral while degrading the resolution from 
0.0125 cm -1 to about 0.015 cm -1 between the several functions. The apodiz- 
ing functions widen the spectral features and cause the peak absorptances of 
the lines in the spectra to decrease. The apodizing function chosen in this 
study eliminated false spectral features greater than the noise level and had 
a resulting resolution of about 0.015 cm -1 . This was found to be acceptable 
within the measurement requirements §1.12. 

If an interferogram I(x) is sampled at intervals Ax features at other 
wavenumbers may be superimposed on real features unless the value of Ax 
is properly chosen. This “aliasing” is due to the periodicity of the Fourier 
transform and is described by Thorne (1988). Thorne shows that the inter- 
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ferogram must be sampled at least every 

A* = 5 ^- (1-30) 

*Smin 

if the minimum wavenumber is to avoid this sort of overlapping. Special 
optical filters must also be used to prevent false lines horn appearing in the 
spectra. The ATMOS filters were carefully chosen to avoid aliasing [Farmer, 
1987]. 

Absolute line positions, are ultimately limited by the accuracy of the He- 
Ne laser wavelength used to calibrate the instrument. If the recorded value 
of the path difference x * is related to the true path difference * by 

x* = (1 + a)x (1.31) 

where |a| <C 1 then the interferogram (1.7) is 

/(l) = 1 (its) = £L *' cos 2 ^' (its) t 1 - 32 ) 

This shows how the interferogram is related to the path difference recorded by 
the instrument. These interferograms are then processed by (1.9). However, 
since the path difference is uncertain, we obtain a spectrum 

/'(?) = f dx * cos 27TCX* /(®) 

J — It 

= / dc'J(c') f dx* cos2xc®* cos2irc' — ^ — . (1.33) 

J — oo J —L 1 -(- tt 

The instrument function is then 

»a(c, o = L [sine 27T (q + (1 + a)c') L + sine 2t (c - (1 + a)?') L] . (1.34) 

From this result the uncertainty in the laser frequency causes the path lengths 
to be in error which causes features at frequency ? in the convolved spectra to 
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be shifted by — etc. Typical errors in the He-Ne frequency are given by c(a) ~ 
0(1O -8 ) [Norton, 1989] and line positions determined from the spectra should 
be uncertain to this degree. It is shown in §11.4 that these uncertainties do 
not affect the accuracy of the winds from Doppler shifts provided that a set 
of instrumental lines with atmospheric counterparts can be measured in the 
spectra. 

Due to the digital nature of the interferogram data and the presence of 
noise in the detector signals the location of zero path difference (ZPD) in the 
interferogram can be in error. The point of zero path difference is defined as 
that value of x for which all frequencies interfere constructively. For an ideal 
interferogram this point is taken as the origin of x. Errors in locating ZPD 
can introduce apparent asymmetries in the interferograms which may cause 
frequency dependent calibration errors in the resolved spectra /'(c). To see 
this, suppose that an estimate of ZPD is in error by an amount £. Then the 
measured interferogram (1.7) is 

J(x) oc f dc ; cos 27 tc'(x + () Z(c'). (1.35) 

J — OO 

Using the form of this interferogram which is sampled from — L < x < L we 
have 


f L 

/'(c) oc J dx cos 2 ttcx /(x) 

= f dq r Z(c') f dx cos 2irc® cos 27tc'(s + £) 
J — oo J—L 

= r dc'x(c,c')/(c'). 

J —oo 


(1.36) 


The instrument function contains the phase error and can be evaluated using 
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standard integration 


*(s>0 


dx cos 2tt?x cos27rc'(x + £) 

= L [sine 27r(c — $')L + sine 2?r(c + c ')L\ cos 27rc'£ 
~ L sine 2ir(s — ?’) L • cos 2 



(1.37) 


The instrument response function »(?, ?') is modulated by the term cos 
For values of <;'£ <t. 1 the interferogram will peak when ?' = ? and should 
produce symmetric line shapes. For large frequencies where ~ 1 the 
instrument function will not peak at ? and the phase errors can produce 
significant distortions in the line shapes which may cause systematic errors 
in the positions of spectral features. 

An estimate of the magnitudes of these systematic fluctuations in line 
positions set a limit of the accuracy of the relative speed measurements. The 
sampling interval Ax of two fringes of the He-Ne laser used for calibration 
was used to obtain the interferograms [Norton, 1989]. The wavelength of 
the laser is about 633 nm which corresponds to Ax ~ 1.27 x 10 -4 cm. The 
precision of determining ZPD was estimated to be 1000 times better than this 
value [Norton, 1989] which corresponds to ?(£) ~ O(10“ 7 ) cm. For typical 
infrared transitions the distortions in the instrument function are negligible. 
With this degree of precision distortions the phase errors are negligible and 
we can therefore expect the calibration of the instrument to be extremely 
stable. 



Figure 1.1.1 Simple Michelson interferometer. The ATMOS interferometer 
has a much more complex optical arangement [Farmer, 1987], however, the 
basic interferometric principles are the same as in this simple design. 



Figure 1.1.2 Instrumental response function sine®. The amplitude of the 
oscillation decreases as |®| increases. 
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Figure 1.1.3 Absorption line from ATMOS spectra without apodization. No- 
tice the sine function behavior of the instrumental response function away 
the line center where it is a maximum. 



Figure 1.1.4 Same absorption above with apodizing functions 1, 2 and 3 
[Norton, Beer, 1976). The figures have been displaced vertically for clarity. 
Notice that the ‘ringing’ in Figure 1.1.3 has been eliminated at the expense 
of instrumental resolution. The peak absorptance of the line is also smaller. 
Apodization ‘2’ was used in this study. 
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1.2 Measuring Doppler Shifts from ATMOS Spectra 


Due to the geometry of the experiment the measurements are confined 
to the components of atmospheric motion along the line of sight. In order 
to determine the components we need to measure to total relative motion 
between the gas near the tangent points and the orbiter and then subtract 
from these results motion due to the shuttle’s orbit and the earth’s rotation. 
The residual speed is then the component of wind relative to the earth’s 
rotating frame along the tangent ray path. In order to obtain the relative 
speed it is useful to observe the actual positions of absorption lines in the 
ATMOS spectra and compute the Doppler shifts of their wavenumbers from 
a rest frame position. The observed position Cj of a feature in the moving 
frame is related to its position cj- in the rest (non-moving) frame by the 
non-relativistic Doppler shift 

5j=(l + ;)ijS(l+m- (1.38) 

Here v is the relative speed between the orbiter and the absorber, c is the 
speed of light, and /3 = vjc which is O(l0 -8 ) for orbital speeds. In order to 
obtain winds to 5 ms -1 we need to be able to determine c to about 0.00005 
cm -1 . Equivalent precisions can be obtained by using many absorption lines. 
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1.3 Quantum Mechanics of Infrared Spectroscopy 


The ATMOS interferometer is designed to measure the absorption of solar 
radiation in the infrared. In this section we examine the physical processes 
involved and describe the quantum mechanical basis for the spectra. 

The energy of a molecule depends on the locations and motions of its 
nuclei, center of mass, orientation, and electronic state. Let the locations 
of all the nuclei relative to the molecule’s center of mass R™ be described 
by the vector R. The electronic energy of the molecule is then described by 
f?(R). This energy is evaluated with stationary nuclei. However, due to the 
uncertainty principle 

"«„,) <PL.) > J (1-39) 

no real molecule can have fixed nuclei t since the more accurately their po- 
sitions R* are determined the larger is the uncertainty in their momenta 
P*. The nuclei therefore “oscillate” about an equilibrium configuration R eq 
in accordance with the uncertainty principle. Also, the molecule cannot be 
considered to have a fixed center of mass R OT and from the uncertainty prin- 
ciple it has some translational momentum. The molecule also cannot have a 
fixed orientation since the uncertainty in its angular momentum would then 
be infinite. All these three types of motion must be considered as part of the 
wave function. 

The nuclear motions can be treated within the context of the Born- 
Oppenheimer approximation [Berry, Rice, Ross, 1980]. The lighter electrons 
are considered to move much more rapidly than the heavier nuclei. For every 
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configuration of the nuclei during their slower scale motions, the electrons 

A 

are considered to attain an “equilibrium” state. The Hamiltonian 7f of the 

A 

molecule can be assumed to be the sum of the nuclear kinetic energy 7 f n and 
the energy of the electronic state 

H = ‘H n + 'H t . (1.40) 

Let the positions of all the electrons relative to the center of mass R. cm be 
given by the vector r. Then Born-Oppenheimer approximation states that 
the wave function ^(r, R) of the molecule is separable into a nuclear and 
electronic wave function 

^(r,R)=«R)^.(r,R). (1.41) 

The electronic part satisfies the relation 

7f e y> e (r, R) = E( Rty e (r, R). (1.42) 

By using the Schrodinger equation 

(1-43) 

it can be shown [Berry, Rice, Ross, 1980] from (1.41-3) that 

V> e (r, R)74V>„(R) + V’n(R)WnV'e(r, R) = V>c(r, R) [E - E{ R)] i> a (R). (1.44) 

The second term on the left hand side is neglected since there are only 
small interactions between electronic and nuclear motions. This is the Born- 
Oppenheimer approximation and (1.44) can then be expressed as 


[«„ + £(R)] V»(R) = £<*n(R). 


(1.45) 
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The quantity 12 (R) is the effective potential energy field caused by the elec- 
trons in which the nuclei move. 

Now assume at that the energy is separable into three parts, translation 
of the center of mass, internal vibration, and rotation of the molecule’s orien- 
tation. (Even though this is not quite accurate and there are corrections to 
this assumption, it will suffice for the following discussion). Then the nuclear 
Hamiltonian can be written as 

W n = 7ftr^ + W T ib + Wrot. (1-46) 


Let the total mass of the molecule be M, the individual nuclear masses be 
77i j, the center of mass momentum operator be P™ 1 , the nuclear momenta 
relative to the center of mass operator in the “molecular coordinate frame” be 
Pi, the angular momentum of the molecular coordinate frame relative to the 
observer’s frame operator be L, the molecular moment of inertia be I, and the 
set of angles O needed to specify the orientation of the molecular coordinate 
system in the observer’s frame. The moments of inertia are evaluated by 


I x = 53 m,( 2 /? + 2, 2 ), 

I 

Iy = X>i(z? + Z?), 
i 

i* = + y^). 


The the Hamiltonians in (1.41) can be expressed as 


7^truu 


n ri b = 


(pcm)2 

2M ' 



+ £(«•), 


(1.47) 


(1.48) 
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The wave function is 


. LI L* L\ 

Htoi 2I X + 2J V + 21 z ' 


V>n = lj> truuV’intV'rot 


and the energy becomes 


E = Etim, + E y ib + E iat . 


( 1 . 49 ) 


( 1 . 50 ) 


The three eigenvalue equations for the translational, vibrational and rota- 
tional energy levels are 


^ftrmn»V , tran»(®-cxn) — 


( P 6 ”) 2 
2M 

■Etr*n*V’tran»(Rcm)j 

?S + £(R) 




‘WvibV’Tib(R-) = 


V’rib(R) 


= Eyi^lfiyi^R), 


( 1 . 51 ) 


^rotV’rot(O) — 


* n 


2 /« 


2L 


L l 

21* 


V’rot(O) 


— ■^rotV’rot(O)" 


The translational energy levels form a continuum when the length scale of 
the motion is large. The occupation of translational energy levels is related 
to the temperature of the gas and, els we shall see, the Doppler shifts due 
to these thermal motions cause the absorption lines to have a width char- 
acteristic of their average motion. The vibrational energy levels Eire discrete 
and are well approximated by the harmonic oscillator model near the minima 
of their vibrational potential wells. Infrared transitions between vibrational 
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and rotational levels are observed in the ATMOS spectra for many atmo- 
spheric gases. The rotational levels have quantized angular momentum due 
to the azimuthal periodicity of the wave function and accompany vibrational 
transitions. 


1.4 Normal Coordinates and Vibrational Energy for 

Triatomic Molecules 


There were about 40 molecular species whose spectra were measured with 
the ATMOS interferometer. In this study C0 2 , H 2 0 and N 2 0 were used to 
obtain wind speeds. We will therefore concentrate on the quantum mechanics 
of vibration of triatomic molecules which describe the spectra of these three 
molecules. 

A molecule with N atoms has 3 N independent nuclear coordinates in 
general. However, not all of these coordinates are required in describing the 
internal molecular energy E Y n,. In the absence of external fields the center 
of mass coordinates R cm and the orientation Cl do not affect the vibrational 
energy. For linear molecules such as C0 2 and N 2 0, two angles (spherical 
coordinates) are needed to specify the orientation of the symmetry axis, 
while for non-linear molecules such as H 2 0 three angles (Euler angles) are 
needed to specify the complete orientation. Thus there are n = 3 N — 5 or 
3AT — 6 vibrational degrees of freedom. For linear triatomic molecules this 


27 


corresponds to n = 4 degrees of freedom while for non-linear molecules we 
have n = 3 degrees of freedom. 

In a triatomic molecule the vibrational energy is a function of three or 
four independent generalized coordinates q a , a = l,..,n [Herzberg, 1945]. 
The cartesian coordinates R,- of each nucleus i are related to the generalized 
coordinates by 

R , = R.-(q). (1.52) 


There is a preferred set of these coordinates q in terms of which the vibra- 
tional motion can be expressed as simply as possible, the normal coordinates 
[Herzberg, 1945]. In this approach we assume a well-defined equilibrium ge- 
ometry about which the nuclei perform small oscillations. This equilibrium 
geometry is called a molecular frame. The nuclei are assumed to perform 
oscillations in a potential energy field F(q) = E (R(q)) about an equilib- 
rium position given by q eq . Since this equilibrium position is defined to be a 
potential energy minimum in the context of the Born-Oppenheimer approx- 
imation, then (dV/dq ^) q .« = 0. If we define our normal coordinates to have 
a zero value at equilibrium, q eq = 0, then the potential energy function can 
expressed as 




+ 


(1.53) 


from a Taylor series expansion about the equilibrium point. For sufficiently 
small oscillations all of the extra terms in the expansion (1.53) can be dropped 
as well as V(0) since it is a constant for the equilibrium geometry. This 
harmonic approximation is only useful for the states near the bottom of the 
potential well and additional anharmonic terms need to be included for a 
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more accurate solution elsewhere [Herzberg, 1945]. 

If we define the force constants [Berry, Rice, Ross, 1980] 



then the potential energy in the harmonic approximation is 

^(q) = 

The kinetic energy of the molecule is 


(1.54) 


(1.55) 


X = i5>|Rj 


(1.56) 


It can be shown [Fetter, Walecka, 1980] that, in terms of the normal coordi- 
nates, 




(1.57) 


where the reduced mass matrix is defined by 


m a<T . 



(1.58) 


which becomes a real, symmetric, constant matrix for the equilibrium geom- 
etry. 

The Hamiltonian operator H r a, could be obtained from (1.48) and (1.55), 
however the Schrodinger equation would be very difficult to solve due to the 
coupling k aa t between terms. It can be shown [Fetter, Walecka, 1980] that a 
set of normal coordinates can always be determined so that the matrices k 
and m are diagonalized. In this case the potential energy is given by 


V = 



(1.59) 
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and the kinetic energy is 

r=i? m '* = ?£ t 1 - 60 ’ 

where we use p„ = m a q a . Then the Hamiltonian operator becomes 


7Aib = £?W (l- 61 ) 

a 

where 

'Hrih, tr = 2 m ^ 2^^*’ (1.62) 

The vibrational Hamiltonian can be written as a sum of noninteracting har- 
monic oscillators and the wave function can be written as the product 


^rib(q) = n^o-tea). ( 1>63 ) 

<r 

The energy is then given by the sum 


£ Tib = £ E yih<a (1.64) 

<r 

where 

= ■Erib,fft^Tib,»'. (1.65) 

The operator (1.62) is the Hamiltonian of a simple quantum mechanical 
harmonic oscillator and the eigenfunctions V^lb , v (q a ) are well known. The 
energy eigenvalues are given by 

Erib,<r = {v<r + ( 1 . 66 ) 

where 

(1.67) 
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and 

v„ = 0, 1,2, • • • ; <r = 1, •••,«. (1.68) 

For small displacements the oscillations about the equilibrium configurations 
reduce to those of a set of independent harmonic oscillators of frequency 
Any small amplitude motion is then a linear combination of the harmonic 
oscillators. 

The determination of the types of motion corresponding to each normal 
mode can, in general, be determined from the masses of the nuclei and the 
equilibrium configuration. One method [Berry, Rice, Ross, 1980] for deter- 
mining the normal mode motions is to examine the molecule in the molecular 
frame. Here the sum of the momenta of the nuclei is zero 

E P i = °- ( L69 ) 

t 

Since the coordinate system rotates with the molecule the total angular mo- 
mentum must also vanish 

E L < = °- (1.70) 

t 

Since the particles are in simple harmonic motion the momentum components 
are a quarter of a cycle ahead of the displacements 

P, = ^ = tuvR, (1.71) 

so that (1.69) becomes 

T. ro,Rj = 0. (1*72) 

t 

Since each nucleus i oscillates about an equilibrium position R* q each angular 
momentum is 


L t - = m,R* q xRi 


(1.73) 
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so that 

^mjR' q xRj = 0. (1.74) 

t 

Combining (1.72) and (1.74) a set of relative coordinate displacements can be 
determined for each type of molecular configuration. Figure 1.4.1 and 1.4.2 
show the various vibration modes for the COj and HjO molecule respectively. 

To evaluate the relative vibrational frequencies of the various modes re- 
quires a knowledge of the potential energy surface E(R) and this is not always 
available. However the relative magnitudes of the frequencies can occasion- 
ally be established. Consider the C0 2 molecule in the following example. 
This molecule has four modes, an asymmetric stretching (i/ 3 ), symmetric 
stretching (vi), and two degenerate bending modes (i/ s ). For a given dis- 
placement, the bending modes require less stretching of the bonds and hence 
smaller changes in the bond energy than the stretching modes. Therefore 
the bending modes should have lower frequencies. Also since the C atom 
is lighter than the O atom, and since the asymmetric mode requires larger 
bond stretching than the symmetric mode, the asymmetric stretching mode 
should have a higher frequency. Thus the frequencies should be ordered as 

c(asymm stretch) > ?(symm stretch) > c(bending). (1*75) 

This is in fact true for the COj molecule where it is found spectroscopi- 
cally that c(asymm stretch) = 2349cm -1 , ?(symm stretch) = 1388cm -1 , and 
?(bending) = 667cm -1 . These are referred to as the v 3 , V\, and 1/2 vibration 
modes respectively. 

Each vibrational mode can be excited to many energy levels. In the 
harmonic oscillator an infinite set of energy eigenvalues axe spaced at inter- 
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vals of fiu>. In real molecules, the higher vibrational energy levels stretch 
the molecule. At farther separations the molecule breaks apart because the 
quantum mechanical binding energy of the electron cloud becomes weaker 
than nuclear repulsion [Berry, Rice, Ross, 1980]. The bound energy levels 
become more closely spaced as the ibational quantum number increases since 
the average curvature of the potential well represented by the k„ well is re- 
duced. This due to deviations of the binding energy from a pure parabolic 
potential and is shown in Figure 1.4.3. 

Each mode of oscillation may have more than one quantum of vibrational 
energy. In fact, different modes of oscillation may be present simultaneously. 
The modes are classified by the number of vibrational quanta associated with 
the state [see Herzberg, 1945]. As an example, the first excited asymmetric 
vibrational mode of COj is labeled the ‘i/3’ which stands for one quantum of 
the n = 3 vibrational state. Two quanta of the symmetric stretching mode 
is labeled Combinations are designated by a sum. For the vibrational 

state with one quantum in the symmetric stretching n = 1 and two quanta 
in the bending n = 2 modes is labeled V\ + 2j/ 3 . Further elaboration on 
rotational- vibrational states and nomenclature for many molecules is given 
by Herzberg (1945). 
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Figure 1.4.1 Normal mode vibrations for the COj molecule. 



Figure 1.4.2 Normal mode vibrations for the HjO molecule. 
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1.5 Rotational Energy Levels 


In the case of diatomic molecules (such as CO) and linear triatomic 
molecules (such as COj and N 2 O) the rotational energy can be approximated 
reasonably well by the rigid symmetric top model 


I 0 = I m = I y ^ l , 


(1.76) 


where the z axis is defined to be the bonding axis. Then the rotational 
Hamiltonian can be expressed as 



2 I, + 2 U ij ■’ 


(1.77) 


where L 2 is the total angular momentum operator (squared) and the z com- 
ponent of the angular momentum operator is L z whose eigenfunctions are 
the spherical harmonics Yi m (0,<f>) [Berry, Rice, Ross, 1980], 


VY lm {0, <f>) = l{l + l)h 3 Y lm (0, <t>), (1.78) 


L t Y lm {0, <f>) = rnhY lm (0 , <f>). (1.79) 


In the case of the rigid rotor the energy levels axe 

(1.80) 

for l — 0,1,2,... and m = —I, ...,/. 

If the atoms are considered as mass points then there is no z component 
of the angular momentum and the levels are 21+1 degenerate. However, in 
a linear molecule the contributions to I z from the electrons and the non-zero 


E Ioi — 


l v +1 K m > 

; h m 


(K); 
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nuclear radii are both very small, I x <C J 0 , and the coefficient of L\ in (1.80) 
is much larger than that of L 3 . For example, the energy associated with 
L\ is very large when L z ^ 0 in all states except the ground states when 
it varnishes [Goody, 1964; Berry, Rice, Ross, 1980], In any case, this second 
term can be included into the electronic energy since it is a constant and the 
rotational energy can be described by 


„ AJ + 1)* 1 

2h~ 


(1.81) 


for J = 0, 1,2, • ••. The quantum number J is used instead of l in describing 
molecular angular momentum and in the case of linear molecules A h is iden- 
tified as the z component of angular momentum where A is an integer. For 
a linear molecule A = 0, 1,2, • • • and therefore the total angular momentum 
must be equal to or greater than A, J — A, A + 1, A + 2, • • •. Most of the 
infrared rotational- vibrational spectra described in this study are associated 
with molecules that are in the electronic ground state which corresponds to 
A = 0. In the case of the asymmetric top where all moments of inertia are 
unequal such as the case for H 2 O there is no axis about which the angular 
momentum is conserved and the solutions cannot be described by spherical 
harmonics. The solutions still have 21+1 levels for each value of Z, but they 
are complex and there are no degeneracies [Herzberg, 1945; Goody, 1964]. 
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1.6 Rotation- Vibration Spectra 


The moment of inertia Jo is related to the bond lengths between atoms. 
When a molecule undergoes vibrational oscillation, the anharmonicity in the 
potential well forces the average equilibrium bond length for each vibrational 
state to increase. This increases the moment of inertia and decreases the 
energies of the rotational levels (1.81). The modified rotational energy levels 
are described by [Berry, Rice, Ross, 1980] 


Eyib- 


rot 


Eyib + -Erot 

Ey ib + J{J + 1 '}hcB„ 

J{J + l)hc - a e (u + ^) + • • • 


(1.82) 

(1.83) 


where B v is the rotational constant for the vibration level v, a e is a constant, 
and B e = h/ Air cIq is the constant which characterizes the energy of the 
rigid rotator. The energy levels (1.83) describe linear molecules such as CO, 
C0 3 and N 3 0. The energy levels are much more complicated for nonlinear 
molecules such as HjO and CH 4 and cannot be described by (1.82-3). 

If the energy of the incoming photon is equal to the difference between 
two energy levels of a molecule, a transition between the two states may 
be induced provided the molecule can interact with field as a whole. The 
probability of a transition between two states is proportional to |/x„ in . | 2 where 
/i nn . is the dipole between the two states n and n' 




(1.84) 
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The dipole operator /i is usually given by er, e being the charge of the elec- 
tron and r the position vector. There are additional terms to the transition 
probability involving quadrupole moments, magnetic dipole moments, and so 
on, however rates of these types of transitions are usually of an order of mag- 
nitude or more smaller than those for electric dipole transitions [Steinfeld, 
1985]. For the rigid rotator it can be shown [Goody, 1964; Berry, Rice, Ross, 
1980] that the probability for a dipole transition is non-zero if A J = ±1. For 
the harmonic oscillator the dipole selection rule is Ai/ = ±1. A molecule can 
undergo transitions in which both the vibrational and rotational quantum 
numbers change according to the selection rules. I will be concerned here 
with rotation-vibration transitions in the ground state since these are the 
dominant type of transitions observed in the ATMOS spectra. 

In (1.82-3) we see that the rotation-vibration energy is proportional to 
the constant he. We can define a new quantity which describes the rotational 
part of the energy 


c rot 


Eiot 

he 

J(J +1) ^i?e — a e (v + 2 ) + * * *j • 


(1.85) 

( 1 . 86 ) 


The unit of measurement of this quantity is cm -1 . The differences between 
energy levels Ae described in the manner of (1.85-6) is in terms of wavenum- 
bers ? which is the number of wavelengths of radiation within a centimeter 
of length in a vacuum. If A is the wavelength in vacuum then the frequency 
/ is related to it by 


/ = 


c 

A' 


(1.87) 
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For a transition of energy 

difference A E the corresponding frequency of ra- 

diation is 

A E 

(1.88) 


f = —j— = c Ac = c?. 

Explicitly, 

f <xq = 1/A. 

(1.89) 


Therefore frequency and wavenumbers are proportional to each other. 

In an absorption or emission process the energy of each photon involved 
must equal the difference between the energy levels. For linear molecules, 
if the lower energy level of the vibrational mode is given by the quantum 
numbers x/, J and that of the upper energy level by v',J' then the energy 
difference, in wavenumbers, is 

< = (*'- v)^ + J'( J' + 1 )Bw - J(J + 1 )B V (1.90) 


where 


= 


he 


(1.91) 


for the vibrational mode. For absorption processes in which A v = +1 there 
are two sets of rotation lines, those with A J = +1 and those with A J — —1. 
A J = 0 is prohibited in dipole transitions in the ground electronic state since 
there is no net dipole moment between the states. The set of transitions 
where A J = +1 is called the R branch and is given by J' = J + 1. The line 
positions for this branch can be written in the form 


Sr — Se + (*7 + l)(5 v +i + B v ) + ( J + lY(B v+ i — B„). 


(1.92) 
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The case where A J = — 1 is called the P branch and is given by J 7 = J — 1. 
The line positions are 

Cp = Se — + B v ) + J 2 (B v +\ — B v ). (1.93) 

The quadratic terms involve the differences due to vibrational stretching be- 
tween the rotational constants. In an infrared transition where the electronic 
states remain the same for reasons explained earlier the rotational constant 
decreases with increasing the vibrational quantum number v . From (1-92) it 
can be seen that a value for J is eventually reached where the frequencies be- 
gin to decrease with increasing J. This causes a band head - such as the solar 
CO band heads observed in the ATMOS spectra §11.3 - and occurs in the 
R branch in rotation-vibration spectra since B v+ x < B„ [Berry, Rice, Ross, 
1980]. The B v values, however, are usually very close to one another and it 
is difficult to observe band heads in spectra taken at ordinary temperatures. 
Typical values for the rotational constants is of the order of several cm -1 for 
atmospheric gases such as COj and N 2 O. The lines are spaced far enough 
apart that lines of the same vibration band do not interfere with each other. 
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1.7 Radiative Transfer 


The energy of a molecule in the ground electronic state consists of trans- 
lational, vibrational and rotational energy. The vibrational and rotational 
energy are associated with internal degrees of freedom while the transla- 
tional motion is related to the center of mass coordinates. The molecules 
can transfer energy among themselves via collisional and radiative processes. 
Molecules can lose or gain translational energy due to elastic collisions with 
other molecules. The internal and translational energy can also be exchanged 
in inelastic collisions. When a gas is subjected to radiation of frequencies 
comparable to differences in energy levels of its molecules energy can be ab- 
sorbed and emitted. The .populations of the various energy levels is therefore 
related to the balance between collisional and radiative processes. 

In the lower atmosphere below 70 km the collision rates are larger than the 
radiative rates due to the higher molecular density [Gao, 1988]. The popu- 
lations of the the rotational- vibrational levels are then dominated by molec- 
ular collisions and are described by the Boltzmann distribution [Thorne, 
1988], A quantity called the temperature T can be defined which specifies 
the mean kinetic energy | kT of the colliding particles and the particles follow 
a Maxwellian velocity distribution. This is called local thermodynamic equi- 
librium (LTE) since the molecules cannot be considered to be in complete 
thermodynamic equilibrium due to radiative processes between themselves, 
the rest of the atmosphere, space and solar radiation. The molecules can be 
considered to be in local thermodynamic equilibrium with the total radiation 
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of the atmosphere which is described by the Planck blackbody distribution 
at the kinetic temperature temperature T [Thorne, 1988]. In the following 
derivations we assume that the condition of LTE is adequate and this was 
found to be approximately true by Gao (1988) for the i/ 3 rotation-vibration 
transitions of COj at altitudes between 70 and 130 km. 

SPONTANEOUS EMISSION 


Spontaneous emission at a frequency / 0 occurs naturally from an excited 
state a to a lower energy state s' . It is induced by zero point oscillations in 
the radiation field [Cohen-Tannoudji, 1977]. The rate of this process 
is related to the electric dipole matrix element = (s|fz|a') between the 
states by 




16t r 3 / 0 3 
3e 0 he 3 




(1.94) 


[Thorne, 1988]. The mean lifetime between the states is then given by 


T..i = A 


- A - 1 


(1.95) 


STIMULATED EMISSION AND ABSORPTION 

Besides the spontaneous emission of radiation from an excited state to 
a lower energy state there can also induced transitions in the presence of 
external fields. Suppose a lower energy state has the energy hf,i and the 
upper state has the energy hf t . An electromagnetic wave with frequency 
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fo = ft — f t > can be described by an electric field 

£ = £o cos(27r f 0 t) = £ 0 cos(2tt f g t — 2irf t >t) 

= So [cos(27r f t d) cos(27r/ # t) + sin(27r/,i£) sin(27r/,£)] . (1.96) 

Each element of the sum (1.96) has a component that oscillates with the 

frequency of each state. It is these resonances with the dipole moment of 
each state that eventually induces transitions between the two states [Aller, 
1963]. The photon of frequency f 0 can drive a transition between the two 
states, with an efficiency which depends on the strength of the dipole between 
the states. The electric field can induce both absorption and emission of 
a photon and the rates of both of these processes are proportional to the 
average photon intensity I near the transition frequency v 0 [Thorne, 1988] 

7»b* = (1.97) 

Tenu = (1.98) 

respectively. The coefficients At-*,', £,*_♦, and are commonly referred 

to as Einstein coefficients and are used to describe the rates of the various 
emission and absorption processes. 

RADIATIVE BALANCE 

If n t > and n, are the density of molecules in the lower and upper state 
respectively then the relationship between the absorption and emission pro- 
cesses in radiative equilibrium is given by the relation 


= n,A t -> t i + n, ■£?,_,» J 


(1.99) 
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since in equilibrium the rate of photons emitted is equal to the rate of photons 
being absorbed. Then the photon intensity can be solved for 


7 = 


/B„. 

(iiji / ji, Si-tj' ) 1 


( 1 . 100 ) 


The ratio of the populations n' t /n t is undetermined. In local thermodynamic 
equilibrium at temperature T the ratio of the populations in each state is 
given by Boltzmann’s equation [Steinfeld, 1985; Thorne, 1988] 


n »* _ 9>< hfa/kT 

n, g , 

where are the degeneracies of each state. Then 

j A. t _ t > /■£?,_,» 

However from Planck’s Law [Aller, 1963] 


2hf 3 / c 3 

* e hf/kT _ l 


( 1 . 101 ) 


( 1 . 102 ) 


(1.103) 


for a blackbody radiator at the temperature T. Since If varies slowly over 
the width Ai/ of the line, 7 can be replaced with 7/o 
2hf 3 /c* 


e hfo/kT _ ! (g t ,B t ,^ l /g.B t -. t .) e h J°/ kT - 1 ‘ 


(1.104) 


Since this must hold for all T it can be shown that the Einstein coefficients 
are then related to each other by 


A, 


•' = 



(1.105) 




(1.106) 
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and depend only on the electric dipole matrix element and the energy differ- 
ence between the states. 

The case where LTE is not valid is discussed by Gao (1988) and Thorne 
(1988). At altitudes in the atmosphere high enough where molecular colli- 
sions are less frequent than radiative processes, the populations of rotational- 
vibrational energy levels are dominated by absorption and emission processes. 
The kinetic temperature T is then not adequate to describe the populations 
of the states. The lower atmosphere is the main source of photons present 
for absorption and the blackbody radiation distribution which describes them 
may not correspond to the local kinetic temperature of the gas [Gao, 1988]. 
The atmosphere is then considered to be in a state of non-LTE (NLTE). 
The grand effect of NLTE is to alter the populations of the various occupied 
states. It was found by Gao (1988) that between 70 and 130 km the devia- 
tions of the atmosphere from LTE do not affect the populations of rotational 
levels significantly from the Boltzmann distribution at a kinetic tempera- 
ture corresponding to the temperature of the lower atmosphere. Although 
LTE/NLTE is important for determining the populations of various states, 
it is not important for obtaining wind speeds §§II. 


1.8 Line Shapes 


The shapes of spectral lines depends on the molecules and their physical 
conditions. In this section we examine some of the most important physical 
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processes which describe spectral line shapes. Due to the uncertainty princi- 
ple the finite life times of excited molecular states causes the lines to have a 
natural width (in photon energy) and shape. This width is a characteristic 
fundamental to the molecule and depends only on the frequency of the tran- 
sition and the dipole matrix elements between the states (1.94). The thermal 
motions of the molecules also causes the frequencies of radiation involved in 
the absorption/emission processes to be Doppler shifted from the natural line 
shape as described by (1.1). The distribution of velocities and hence Doppler 
shifts are governed by the Maxwell velocity distribution. Finally, a molecule 
may also undergo collisions with other molecules in the gas. A molecule can 
interact with the radiation field during the period of time between collisions. 
A sum of all frequencies over that period of time can appear to the molecule 
as the correct driving frequency due to the properties of the Fourier trans- 
form. Frequencies well away from the transition frequency may drive the 
transition. All of these mechanisms as well as combinations of them can be 
described by a line profile ^(c) which describes the normalized distribution 
of the photon energies which induce transitions. 

The width A E of an energy level s is implied by the uncertainty relation 

A E t T.>± (1.107) 

where r. is the mean lifetime of the state. If the lower state has an infinite 
lifetime, which is the case for transitions to the ground state, then width of 
the transition, in wavenumbers (1.85), is 

1 

47TC T, ’ 


Ac. > 


(1.108) 
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In general, each transition has a natural lifetime given by 

r,- 1 = (1.109) 

where s' denotes all lower energy states and is given by (1.94), the 

coefficient of the wave function of the state s is then of the form 


V> oc r t 1 / 3 e“ t / 2r * 


( 1 . 110 ) 


since the square modulus of the wave function represents the probability of 
finding the state and r~ 1 e _t / T * is the probability density of a decay between 
times t and t + dt. Consider a transition at frequency / 0 . The decay causes 
the emission of photons about the frequency /<,. The electric field £ due to 
this transition caused by emission decays as [Thorne, 1988] 

£(t ) oc (1.111) 


where w 0 = 27 r/ 0 = 2i rc? 0 . The strength of the electric field decays with 
increasing time and has an oscillating part due to the transition frequency. 
This electric field has a Fourier transform 

£{u) = (27T)- 1 dt £{t)e iwt . (1.112) 

J — oo 

Since £(t) is only defined for t > 0, 

£(u>) <x y i(v — w 0 )J . (1.113) 

The line intensity is proportional to the square of the electric field [Jackson, 
1975] I u = [£(u>)| 2 and it can be shown by integration that 

jf \ j I 

0 *(?-*>)’ + (A<k) 2 * 


(1.114) 
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where Acn = (47rcr # ) -1 is the natural line width in wavenumbers which is 
in agreement with the width predicted by the uncertainty principle (1.108). 
The total line intensity integrated over all wavenumbers is defined as 

Jo=r°° /(?)*. (1.115) 

J — OO 

The line profile for natural broadening is 


<£n(?) = — 

W J ft) 2 + (A 0* 

and this is normalized as <^jv(c) d<; = 1. 

As an example, for an infrared transition to the ground state t, = Aj_i. 

Since f 0 = eft then (1.94) becomes 



A3-.1 


16x 3 gp 

3e 0 h 




(1.117) 


For a typical i/ 3 C0 2 transition we can choose, for exampleoul-m, ?o = 2330 
cm -1 and /ijj ~ 10 -29 Cm aa typical values. Then 1 ~ 3.6 x 10 4 Hz. 
This gives Aftj ~ 7 x 10 -9 cm _1 as the natural width which indicates that 
these lines are very narrow, about six to seven orders of magnitude smaller 
than the resolution of the ATMOS interferometer. 

The thermal motions of the molecules in a gas cause the frequencies of 
the absorbed or emitted photons to be Doppler shifted 


Co ~ go v 
ft c 


(1.118) 


where is the wavenumber in the observer’s frame, ft is the frequency in the 
molecule’s rest frame and v is the relative speed between the molecule and 
the observer along the line of sight. The fraction of molecules moving with 
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speed v and v + dv along the line of sight is given by a normalized Maxwell 
distribution [Thorne, 1988] 


dn _ fakT 
TO V TTTTO 


2 kT 

21 

mv 4 

dv 

CJLp 

7TT71 

2 kT 

l2kT 

_ [ mc J (c - <ro) 3 l 

V rm ^ l 

2kTq> J 


dq. 


From this result the Doppler line profile is given by 


M?) = 


A?D\/5r 


exp 


- (— V 


(1.119) 

( 1 . 120 ) 

( 1 . 121 ) 


The Doppler width is 

A® = C „ (1.122) 

¥ me 2 

As an example, for a temperature of 300 K, Co = 2300 cm -1 , and a COj mass 
of 44 gm/mole, A® = 2.6 x 10 -3 cm -1 . Note that the Doppler broadening is 
much larger than natural broadening 


A?*) A qif. 


(1.123) 


in the infrared regions but is still about 5 times smaller than the resolution 
of the ATMOS interferometer. 

A molecule can also have collisions with other molecules. Let the mean 
time between collisions be re* For simplicity consider a gas with a density 
such that the time between collisions is much greater that the interaction 
time rt, the time for the collision processes. If the time between collisions 
is <* the radiation affecting the molecule will subject to random phase inter- 
ruptions caused by collisions [Rybicki, Lightman, 1979]. Specifically, between 
collisions, the component of the electric field which affects the molecule is 

£'(t) oc e~ iu,ot (1.124) 
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for 0 < t < t* where again u> 0 = 2vf 0 = 2ircso- By using (1.112) the Fourier 
transform of this field is 


Sm |L M _ 1 T 0 )/2 /2] ' 

Again since I'(w,t*) = |£'(w,i*)| 2 then 

'<-■>» ptfSE 3 )’- 


(1.125) 


(1.126) 


The probability dV{t*) that the molecule has a collision between time t* and 
t* + dt* is 

dV{f) = Tq 1 e~ l * l Tc dt* . (1.127) 

The final intensity distribution is given by the average of the individual in- 
tensity distributions (1.126) over the probability function (1.127) 


/( w) = re 1 [°° I' {u>,t*)e- r/Tc dt*. 
Jo 


(1.128) 


This integral can integrated by standard methods to be 

= Io *(t-*>) 2 + (A ?c ) 2 ^' 129 ^ 

where the collisional width, also called the Lorentz width, is given by 


A?c = (27TCT C ) 


-l 


(1.130) 


and J 0 is the total line intensity (1.115). The profile for pressure broadening 
is also a Lorentzian like natural broadening 

1_ A?c 

ir (? - To) 2 + ( A?c) 2 ’ 



(1.131) 
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The mean time between collisions can be estimated by 

1 1 3k T 

rc = (mean free path)/(thermal speed) = — — / \l 

nli V m 


(1.132) 


where n is the density of molecules, S is an average collision cross section, 
m is the molecular mass. From (1.130) 


Acc = ~PZ] 
c 


1 — r— 77= (1-133) 

47r 2 mfcT 

by using the ideal gas equation P = nkT. Since atmospheric densities de- 
crease with height exponentially we expect that collision broadening to be 
important primarily near the surface. As height increases the collisional 
frequencies decrease in a like manner and we expect that the lines axe gov- 
erned by the Doppler shape. Above about 50 km the primary line shapes 
are Doppler with small Lorentz contributions. Below about 20 km pressure 
broadening is important. 

In the regions between 20 and 60 km altitudes the atmospheric absorp- 
tion lines can show both Doppler and Lorentz (pressure) broadening. The 
observed line shape is actually a set of Lorentz profiles averaged over the 
Maxwell velocity distribution [Rybicki, Lightman, 1979]. The line shape is a 
convolution of the Doppler and Lorentz profiles and is called the Voigt profile 

<M?) = <fc>(0 * <£l(0 = / dq’ fo(q' - q) <f>c(s')- (1.134) 

t/*~ oo 

This profile is normalized since the Doppler and Lorentz profiles are normal- 
ized and this can be seen by integration 

f dq = / dq' <£l(c') / dq' <j>{q' - ?) 

J — OO J — OO * — oo 


1 . 


(1.135) 
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The Voigt line shape is governed by two widths, the Doppler width Acd and 
Lorentz width Acc- When pressure broadening becomes dominant Acc 
A cd and the line shape becomes a pure Lorentzian. In the upper atmo- 
sphere where pressure broadening is negligible compared to Doppler broad- 
ening Acd Acl the line center has a Doppler shape while the ‘wings’ - the 
absorption outside of the Doppler width - have the Lorentz shape. 
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1.9 Absorption and Emission Coefficients 


The ATMOS spectra show the absorption of the atmospheric gases along 
the line of sight between the observer and the sun. The ray path passes 
through all altitudes above the tangent point of observation. In this section 
we discuss how the transmitted intensity depends on the amount of absorber 
along the line of sight. We use this information to examine how important 
the contributions of the absorber in regions far away from the tangent point 
contribute the the observed line positions in §11.13. 

The line shapes discussed in the previous section apply to large ensembles 
of molecules undergoing radiative and collisional transitions. Large amounts 
of gas along the line of sight between the radiation source and the instrument 
enable the molecules to absorb significant amounts of the incoming radiation 
due to cumulative effects of absorption. The absorption coefficient is 
defined to describe the loss of beam intensity dl f at wavenumber q as it 
travels over a distance ds of the medium by 

dl t = —ci'If ds. (1.136) 

A spontaneous emission coefficient j t is also defined as the power emitted per 
unit volume. This can be assumed to be isotropically distributed since there 
is no preferred direction of radiation in an emission. If an incident beam 
with parallel rays such as solar rays in the ATMOS experiment is used, then 
the spontaneous emission term j t can be neglected since most of the emitted 
radiation goes out of the field of view. 
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Let a be the upper energy state of a transition, a' the lower energy state 
and n the number of active molecules per unit volume. Then he q n,A,- tt i<f>(q) 
is the power emitted by spontaneous emission per unit volume per unit 
wavenumber where hcq is the photon energy and n,A,-, t i<f)(q) is the rate 
of photons emitted per unit wavenumber per unit volume. This radiation 
caused by this term is isotropically distributed and can be neglected relative 
to the beam absorption. 

The power absorbed per unit wavenumber per unit volume is the differ- 
ence between absorption and stimulated emission since the stimulated emis- 
sion will be polarized with respect to the incoming radiation [Steinfeld, 1976] 
and can be thought of as “negative absorption” 

he q (n,» ^(c). 

Since the intensity / c is the power per unit wavenumber per unit area the 
power absorbed per unit wavenumber per unit volume is the change in in- 
tensity per unit length along the beam 

= hcq ( (1.137) 
as 

Then the absorption coefficient can be determined from (1.136) 

= hcq — 7v2? f «_ # )<£(c). (1.138) 

If the system is in thermodynamic equilibrium with itself, which is the case 
with of an atmosphere in LTE and will be assumed to be approximately true 
for NLTE [Gao, 1988] then from (1.101) 

n »' _ 9 •' hct/kT 

n. g. 


(1.139) 
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By using this in (1.138) and using the Einstein relation (1.106) 

= hcqn.,B.,^.{ 1 - e~ hc ^ kT )<l>( c). (1.140) 

The transmitted intensity neglecting the spontaneous emission term is 


the solution of (1.136) 

dI < „ T 

da ~ 

(1.141) 

which can be solved 

I ; = l f (0)e~ Jo a<dt 

(1.142) 


where L is the total path length the radiation travels through the absorber 
and 7<(0) is the radiation intensity in the absence of any absorption. 

It was shown §1.8 that the line widths are much smaller than the central 
wavenumber <; 0 . Also, at infrared wavenumbers and at typical atmospheric 
temperatures the exponential factor in (1.140) is much less than unity and 
the absorption coefficient can be well approximated by 


Qf — TlgiltC B 

Then the transmitted intensity (1.142) can be written as 

h = A(0)e- fc '“. 


(1.143) 


(1.144) 


The number of molecules per unit area in the lower state is defined by 

rL 


u = f da n,i 
Jo 


(1.145) 


where it is should be recalled that n,< is the density of molecules in the lower 
state. The absorption cross section per molecule is given by 


k t = 


(1.146) 


and the line strength is defined as 


The line strength 
the line profile is 


S * — he 


(1.147) 


is also equal to the integrated absorption coefficient since 
normalized and this can be shown by integrating (1.146) 
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1.10 Optical Depth, Equivalent Widths, and Curve 

of Growth 


The optical depth r< at a given wavenumber c is related to the mean free 
path of a photon [Thorne, 1988]. We define 

r c = k K u (1.149) 


so that (1.144) is 

= 1,(0) e-». (1.150) 

For 1 the medium is “optically thin” while for r t 1 it is “optically 

thick”. For the case of optically thin layers we can expand (1.150) 

= J,(0)(l-r, + •••). (1.151) 


In this case the transmitted intensity I f decreases linearly with an increase in 
the amount of absorber u. As the amount of absorber u increases this linear 
behavior no longer holds true and saturation is reached. The absorption line 
initially saturates at the line center and then outwards towards the “wings”. 
Figure 1.10.1 shows this saturation as the amount of absorber increases. 

The total loss of incident radiation due to absorption by the line is 


/ d<; [I*(0) — /*(tv)] = Power Loss/Area. (1.152) 

J line 

The ratio of this quantity to the incident intensity per unit wavenumber 
defines the equivalent width 

m * 



(1.153) 
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From (1.150) we have 



(1.154) 


The equivalent width is independent of the instrumental resolution [Thorne, 
1988]. Although a decrease in resolution will cause the observed line to 
“spread out” over a greater spectral region, the equivalent width remains 
the same as the ‘true’ line. Thus measurements of the equivalent widths 
of absorption lines can be used to obtain information regarding the total 
absorption of the line and properties associated with it such as temperature, 
pressure, and absorption cross sections u [Gao, 1988]. If the instrumental 
resolution becomes larger and hence the resolved width of the line at, say, 
half maximum absorption increases, the peak absorptances must decrease 
since the equivalent width with remains the same. 

A typical characteristic of an absorption feature is its width Ac as de- 
scribed in §1.9. This is characteristic of Doppler and Lorentz lines, as well as 
the convolution of the two which is called a Voigt line. If we change variables 


x = 


C - go 

Ac 


(1.154) can be written as 


W 

Ac 



For the Doppler line (1.121) 

T x = Su <£d(c) = * e - **. 

Acd -A 

For the Lorentzian (1.131) 


r m 


= Su<f> l(c) = 


Su 1 

7T ACC 1 + X 3 ’ 


(1.155) 


(1.156) 


(1.157) 


(1.158) 
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We see that in both cases the fraction of the equivalent width (EW) to the 
line width W/ Ac is a function of Su/ Ac. If we define 

Su 


a = 


Ac 


then (1.157-8) can be written 


^(Doppler) = —z= e" 
v' 7r 


T*(Lorentz) = 


a. 


a. 


(1.159) 

(1.160) 
(1.161) 


7T 1 + X 1 

For the two cases I computed the integrals (1.150) as a function of a and 
are plotted in Figure 1.10.2. Both curves have the same slope in the region 
where a < 1. This linear region is predicted by (1.156). For small r c we have 

W 

^ 

Ac Jiint 

For the Doppler and Lorentz lines we have 


dx Tm. 


(1.162) 


(W\ 


\Av) 

j 

Doppler 

/ W' 

\ 

\Ai/, 

' Lorenti 


~ a 


1 , i 

~ a— dx = a. 

7T J—oo 1 + X 2 


(1.163) 

(1.164) 


The growth rate is the same for both type of lines in the linear regions. The 
equivalent width is then directly proportional to the amount of absorber. The 
linear growth for both Doppler and Lorentz lines continue until log Su/ Ac ~ 
1/2 which can be seen in Figure 1.10.2. 

As the amount of absorber continues to increase the equivalent width of 
a Doppler line Wd approaches almost ten times the line width A?t> and in- 
creasing the amount of absorber u does not increase the equivalent width sig- 
nificantly beyond this which can be seen in Figure 1.10.2. Since the Doppler 


59 


width A?d depends primarily on the temperature T and not the pressure the 
equivalent width W b will not change significantly, with height. 

For a pressure broadened line, the curve in Figure 1.10.2 shows that the 
equivalent width of a Lorentz line Wq increases indefinitely relative to the line 
width A?c- Since the line width depends directly on the pressure P (1.133), 
as the altitude decreases and the amount of absorber along the line of sight 
remains significant, the increasing pressure causes the equivalent width Wq 
of the lines to grow without limit. In the case of strong absorption lines such 
as the 1/3 COj band, large spectral regions can become opaque due to strong 
pressure broadened saturation. 

The increasing equivalent width limits the altitude ranges where pressure 
broadened lines may be useful for measuring wind speeds since the lines as 
observed by the instrument must not be distorted §1.12. If the equivalent 
width of the line is much smaller than the instrumental resolution, the line 
will yield approximately the same shape as an infinitesimally thin absorption 
line with the same equivalent width. “Flattening” of the line shape near the 
line center due to saturation can not be seen in the spectra and does not 
significantly effect the resolved line shape provided that Acd resolution. 
If the equivalent width becomes significant with respect to the instrumental 
resolution we expect that spectral features will become broadened. As shown 
in the next section, increasing the width of a feature in the processed spectra 
can significantly degrade the precision of determining line positions. Since 
the equivalent width of a Doppler line Wq is limited to about ten times 
the Doppler width A?d, and since the Doppler width is about 100 times 
smaller than the instrumental resolution, Doppler lines can be used over 
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altitude ranges spanning tens of kilometers for obtaining wind speeds without 
having to be concerned with significant increases in the Doppler equivalent 
width. On the other hand, the equivalent width of pressure broadened lines 
increases rapidly with increasing pressure and absorber amount u. Pressure 
broadened lines therefore have smaller ranges over which they can be useful 
for measuring Doppler shifts and lines used for measuring relative speeds 
need to be carefully chosen to avoid the effects of pressure broadening in the 
resolved spectra. 


Figure 1.10.1 A sketch of line saturation as the optical depth changes from 
‘thin’ (r < 1) to ‘thick’. At saturation the line becomes opaque near the 
line center ? 0 . 

Curve of Growth 



Figure 1.10.2 Curve of growth for Doppler and Lorentz lines computed from 
(1.156). The equivalent width Wp of the Doppler lines becomes almost 10 
times its Doppler width Acd as the amount of absorber v increases. The 
equivalent width of the Lorenz lines Wq increases with u without limit rela- 
tive to its collision broadened width A?c- 
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1.11 Measurement Precision Theory 


The purpose of this section is to examine how the measurement precision 
of the wavenumber of an absorption line is affected by the resolved shape 
of the line with a noisy background and digital frequency spacing. In the 
ATMOS spectra the noise level is about a(N ) = 1% of the background 
signal and the wavenumber spacing is 0.0075 cm -1 . For an absorption line 
with a typical resolved width of .015 cm -1 at half maximum, there are only 
3 points to use to estimate the line position. The finite resolution of the 
spectra causes the lines to be “spread out” and the data used to obtain 
line positions correspond to the true spectra convolved with the (apodized) 
instrument function. Therefore the ‘line shapes’ present in the spectra are 
degraded by the instrumental resolution. In this section the ‘width’ of a line 
is the apparent or geometric width as it appears in the processed spectra. 
This instrumental width is limited by the resolution of the spectra. 

Provided pressure saturation is not present, we can approximate the ge- 
ometric shape of the a line in the convolved spectra by a parabola of the 
form 

9(x) = (1.165) 

for —w/2 < x < w/2 where x = 0 corresponds to the line center, A p is the 
geometric peak absorptance of the line and w the geometric width at half 
maximum. Figure 1.11.1 shows an example of this function. In order to 
test our ideas, we allow this function to describe the shape of the absorption 
line as it would be measured in the absence of noise. We can model the 
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actual observed data with this “clean” signal and a random distribution 
of noise added to each point. A random noise distribution N(x) with a 
standard uncertainty <r(N) is added to this “noiseless” signal ^(x) to produce 
an “observed signal” < 7 (x) + N(x). 

Having obtained the data for this line suppose, for simplicity, that a trial 
function of the form 

f(x) = (2 A p /w 3 )(x - yf (1.166) 


is to be fit to the observed signal g(x) + JV(x) by the least squares method. 
We wish to determine the best fit and hence the line position from this fit. 
We sample the signal flr(x) + N(x) at the n points Xj. The measured line 
center y 0 is the value of y which minimizes the least squares sum 


M{V) ~ ^ 2o*{N) • 


(1.167) 


i= i 


The least squares sum M can be approximated by an integration over the 
range of sampled points 


M 


*• [ /(I) “ + 

2 J ^ dxN(x)(f(x) - g(x)) + J dxN 2 (x ) . 


(1.168) 


The second term on the right hand side is essentially zero since N(x) is a 
random function with a mean of zero. The last term can be ignored since it 
is a constant which only depends on the noise distribution <r{N). This gives 

M ~ 2^(N) l dx Wx) - S(I))I • (1 ' 169) 
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The value y 0 is chosen which minimizes M by the condition 

dM I 


dy 


= 0 . 


(1.170) 


1M 


-1 


At the stationary point yo the precision is given from elementary statistics 
[Bevington, 1969] by 

For our model functions 


(1.171) 


VO 


f(x) g{x) = ^ 3 P ( 2 xy + y 3 ) 

(1.172) 

»>r 2 nAp f 4 1 

M= „\N)v,> {y + i wy) - 

(1.173) 


so that 


The minimum occurs at y 0 = 0 which gives 

3 a(N)w 


<r(y 0 ) = 


4n 


(1.174) 


This relates the measurement precision to the geometric line width w, the 
noise-to-signal ratio <r(N), the number of points sampled n, and the geometric 
peak absorptance of the line A p . 

The result for <r(yo) which comes from the integral approximation can be 
tested against a more direct calculation to estimate the uncertainties from 
the least squares fit. Once again introduce the least squares sum 

m = £; M&l ~ (1.175) 

2<r 3 (lV) 

where this is related to the previous example by n = 2 m + 1 and Xj — 
(j /m)(w/2) for — m < j < m. Using the same functional forms /(*) = 
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(' 2A p /w 7 )(x — y ) 2 and g{x) = (2 A p /w*)x 2 , the least squares sum can be 
expressed, after discarding terms involving the random function N(x) for 
the same reasons given above, as 


M = 


/2 A b 




(«i - vf - (^*J + "My)' (1.176) 


2 2 




i 53 { m 2 y 3 ~ 2 j mw y + jw 3 ) (1.177) 


m 3 <r 2 (N)w* 

The sum (1.177) can be shown to be 

^ _ (12m 2 + 6 m)y 4 + (4m 2 -f 6m -f 2)u; 2 y 2 A£ 

3m a 2 (N)w*‘ 


(1.178) 


Once again the measurement uncertainty can be estimated from (1.165) with 
yo = 0 which gives 


r (yo) = 


I 3m <r(N)w 
8m 3 + 12m + 4 A~ 


(1.179) 


The result of this calculation shows the same dependence on the parame- 
ters w, and A p as the result from the integral approximation except for 

a multiplicative form factor f n which only depends on the number of points 


n 


<r(Vo) = fn <r ^ W . (1.180) 

Table 1.11.1 shows how these form factors are related for various values of m. 
The values listed show that the integral approximation is higher by as much 
as 20% for small n than the direct sum calculation with better agreement as 
n increases. 

This method of determining line positions is similar to the procedure 
which is followed by the line-finder program [Norton, 1986] which determines 
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the line positions from the ATMOS spectra. The model for the absorp- 
tion line shapes is somewhat simplified and it is of interest to see how this 
parabola-parabola least squares fit compares to a fit of a parabola to a gaus- 
sian absorption. This pair of functions is chosen since the true absorption 
lines ‘taper off’ at wavenumbers away from the line center similar to a gaus- 
sian shape. Figure 1.11.2 shows such a sample gaussian described by 

g{x) = A p { 1 - c - 41n2 ( x /“) a ). (1.181) 

Again A p represents the geometric peak absorptance and w is the geometric 
width at half maximum. 

Although the gaussian line shape is chosen rather arbitrarily the depen- 
dence of the measurement precision to the geometric width w , the peak 
absorptance A p and the noise amplitude <r(N) remains the same as for the 
parabola-parabola fit. To see this, we fit 

f(x) = (2Ap/w a )(x - y) a (1.182) 

to g{x) for —w/2 < x < w/2 by both the integration approximation and 
the direct sum methods. To solve for the measurement uncertainty we again 
followed the procedure above (1.169-71). The rather involved calculations 
were performed using the symbolic manipulation program MACSYMA at the 
VAX-8600 facility in the Physics Department at the Ohio State University. 
It was found that, by closely following the procedure (1.169-71), the result 
using the integration approximation is 

<r(N)w 
A p 


« 
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(1.183) 


Once again the result is a form factor multiplied by the product <r( N)wA~ x . 
For the direct sum calculation the resulting measurement uncertainty fol- 
lowing a procedure identical to (1.175-79), also evaluated by MACSYMA, is 
given by 



The form factors are shown in Table 1.11.2. Again the integration approx- 
imation introduces larger uncertainties than the direct sum method, how- 
ever, even more interesting is that the values for the parabola-parabola and 
parabola-gaussian direct sum methods give approximately the same answer. 
The former case has slightly smaller form factors. These results suggests 
that the measurement precision may not be very sensitive to the details of 
the geometric line shapes which can be seen by comparison of the direct sum 
methods which are more accurate. 


The calculations above suggest that the measurement precision is pro- 
portional to the noise level <r(N), the geometric line width W , and inversely 
proportional to the peak absorptance A p . The signal to noise ratio (SNR) 
in the ATMOS spectra is approximately 100 in the 1/3 CO* region near 2350 
cm -1 which corresponds to c(N) = 0.01. For the apodized spectra the width 
at half maximum peak absorptance is about 0.015 cm -1 , which can be at- 
tributed to the instrumental resolution and can be seen from a few sample 
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line shown in Figures 1.11.3-4. (The ATMOS spectra is digitally processed 
at intervals of 0.0075 cm -1 and due to a “smoothing” procedure written by 
Norton the drawings of the absorption lines appear as curves.) The width of 
the line is approximately twice the digital spacing which gives n = 3. Table 
1.11.3 shows estimated measurement uncertainties using these parameters 
for lines with various values of A p . 

From this information we can estimate the precisions by which Doppler 
shifts can be measured. From 1.2 the uncertainty in measuring relative speeds 
v can be determined 

<r(v) = <r(A?)/? • c. (1.186) 

The measurement error per line is dependent both on the precision by which 
line positions can be obtained er(?) and the central wavenumber c- The 
measurement precision <r(?) is frequency independent. Therefore, for lower 
wavenumber regions the error in measuring the relative speed a(v) increases. 
Using the data presented in Table 1.11.3 we see that, as an example, for lines 
with a peak absorptance of about 30 % and central wavenumber <; ~2000 
cm -1 the precision of measuring relative speeds is about 25 ms -1 per line. 
Therefore a set of 25 lines can bring the mean precision to 25 ms -1 /\/5 =5 
ms -1 . This measurement precision theory predicts that precisions of obtain- 
ing 5 ms -1 are within the capabilities of the interferometer. Using higher 
frequency regions and stronger lines it is feasible to improve these statistics. 
However, due to the broadening effects discussed previously which may limit 
the altitude ranges for using atmospheric lines to obtain Doppler shifts, a 
wind speed measurements with precisions of 5 ms -1 was determined to be 
within the capabilities of the ATMOS instrument. 
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Figure 1.11.3 A C0 3 absorption line taken from ATMOS spectra. 
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Figure 1.11.4 An H 3 0 absorption line taken from ATMOS spectra. 


Table 1.11.1 Form factors for parabola-parabola least squares fit for various 
values of m. 


m 

n = 2m + 1 

Integral 

Direct Sum 

1 

3 

0.50000 

0.35355 

2 

5 

0.38723 

0.31623 

3 

7 

0.32733 

0.28347 

4 

9 

0.28868 

0.25820 

5 

11 

0.26112 

0.23837 


Table 1.11.2 Form factors for parabola-gaussian least squares fit for various 
values of m. 


m 

n = 2m + 1 

Integral 

Direct sum 

1 

3 

0.62051 

0.35355 

2 

5 

0.48065 

0.32063 

3 

7 

0.40622 

0.28924 

4 

9 

0.35825 

0.26440 

5 

11 

0.32405 

0.24464 


Table 1.11.3 Estimated measurement uncertainty from parabola-gaussian 
least squares fit using direct sum. 


El 

<r(y o )(xl0 4 cm _1 ) 

El 

<r( 2 /o)(xl 0 -4 cm *) 

EH 

5.30 

ES 

1.06 

i 

2.65 

m 

0.88 

i 

1.77 

0.7 

0.76 

0.4 

1.33 

0.8 

0.66 
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CHAPTER II 


EXPERIMENTAL METHOD 


II. 1 Instrumental Calibration 


CALIBRATION MODEL 


To get useful wind speed measurements we must be able to measure 

Doppler shifts precisely for each spectrum and must be able to intercompare 

the spectra. The comparison of relative speeds determined from a series of 

spectra is limited by the calibration stability of the interferometer. As well 

as noise in the spectra, there is the possibility that the calibration contains 

random offsets due to errors in locating ZPD in each interferogram §1.1. 

An estimate of the magnitude of this scale uncertainty provides information 

about the precision of the ZPD locator program as well as setting limits of 

* 

the accuracy of relative speed measurements. Sets of spectra may also have 
slowly changing calibration offsets caused by variations in the environment 
of the interferometer as well as “man-made” calibration offsets introduced 
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during the processing stage [Norton, 1986]. In order to obtain accurate mea- 
surements of Doppler shifts it is necessary to calibrate the spectra accurately. 

Of the four spectral filters used with the interferometer only Filter 3, 
which covers the region 1600-3200 cm -1 , provided spectra with lines due to 
gases trapped or around the instrument by which the calibration stability 
could be analyzed. The strongest lines of the v 3 C0 2 and i/ 2 H 2 0 vibration 
bands from 2300 to 2380 cm -1 and 1600 to 2000 cm -1 respectively appeared 
in spectra taken at tangent heights far above the atmosphere. These lines 
appeared to be related to the environment of the instrument and were due to 
residual gas trapped in the instrument as well as gas floating in space with 
the orbiter. These lines were called instrumental lines and it was found that 
they provided a tool for directly measuring the behavior of the calibration 
system. Figures 2. 1.1-2 show 6ome of these absorption lines as they appeared 
in the ATMOS spectra. In Figure 2.1.1 the v 3 C0 2 band is shown as well as 
the v = 0 — f 1 band head of solar CO. In Figure 2.1.2 some strong lines of 
i/ 2 H 2 0 band are shown together with other weak lines of solar CO. 

In order to test the accuracy of the instrumental calibration the measured 
positions from a set of instrumental lines tj* of each line j in each spectrum 
k were assumed to be related to accepted standard positions by the model 

C* = *(1 + a) + e) + e* . (2.1) 

Here, a is an offset which is used to model any wavenumber errors in the 
calibration of the interferometer. This global offset is assumed to vary slowly 
over large wavenumber regions of the spectra and we will consider it be rela- 
tively constant over regions of several hundred cm -1 . The source of this kind 
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of offset is assumed to be related to environmental effects or arbitrarily intro- 
duced data processing artifacts [Norton, 1986]. The calibration noise due to 
phase errors in the processing stage is given by e*, a frequency independent 
scale error for the spectrum k , and ej, the measurement error of the line j 
in the spectrum k. It was shown in §1.1 that the calibration uncertainties 
due to phase errors should be negligible but they are still included here. It is 
assumed that all of these errors are independent of each other. The position 
of each spectral line j in a given spectrum k is measured by a line-finder 
program developed by Norton (1986) which estimates the positions of ab- 
sorption lines from the spectral data and its precision of C?(10 -4 cm -1 ) is 
discussed in §1.11. The measurement uncertainty <r(e,-) and scale uncertainty 
<r(e) both contribute to the position uncertainty and it is necessary to 
estimate the magnitudes of each of these. 

DRIFT PARAMETER 6 


Before proceeding along these lines it is important to check the calibration 
stability in each set of spectra for each occupation. A reference spectrum r 
was arbitrarily chosen, in this case at the highest tangent height, and a check 
for a calibration drift factor 6 was made. It is assumed that there exists the 
possibility that 

= (1 + «)<; ( 2 - 2 ) 

where k are the other spectra in the set and 6 may change between spectra in 
an occupation, a sort of calibration “drift” factor. Stability will be inferred 
if there are no statistically significant deviations away from 8 = 0. The 
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measured positions of each of the CO3 and HjO lines j in each spectrum k 
were compared to their corresponding positions in the reference spectrum r 
and a value for 8 * determined by 



An average value 

«* = fi), 


(2.3) 


(2.4) 


and the uncertainty <r(8 k ) for each spectrum was determined. The brackets 
stand for an average over the subscript j - i.e. all of the lines used in spectrum 
k. 


Estimated values of 8 k for the five occultations are shown in Figures 
2. 1.3-4 for CO2 and HjO instrumental lines respectively. The figures show 
the distribution of the drift factor of each spectrum as a function of tangent 
height. The values for 8 have a scatter of about 5xl0 -8 and about 8x10"* 
for the CO3 and H3O lines respectively. There does not appear to be any 
‘trend’ to the data and no apparent calibration drift. Figures 2. 1.5-6 show 
the respective uncertainties <r(8 k ). The uncertainties are of the order as the 
scatter of the 8. This indicates that the scatter of the drift factors 6 are due 
to noise in the data and are not caused by any real drift. Tables 2.1. 1-2 list 


* = <«*>*. 


(2.5) 


the average over all spectra k and the standard error <r(6) for CO3 and 
H3O lines respectively for the five occultations. From these results there is 
no evidence of a measureable drift of the frequency calibration during an 
occultation since the average values of 8 found are of the same order as 
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the measurement uncertainty in all of the spectra k. This indicates that 
the calibration does not “drift” between spectra during an occultation. The 
fractional change of a line position predicted in Table 1.11.3 is of the same 
order as these fluctuations. As an example, a line with a peak absorptance 
of 40 % at a frequency of 2330 cm -1 can be predicted to have a fractional 
measurement error of 1.33 x 10 -4 /2330 = 5x 10~ 8 . This corresponds to wind 
speed errors of about 15 ms -1 per line. 

GLOBAL CALIBRATION OFFSET a 


It is useful to compare the instrumental line positions to a set of accepted 
standards in order to estimate a for each occultation. Assume that such a 
set of standards Cj exist. Sets of COj and HjO standards recommended by 
Brown (1986) were used and are listed in Tables 2. 1.3-4 respectively. From 
(2.1) we can write 

a = 1 + error. (2.6) 

For each line j in each spectrum k a value 

= «?/& - 1 (2.t) 

is computed. The average value over all lines 

= <«;>, (2.8) 

was determined for each spectrum and is displayed in Figures 2.1.7-11 for all 
flve occultations using CO3 and H3O instrumental lines. The data show a 
variation from spectrum to spectrum of about 5 x 10 -8 . Notice that there 
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does seem to be a significant offset on the average between the values of a 
computed for the CO? and H 2 O lines. Figures 2.1.12-13 show o-(a fc ) for the 
COj and HjO lines respectively. The uncertainties for the CO 2 lines are 
about 2-3 xlO -8 which is of the order of the average fluctuations of a*. For 
H 2 O the values for <r(a fc ) are about 4xl0~ 8 and they are significantly larger 
than the C0 2 values but are also consistent with the average fluctuations of 
a k . This information indicates that the fluctuations of measured values of 
a k are due noise in the spectra rather than calibration fluctuations. The null 
result above for the calibration drift factor S also confirms this. 

Since a k was found to vary randomly due to instrumental noise we deter- 
mined the mean value for a in each occultation 


a = ( a k ) k . (2.9) 

The mean values a and mean standard error <r(a) are given in Table 2. 1.5-6 
for the five occultations as computed from the u 3 CO a and v 3 H a O bands 
respectively. Table 2.1.7 lists the differences between the two sets and it is 
evident, since the offset is larger than the mean errors, that the difference 
in their offsets Aa = aco a — £*H a o are statistically significant. Figure 2.1.14 
shows the differences A a k for all occultations. The general trend of the 
difference of the offsets Aa k appears to be between about 3-5 xlO -8 . The 
average difference between the instrumental global offsets a is +1.38 ±0.41 x 
10 -8 . This indicates that the discrepancy between the positions determined 
from the instrumental lines and the standards is significant. Since, according 
to Brown (1986), the CO 2 standards are “good” (accurate to 10 -4 cm -1 ) the 
H 2 O standard frequencies may be too low by this factor. This also suggests 
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that our assumption that the “instrumental” lines are at rest with respect to 
the interferometer may not be correct. This question of whether or not the 
“instrumental” lines are at rest is examined in §11.10. 

Comparing the differences in A a of the instrumental lines to the standard 
positions is like ‘comparing apples and oranges’. Either set of standards may 
be ‘correct’ but there is no way of determining it from the data. However, 
we do have insight into the magnitudes of the discrepancies. It is shown 
in §11.5 that any absolute calibration offsets do not affect the accuracy in 
which Doppler shifts can be determined, provided they are consistent from 
spectrum to spectrum and that a set of instrumental ‘absolute’ standards are 
available, such as the instrumental COj and HjO lines. 

The systematic relative velocity An shift caused by the global offset a 
for the two sets of instrumental lines is listed in Tables 2. 1.5-6. The velocity 
shift is given by Aw = a c where c = 299 792 458 ms -1 is the speed of 
light. The mean errors in the relative velocity shift o-(Aw) = c<r(a) are 
between 1 and 3 ms -1 . The mean values Aw are an order of magnitude or 
more larger than these uncertainties and therefore there were statistically 
significant calibration offsets between occultations. The sunset occultation 
SS09 was found to have a particularly large offset of about 80 ms -1 compared 
to the other occultations. These offsets were probably introduced during 
the processing of the spectra and are artificial in origin [Norton, 1986] - 
specifically, the calibration scales were multiplied by certain factors during 
testing of the spectra. It is evident that the wavenumber calibration of 
the spectra need corrections, to agree with calibration standards as well as 
consistency between occultations, for wind speeds to be measured. When 
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instrumental lines are present in the spectra this calibration can be achieved 
since the instrumental lines mark the rest frame positions. 

The frequency scales of all subsequent spectra were shifted by the correc- 
tion factor 1 — a to compensate for these offsets for each occultation. The 
positions of the instrumental u 3 COj lines were thus forced to agree with the 
standards. We can say that the spectra were calibrated with respect to the 
COj standards. The COj correction factors used are given in Table 2.1.5. As 
discussed above, the absolute calibration of the instrument is only as good 
as the 1/3 line standards themselves. We accepted these standards as the 
‘best available’ and calibrated the instrument relative to these. However, it 
is shown in §11.4 that the accuracy of the relative speeds is independent of 
the set of standards used provided that instrumental lines are available in 
the spectra for calibration. 

After these correction factors were applied to the frequency calibration 
it was then possible to make new estimates of the CO 2 and H 2 O standards 
from the instrumental lines. All of the instrumental lines from n = 125 
spectra from the five occultations were used. Tables 2 . 1.8-9 list the new 
positions from the ATMOS instrumental lines as compared to Brown’s values, 
A( = Qnu — (Brown for both regions. After calibration, QH a o was determined 
and the instrumental positions of the H 2 O lines differ from Brown’s standards 
by a factor a = — 1.25 dh 0.36 x 10 -8 which agrees with the estimate given 
previously. This can be seen by the multitude of minus signs in the differences 
column in Table 2.1.9. The average global offset for the CO a lines vanishes 
as required. 

Also listed in Tables 2. 1.8-9 are the precisions of the new positions. It 
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is worth noting that Brown’s standards are known to 10 -4 cm -1 accuracy 
while the precisions of the ATMOS values are almost an order of magnitude 
smaller than this. This suggests that there is an improvement of a factor of 
almost five in the line position precisions, particularly in the 1/3 CO 3 lines. 
However, most of the differences between the instrumental lines and Brown’s 
standards are within 10 ~ 4 cm -1 . Hence, this agreement suggests that the 
ATMOS interferometer can be used for determining precise relative positions 
for other lines based on Brown’s 1/3 C0 2 values. 

From this analysis we found that the instrumental frequency calibration 
remains stable during the course of each occultation. Provided that scale 
errors do not exist, and we show in the next section that this appears to be 
true, we can rely on the frequency calibration of the interferometer to mea- 
sure Doppler shifts and positions of lines relative to the overall calibration 
of the instrument. Measurement of the fluctuations in a k and fi k also hint 
that the measurement precision theory may be reliable and a more quantita- 
tive treatment is given in §11.11. After calibrating the spectra is was found 
t&at the measured positions of the instrumental CO 3 and H 3 O lines were 
significantly improved over the precisions of the standards. Typical precision 
of 2-3 xlO -6 cm -1 were found. This indicates that the differences between 
the values determined from ATMOS and the standard positions provided by 
Brown are statistically significant. 

The absolute error in estimating Doppler shifts can now be estimated 
since instrumental lines were available in the spectra for calibration. The 
rest frame positions ? of these lines were determined from the spectra to 
precisions of about 3 x 10 -B cm -1 . For 20 of these lines and using an average 
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frequency of about 2350 cm -1 the absolute uncertainty in relative velocity v 
is determined by 

3 x 10 -6 c 

absolute uncertainty = — j= < 1 ms -1 . (2.10) 

2350 y 20 

The instrumental lines allow us to calibrate the instrument to obtain absolute 
Doppler shifts to < 1 ms -1 provided the number of spectra n and the number 
of lines N are large enough so that the factor y/nN is large enough to reduce 


the mean uncertainties. 
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Figure 2.1.1 Instrumental 1/3 CO a absorption lines and 0 CO v = 0 — ♦ 1 band 
head observed in solar spectra. The CO a instrumental absorption lines 
are used for calibration. The 0 CO lineB tend to interfere with atmospheric 
and instrumental lines causing unwanted line distortions. 



Figure 2.1.2 Instrumental 1/3 H a 0 absorption lines intermixed 0 CO lines 
observed in solar spectra. The stronger lines as well as some weaker lines 
correspond to the instrumental lines. The water molecule is not linear and 
its spectra does not display evenly spaced rotational levels like CO a . Water 
lines must be carefully checked for possible interferences with 0 CO lines. 
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Figure 2.1.3 The test for calibration drift S k using v 3 C0 3 instrumental lines 
for the five occultations. The data points appear to be randomly scattered 
about the origin with a variance of about 5 x 10 -s and this is due to mea- 
surement error associated with locating line positions. 
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Figure 2.1.4 The test for alibration drift S k using i/ 3 H 3 0 instrumental lines for 
the five occultations. The data points appear to be randomly scattered about 
the origin with a variance of about 8 x 10 -8 and this is due to measurement 
error associated with locating line positions. 
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Figure 2.1.5 Measurement errors <r(6 k ) of the S k values in Fig. 2.1.3 using 
CO] instrumental lines for the five occultations. The scatter per measure- 
ment is of the same order of magnitude as the scatter of values in 2.1.3. The 
scatter of points in 2.1.3 are due to random errosr in the line positions. 



Figure 2.1.6 Measurement errors <r(6 k ) of the S k values in Fig. 2.1.4 using i/j 
H a O instrumental lines for the five occultations. The scatter per measure- 
ment is of the same order of magnitude as the scatter of values in 2.1.4. The 
scatter of points in 2.1.3 are due to random errors in the line positions. 
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Table 2.1.1 Check for instrumental drift for the five occultations obtained 
from i/ 3 C0 2 lines. 


Occult ation 

£(xl0 -8 ) 

<r(£)(xl0" 8 ) 

Ref. Spectrum 

SR02 

-0.72 

1.96 

A50430113052 

SS06 

-0.62 

2.15 

A50430230523 

SS09 

-0.88 

2.30 

A50501172502 

SS11 

+0.58 

3.01 

A50501202822 

SS13 

-0.29 

3.11 

A50501233132 


Table 2.1.2 Check for instrumental drift for the five occultations obtained 
from i/j HjO lines. 


Occultation 

£(xl0~ 8 ) 

<r(£)(xl0~ 8 ) 

Ref. Spectrum 

SR02 

-4.70 

3.55 

A50430113052 

SS06 

+0.80 

3.95 

A50430230523 

SS09 

-1.02 

4.12 

A50501172502 

SS11 

+1.20 

3.57 

A50501202822 

SS13 

+1.14 

3.94 

A50501233132 
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Table 2.1.3 v$ COj standards given by Brown (1986) used in primary instru- 
mental calibration. 


No. 

Position (cm *) 

Rot 

No. 

Position (cm 1 ) 

Rot 

1 

2330.5571 

P22 

16 

2364.1053 

R20 

2 

2334.1563 

P18 

17 

2365.3876 

R22 

3 

2335.9140 

P16 

18 

2366.6451 

R24 

4 

2337.6580 

P14 

19 

2369.0853 

R28 

5 

2339.3722 

P12 

20 

2370.2679 

R30 

6 

2341.0620 

P10 

21 

2371.4255 

R32 

7 

2342.7272 

P 8 

22 

2372.5581 

R34 

8 

2345.9841 

P 4 

23 

2373.6656 

R36 

9 

2354.4333 

R 6 

24 

2374.7481 

R38 

10 

2355.8894 

R 8 

25 

2375.8055 

R40 

11 

2357.3207 

R10 

26 

2376.8377 

R42 

12 

2358.7272 

R12 

27 

2377.8448 

R44 

13 

2360.1089 

R14 

28 

2378.8268 

R46 

14 

2361.4659 

R16 

29 

2379.7836 

R48 

15 

2362.7980 

R18 





Table 2.1.4 i/j H 2 O standards given by Brown (1986) used in primary 
instrumental calibration. 


No. 

Position (cm x ) 

Vib-Rot 

No. 

Position (cm -1 ) 

Vib-Rot 

1 

1616.7116 

110 101 

12 

1756.8188 

716 625 

2 

1627.8275 

202 111 

13 

1771.2875 

331 220 

3 

1647.4041 

422 413 

14 

1775.6341 

726 615 

4 

1669.1683 

413 322 

15 

1790.9518 

827 716 

5 

1669.3929 

313 202 

16 

1799.6156 

431 322 

6 

1688.3785 

725 716 

17 

1802.4797 

918 827 

7 

1704.4534 

432 423 

18 

1825.2016 

634 523 

8 

1714.0337 

625 616 

19 

1844.1806 

441 330 

9 

1715.1551 

606 515 

20 

1867.8527 

542 431 

10 

1718.6117 

322 211 

21 

1889.5694 

643 532 

11 

1730.0550 

615 524 

22 

1942.5161 

652 541 
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Table 2.1.5 Estimated global offsets a for the five occultations from 1/3 C0 3 
lines. Lines 1-20 were used for this analysis. 


Occultation 

q(x10 -8 ) 

O’(a)(xl0 8 ) 

Av (m/s) 

a( Av) (m/s) 

SR02 

+2.66 

0.39 

+7.97 

1.17 

SS06 

-0.55 

0.48 

-1.65 

1.42 

SS09 

-27.30 

0.43 

-81.84 

1.29 

SS11 

-7.13 

0.63 

-21.38 

1.89 

SS13 

-2.96 

0.52 

-8.87 

1.56 


Table 2.1.6 Estimated global offsets a for the five occultations from i/ 3 H 3 0 
lines. 


Occultation 

a(xl0 -8 ) 

cr(a)( xlO -8 ) 

Av (m/s) 

<r(Av) (m/s) 

SR02 

-0.53 

0.74 

-1.59 

2.22 

SS06 

-2.88 

0.88 

-8.63 

2.64 

SS09 

-26.32 

0.78 

-78.91 

2.34 

SS11 

-8.56 

0.77 

-25.66 

2.31 

SS13 

-3.87 

0.67 

-11.60 

2.01 


Table 2.1.7 Differences between global offsets a for the five occultations. 


Occultation 

Aa(xl0 -8 ) 

*(Aa)(xlO- 8 ) 

SR02 

+3.19 

0.84 

SS06 

+2.33 

1.00 

SS09 

-0.98 

0.89 

SS11 

+1.43 

0.99 

SS13 

+0.91 

0.85 
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Table 2.1.8 v 3 C0 2 standards given by Brown (1986) as compared to ob- 
served instrumental positions after primary calibration for n = 125 spectra. 
All frequencies in cm -1 where applicable. 


No. 

SBrown 

Qnit 

Ac (xlO -4 ) 

*(&ut) (xlO -4 ) 

1 

2330.5571 

2330.55709 

-0.1 

0.20 

2 

2334.1563 

2334.15633 

+0.3 

0.19 

3 

2335.9140 

2335.91938 

-0.2 

0.18 

4 

2337.6580 

2337.65799 

-0.1 

0.20 

5 

2339.3722 

2339.37219 

-0.3 

0.21 

6 

2341.0620 

2341.06197 

-0.3 

0.21 

7 

2342.7272 

2342.72719 

-0.1 

0.24 

8 

2345.9841 

2345.98410 

0.0 

0.34 

9 

2354.4333 

2354.43342 

+1.2 

0.26 

10 

2355.8894 

2355.88941 

+1.0 

0.22 

11 

2357.3207 

2357.32067 

-0.3 

0.22 

12 

2358.7272 

2358.72723 

+0.3 

0.19 

13 

2360.1089 

2360.10892 

+0.2 

0.18 

14 

2361.4659 

2361.46586 

-0.4 

0.17 

15 

2362.7980 

2362.79803 

+0.3 

0.19 

16 

2364.1053 

2364.10525 

-0.5 

0.20 

17 

2365.3876 

2365.38761 

+0.1 

0.20 

18 

2366.6451 

2366.64511 

+0.1 

0.20 

19 

2369.0853 

2369.08528 

-0.2 

0.25 

20 

2370.2679 

2370.26791 

+0.1 

0.29 

21 

2371.4255 

2371.42548 

-0.2 

0.31 
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Table 2.1.9 t/ 2 HjO standards given by Brown (1986) as compared to ob- 
served instrumental positions after primary calibration for n = 125 spectra. 
All frequencies in cm -1 where applicable. 


No. 

(Brown 

(init 

A? (xlO -4 ) 

(xlO -4 ) 

1 

1616.7116 

1616.71148 

-1.2 

0.20 

2 

1627.8275 

1627.82747 

-0.3 

0.26 

3 

1647.4041 

1647.40405 

-0.5 

0.24 

4 

1669.1683 

1669.16834 

+0.4 

0.34 

5 

1669.3929 

1669.39287 

-0.3 

0.21 

6 

1688.3785 

1688.37845 

-0.5 

0.38 

7 

1704.4534 

1704.45341 

+0.1 

0.21 

8 

1714.0337 

1714.03360 

-1.0 

0.30 

9 

1715.1551 

1715.15508 

-0.2 

0.21 

10 

1718.6117 

1718.61159 

-1.1 

0.22 

11 

1730.0550 

1730.05499 

-0.1 

0.34 

12 

1756.8188 

1756.81878 

-0.2 

0.24 

13 

1771.2875 

1771.28742 

-0.8 

0.25 

14 

1775.6341 

1775.63409 

-0.1 

0.42 

15 

1790.9518 

1790.95184 

+0.4 

0.31 

16 

1799.6156 

1799.61557 

-0.3 

0.30 

17 

1802.4797 

1802.47976 

+0.6 

0.44 

18 

1825.2016 

1825.20154 

-0.6 

0.23 

19 

1844.1806 

1844.18061 

+0.1 

0.18 

20 

1867.8527 

1867.85279 

+0.9 

0.33 

21 

1889.5694 

1889.56936 

-0.4 

0.27 

22 

1942.5161 

1942.51614 

+0.4 

0.36 





Figure 2.1.8 a k determined from COj (solid line) and H 3 O (dotted line) 
instrumental lines for occultation SS06. The a* values for COj are generally 
greater than the H 2 0 values. 
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Figure 2.1.9 a k determined from COj (solid line) and HjO (dotted line) 
instrumental lines for occultation SS09. The a k values for COj are generally 
greater than the H 2 O values. 




Figure 2.1.10 a k determined from COj (solid line) and HjO (dotted line) 
instrumental lines for occultation SSll. The a* values for CO a axe generally 
greater than the HjO values. , 
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Figure 2.1.11 a* determined from C0 3 (solid line) and H 3 0 (dotted line) 
instrumental lines for occultation SR13. The a* values for CO] are generally 
greater than the H 3 0 values. 



Figure 2.1.12 <r(a k ) from C0 2 instrumental lines for the five occultations. 
The mean uncertainties are typically about 2 x 10"* which is typically the 
magnitude of the fluctuations of the a* in Figs. 2.1.9-13 for C0 3 . 
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Figure 2.1.14 A(a fc ) between CO a and H a 0 instrumental lines for the five 
occultations. The luge fluctuations are due to measurement errors, how- 
ever, the average of the differences tends to be larger than zero. The mean 
difference is +1.38 ± 0.41 x 10 -8 for all occultations. 
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II. 2 Scale Uncertainty a(e) 


In the previous section we have calibrated the spectra using about 30 1/3 
instrumental lines of CO3. Global calibration errors are not the only source 
of errors in the measurement of line positions. The scale uncertainty cr(e) sets 
the ultimate limit to the precision by which wind speeds may be obtained. 
To estimate this scale uncertainty, consider (2.1) after correcting all spectra 
by the calibration factor 1 — a §11.1. We have 

c} = fj + e) + e*. (2.11) 

which can be rearranged 

<*=<*-<,- «?■ ( 2 - 12 ) 

Averaging this over all lines j in each spectrum k gives 

«*=<<?- ti) s - <«;>,. (2-i3) 

For a sufficiently large sample size we expect that measurement errors due 
to noise vanish 

0. (2.14) 

If we impose this condition then (2.12) becomes 

«* = (2-i5) 

The COj and HjO line positions determined from the ATMOS spectra 
in Tables 2. 1.8-9 were used for the standard positions q, for both regions. In 
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Figures 2.2. 1-2 values for 6* from 125 spectra are plotted for both spectral 
regions respectively. The estimated values of the scale errors e* appear to 
scatter between 5 and 10x10"® cm" 1 . Also plotted is the mean uncertain- 
ties in measuring the scale error <r(e*). We see that the error in the e* is 
approximately the same as their scatter between spectra. Values of the scale 
error can be determined by defining the average of the scale errors for each 
occupation 

(2.16) 

for each occupation. Tables 2.2. 1-2 list e and a(e) for the C0 3 and H 3 0 
lines respectively for each of the five occupations. It is evident that |e| < 
<x(e) in both regions. This indicates that the scale uncertainty is less than 
5 x 10" B cm _1 . It should also be noted that ir(e) ~ <r(c*). This indicates 
that the measured values e* may in fact be due to noise instead of being 
due to a true scale uncertainty. More precisely stated, the approximation 
(e*) may not vanish in each spectrum but may have a random fluctuation 
due to the finite sample size of about 20 lines. To test this we examined the 
correlations between the estimates for e* in each spectrum k for the H 3 0 and 
C0 3 regions. The correlations 

C k ~ ( 2 - 17 ) 

were determined in each spectrum and displayed in Figure 2.2.3 for all the 
spectra k. The correlations appear to be evenly distributed about zero. This 
suggests that there a correlation does not exist. The total correlation 


C = (C\ 


(2.18) 
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was then determined and values for C and <r(C) are listed in Table 2.2.3. 
The correlation overall correlation C is of 0(3 x lO -10 ) cm -3 and and the 
mean uncertainties <r(C) is also 0(3 x lO -10 ) cm -1 . The upper limit of 
the correlations is therefore determined by the uncertainties o(C). Since the 
correlation is determined from the positions of two lines, the upper limit of the 
scale uncertainties is approximately determined by o-(C') ~ 2<r 3 (e). Using the 
values of <r(C7 ) in Table 2.2.3 this gives an upper limit of the scale fluctuations 
of about 1.3 x 10 -5 cm -1 . This means that the frequencies obtained by the 
interferometer are good to at least 10 -6 cm -1 . The results of the correlation 
analysis do not show any significant evidence of a calibration uncertainty. 
Since the correlations effectively vanish the values of e* plotted in Figures 
2.2.1-2 are due to measurement errors e* and not a true scale uncertainty. 
As previously stated, due to the number of lines sampled, N = 20, there may 
still be significant fluctuations in the average of the measurement errors of 
all the line positions obtained in spectrum k 

«* = <e‘). (2-W) 

of the order 

a(e k ) = N- X ' 2 <r{e). (2.20) 

Using typical values for <x(e) listed in Table 1.12.3, for example 2.65 x 
10 -4 cm -1 for peak absorptances of about 20% and using N = 20 we get 
an estimation for <r(e fc ) = 5.9 x 10 -4 cm -1 . This agrees with the values esti- 
mated for cr(e*) in Tables 2.2. 1-2. 

This suggests there is no evidence for the presence of a scale uncertainty. 
The estimate of the scale uncertainty in §1.1 predicted by the magnitude 
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of the phase errors is several orders of magnitude smaller that the estimates 
£r(e*). From here on it is assumed that e* = 0 and that all fluctuations are due 
to measurement uncertainty. This result simplifies the analysis significantly. 
The scale uncertainty appears to be negligible and the relative accuracy of 
the wind speed measurements is not limited by this effect. 
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Figure 2.2.1 e* for the five occultations determined from COj lines (points) 
and <r(e*) (lines). The scatter of the data about zero is of the same order 
of magnitude of the mean uncertainties. This indicates that the scatter in 
points is probably due to noise instead of a true scale fluctuation. 
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Figure 2.2.2 e* for the five occultations determined from HjO lines (points) 
and ?(e*) (lines). The scatter of the data about zero is of the same order 
of magnitude of the mean uncertainties. This indicates that the scatter in 
points is probably due to noise instead of a true scale fluctuation. 
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Figure 2.2.3 Correlation C k of the CO 2 and HjO e* determined from all 
five occultations. The scatter is within 0.5 xlO -8 cm -3 . The scatter of this 
correlation about zero means that the scale uncertainty <re)is less than 2 x 
10 -5 cm -1 . 
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Table 2.2.1 e determined from C0 2 lines. 


"'cm' 1 ) 

<r(e)(xl0 -5 cm x ) 

.083 

4.533 

.114 

4.920 

.090 

5.269 

.095 

6.767 

.094 

7.184 


Table 2.2.2 e from */ 2 H 2 0 lines. 


Occultation 

e(xl0 'em" 1 ) 

dr(e)(xl0~ 5 cm x ) 

SR02 

-3.356 

6.521 

SS06 

-1.766 

6.820 

SS09 

+3.964 

7.096 

SS11 

-0.258 

6.125 

SS13 

+0.599 

6.721 


Table 2.2.3 Correlations of e* between H 2 0 and C0 2 regions 


Occultation 

C(xl0 _lo cm" a ) 

<r((7)(xl0 10 cm" J ) 

SR02 

+0.838 

6.507 

SS06 

-2.168 

6.074 

SS09 

-12.62 

7.267 

SS11 

+7.728 

5.809 

SS13 

-6.488 

8.475 

all 

-3.362 

3.346 
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II.3 Solar and Ratioed Spectra 


Figure 2.3.1 shows a condensed full Filter-3 solar spectrum collected at 
tangent heights above 120 km where there is not any significant atmospheric 
absorption. The spectral features are condensed due to the large range of 
wavenumbers 1600 < ? < 3200 cm" 1 shown in a small space on the page. 
Many spectral features are evident, particularly the solar ‘O’ CO fundamental 
(Av = 1), - e.g. ti = 0-M) = l, ti = l-ni = 2, etc. - rotation- 
vibration bands with wavenumbers less than 2330 cm -1 . Figure 2.3.2 shows 
a ‘magnification’ of the © CO band heads near 2300 cm" 1 along with the 
i/ 3 COj instrumental lines. Instrumental water lines protrude below 2000 
cm -1 . The lines in the spectrum above the v 3 CO a region c > 2380 cm" 1 
are other solar lines. The 0 CO bands can distort other spectral features 
which may be of interest. The instrumental C0 2 lines cannot be used for 
calibration in the region where the v 3 COj lines overlap with the © CO lines 
since interferences distort the line shapes. Distortions can cause systematic 
errors in determining their positions ?. The lines used in this report were 
carefully selected to avoid these problems. 

The spectral regions where the absorption bands used to obtain wind 
speeds in the atmosphere are located in Figure 2.3.1. All of the lines in 
the solar spectra appear to be superimposed on a ‘background signal’ S. 
The background signal appears to vary with wavenumber <;. This happens 
because the detector sensitivity and the filter transmittance vary with the 
wavenumber [Shaw, 1989]. The background signal as well as the instrumental 
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and solar lines superimposed on it serve as the input radiation intensity /^(O) 
described in §1.9 (1.144). The solar features are present in all of the ATMOS 
spectra, with or without atmospheric absorption. 

The positions and shapes of the solar features and instrumental lines re- 
mained essentially unchanged during each occultation. It was desirable to 
‘factor out’ these lines since they distort atmospheric features. The solar fea- 
tures serve as the background J t (0) in (1.144) for the atmospheric absorption 
exp[— k c u] as described in §1.9. The signal J t (0) was estimated by averaging 
all of the solar spectra taken above the atmosphere for each occultation [Nor- 
ton, 1986]. This produced solar spectra / c (0) with a very high SNR ~ 500:1 
[Norton, 1986]. The atmospheric absorption was then estimated by taking 
the ratio of the spectra containing atmospheric features J t to the solar av- 
erage I 5 ( 0). These were termed ‘ratioed’ spectra and from (1.144) they are 
defined to be 



The ratioed spectra made by Norton (1986) were used for the majority of 
the atmospheric studies conducted by the ATMOS team. 
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3200 

Figure 2.3.1 A complete Filter-3 solar spectrum. Instrumental lines of COa 
and HaO are superimposed on background solar spectra. Solar CO lines are 
present below ~ 2300 cm -1 and the variable output of the detector is evident. 



:220 .00 WflVEfO 2380.0 

Figure 2.3.2 Av = 1 0 CO band heads and instrumental 1/3 COa lines ob- 
served in the solar spectra. This covers a slighly large wavenumber range 
than Figure 2.1.1. 
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II.4 Estimating the Noise Level 


The interferograms included noise in the instrumental signals caused by 
fluctuations of the light source, photon noise caused by random arrival of 
photons from the source, thermal fluctuations in the detector, and external 
vibrations of the instrument [Norton, 1986; Thorne, 1988], The noise in the 
interferograms I(x) is transferred into noise in the spectra /'(c) through the 
transform (1.9). In §1.11 it was shown how this noise affects the precision by 
which the positions of the absorption lines can be measured. 

The ratioing process does not affect the signal to noise ratio SNR* of 
the spectra significantly. In the non-ratioed spectra the noise A/"* in the 
signal J*(0) has random fluctuations er(A/’*). The ratioing process produces 
a constant background signal of unity by definition (2.20). The noise TV* in 
the ratioed spectra iZ* then has fluctuations 

•w = T§- < 2 - 22 > 

In regions where solar absorption features are not present the background 
signal S* and solar signal /*( 0) are the same. The SNR of the ratioed Bpec- 
tra is dependent on the background signal in the denominator of (2.21). In 
regions where solar absorption features are strong the denominator of (2.21) 
may become large and the noise increases. To get the best results in measur- 
ing line positions in the ratioed spectra iZ* it is desirable to avoid spectral 
regions where the noise level caused by solar features is high. An estimate of 
the SNR* of the ratioed spectra iZ* can be determined by examining ratioed 
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spectra above tangent heights where atmospheric absorption is significant. 
Figure 2.4.1 shows one such spectra. We see that between 1600 and 2000 
cm -1 , for example, the SNR is lower than near 2400 cm -1 . By visually com- 
paring the wavenumber variations in fluctuations of the ratioed spectrum in 
Figure 2.4.1 to the wavenumber variations of the background signal in the 
unratioed solar spectrum in Figure 2.3.1 the correlation between the two is 
evident. This is predicted by (2.21). 

Since scale errors e* are negligible, measurement uncertainties of the line 
j in the spectrum k can be determined from (2.10) 



(2.23) 


By using these values, from the theory §1.11 


A, M‘j) * < r (( A ii c i) k ) = /»®W •. 


(2.24) 


where <r(JV) = 1/SNR, is the noise to signal ratio and is described by (2.21) 
in the ratioed spectra jR,. By measuring the positions of instrumental lines 
of HjO and C0 2 in the unratioed spectra we can estimate the signal to noise 
ratios for both regions in the ratioed spectra by using (2.23). We expect the 
relation (2.21) to hold true. From examination of Figure 2.3.1 the SNR of 
the H a O region should be almost a factor of two smaller than the C0 2 region. 

The line noise will be defined here to be f n <r(N) w, f n being the form 
factor, <r(N) = 1/SNR the signal fluctuations and w the geometric width 
of the line as described in §1.11. For a given spectral region the signal fluc- 
tuations <r(N) are approximately a constant and since / n and w are the 
same for similar lines, the line noise defined in (2.23) is approximately a con- 
stant. The line noise over all lines j in each region were determined using 
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(2.23) and are listed in Tables 2.4.1-2 for each occultation individually and 
all collectively for the CO3 and H3O instrumental lines respectively. The 
line noise for the C 0 3 region is almost a factor of two larger than that the 
H3O regions in agreement with (2.21). For the C 0 2 instrumental lines the 
line noise of about 0.47 x 10 -B cm -1 is in close agreement with the estimate 
(2.23) predicted with rough estimates of the line parameters f n = 0.35355, 
ar(N) = 0.01, and W = 0.015 cm -1 of 0.54 x 10 -4 cm -1 . It appears that the 
measurement uncertainty theory can be used to describe the precision of the 
line-finder program quite accurately. 
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Figure 2.4.1 The noise level in the spectra after the ratioing process. This 
spectrum was taken at a tangent height where atmospheric absorption is 
absent. This represents the background noise level of the spectra. Note the 
variations in the spread of the noise with wavenumber. 



Table 2.4.1 Estimates of noise level from CO 2 instrumental lines. 


Occultation 

mean (xl0 -8 cm -1 ) 

SR02 

4.94 

SS06 

4.55 

SS09 

4.42 

SS11 

4.67 

SS13 

4.79 

all 

4.67 


Table 2.4.2 Estimates of noise level from HjO instrumental lines. 


Occultation 

mean (xlO 5 cm -1 ) 

SR02 

8.30 

SS06 

8.01 

SS09 

7.79 

SS11 

8.10 

SS13 

8.09 

all 

8.10 
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II. 5 Estimating Relative Speeds Along the Line of 

Sight 


The method of estimating individual line positions and the calibration 
stability of better than 2 x 10 -B cm -1 determined in §11.3 imply that wind 
speeds can be obtained to a precision of 5 ms -1 or better. It was also pointed 
out by Shaw (1985) that these kinds of precisions could be achieved using 
Filter 3 spectra provided that the SNR was large enough. From the data in 
Table 1.12.3 we can estimate the sample size needed to attain the precision 
of 5 ms -1 or better. From (1.186) the uncertainty per line j of the velocity 
is given by 

<r(vj) = <r(f,)|-. (2.25) 

v 

From the measurement theory (1.180) the uncertainty is proportional to A~*. 
For a line with a peak absorptance A p ~ 20% and a position of ~ 2350cm -1 
we see that <r(vj) ~ 35 ms -1 . Using N lines the precision of the relative 
velocity is then 35 ms -1 /\/V so we need about 50 lines to attain 5 ms -1 . 
However, for stronger lines, using a peak absorptance Ap ~ 70% for example, 
this sample size can be considerably smaller. By using a value of o’(c'j) — 
0.76 x 10 -4 cm -1 we find that <r(vj) ~ 10 ms -1 which requires a sample 
of size of only 4 lines to attain the desired precision. In practice the peak 
absorptances of the lines used are primarily in the 20-70% range and a sample 
size between 20 and 30 lines is usually available. Therefore we can predict 
that relative speeds can be determined to at least 5 ms -1 and most likely 
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better using single bands of atmospheric lines. 

The accuracy by which line Doppler shifts can be determined is inde- 
pendent of any wavenumber dependent calibration errors, provided that line 
positions in the rest frame can be measured by the instrument. To see this, 
the wavenumbers measured by the instrument can in principal be shifted 
from the “true” wavenumber Ctrue by a small amount and this offset can be 
represented by 

Out = (1 4" ®) Ctrue (2.26) 

where a is may be regarded as constant over small spectral regions. Note 
that this a is not the same factor defined in §11.1. If the “true” rest frame 
position is given by Strue then the corresponding position as measured by the 
instrument will be at shut = (1 + a) <ftn»e- For Doppler shifted lines then 

<nntt = (1 4" o)?true = (1 4" o) ^1 4“ Strut 

= (l + ^) Shut (2.27) 

which is independent of the shift a. This is because the form 14- a is the same 
for both the rest frame and moving frame positions. This result indicates 
that accurate relative speed measurements can be made provided that the 
rest frame positions Shut are known and that the calibration of the ATMOS 
interferometer remains stable throughout the course of observation. 

There are typically about 30 lines of the 1/3 COj band which are visible in 
the spectra and can used to measure Doppler shifts. Relative speeds can be 
obtained to 5 ms -1 with these lines provided that the lines are not too weak 
or saturated. From (2.24) we see that the uncertainty of measuring Doppler 
shifts is wavenumber dependent. As an example, for the CO 3 band near 
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667 cm -1 , many more lines are needed. For a wavenumber of about 1/4 of 
the v 3 lines the uncertainty er(vj) will be about 4 times as large. We would 
need 16 times as many lines. From a practical point of view, relative speed 
measurements should be confined to the higher wavenumber regions of the 
spectra. 


II. 6 Relative Speed Measurements Using i/ 3 CO 2 
Lines Between 2330 and 2372 cm -1 


The ATMOS instrument was found to have a stable frequency calibra- 
tion as well as an excellent SNR of about 100:1. This enables precise relative 
speeds to be obtained from a small number of lines in the mid-infrared re- 
gions [Shaw, 1985]. Below about 120 km atmospheric absorption lines of the 
i/ 3 C0 2 absorption band become visible in the spectra between 2330 and 2370 
cm -1 . Figure 2.6.1 shows several of these lines along with their instrumental 
counterparts taken from an unratioed spectrum for illustration. These at- 
mospheric lines were used to measure relative speeds from tangent heights of 
about 120 km to about 80 km where overlapping with other C0 2 bands and 
line saturations occur. A typical overlap of C0 2 lines from different bands 
taken in successive spectra is shown in Figure 2.6.2. As the tangent height 
decreased the lines begin to interfere with each other. The spectra are dis- 
placed vertically for clarity. The v 3 lines used in this study were chosen to 
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be those lines which displayed minimal interferences. These were considered 
to be “good” lines. Another criterion used was that the peak absorptances 
range from 10-70% because resolved line weaker than these became distorted 
by the background noise while stronger lines become saturated. This caused 
the curvature of the line near the line center to become flat which increased 
the measurement error (the geometrical or spectral width w of the line in- 
creases dramatically). 

With the above criteria for choosing absorption lines, the instrumental 
lines given in Table 2.1.8 were used as the standards in (2.26). The relative 
speed for v k for each spectrum k and line j was measured using (1.1) 



From this the relative speed v k for each spectrum k was determined 

•* = <»‘)j (2-29) 

and its corresponding mean deviation <r(v k ) determined. This was performed 
for each of the occultations and the resulting relative speed profiles were 
determined. Figures 2.6.3-4 show the relative speed versus ray tangent height 
using the i/ 3 CO3 atmospheric lines for the sunrise and sunset occultations 
respectively. The relative speed profiles are typically of the order of 6000 
ms -1 . This is due to the orbital motion of the satellite and and this motion 
is estimated and subtraced from the profiles to obtain zonal wind in the 
atmosphere around the tangent point §11.12. 

The relative speed profiles show a variability with tangent height. The 
‘wavelike’ structure of the relative speed profiles is immediately recognized, 
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particularly in Figure 2.6.4 where the four sunset profiles demonstrate a 
remarkable similarity with regards to altitudes of peaks and valleys. Large 
shears, especially between 105 and 120 km are apparent. The sunrise profile 
in Figure 2.6.3 demonstrates similar features. Figure 2.6.5 displays the mean 
error <r(v*) for all of these occupations. Note that <r(t/ fc ) increases at large 
heights because the lines become very weak. At the lowest tangent heights 
the increase in ir(t; fc ) is due to overlapping of lines and saturation. The curves 
all follow the same general trend and, as proposed above, precisions of 5 ms -1 
are indeed found. Between 80 and 110 km these lines achieved a maximum 
precision of ~2.5 ms -1 . Since the calibration of the instrument is stable 
these precisions indicate that the features observed are real. Also evident in 
Figure 2.6.4 are relative differences of over 100 ms -1 between SS06 and the 
rest of the sunset profiles. The relative speed profiles for SS09, SSll, and 
SS13 showed a trend of increasing relative speed with respect to the SS06 
profile as the mission progressed. This is due to the changing geometrical 
orbit of the instrument and its component of velocity along the line of sight. 
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Figure 2.6.1 Pairs of instrumental and atmospheric lines. The atmospheric 
lines are red-shifted (towards the left) from the instrumental lines as the 
tangent height decreases from 112 (a) to 90 km (c) in a sunset occultation. 
The spectra are displaced for clarity. Note the increasing peak absorptions 
of the atmospheric lines as the tangent heights decreases. 



Figure 2.6.2 Overlapping of lines from ATMOS spectra. Notice how the line 
shapes become distorted by other spectral features as the tangent height 
decreases. Lines must be examined for interferences before including them 
in wind measurements. The spectra have been slightly displaced vertically 
for clarity. 
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Figure 2.6.3 Relative speed profile determined using v$ C0 3 atmospheric 
lines for the sunrise occultation. There appear to be rather large shifts in 
the relative speeds near 80 and 105 km. 



Figure 2.6.4 Relative speed profile determined using v$ CO a atmospheric 
lines for the sunset occultations. The four occultations (6) SS06 • (9) SS09 
- (11) SSll - (13) SS13 are displaced due to changes in the orbital motion 
§11.12. The appear to have similar relative speed shears near 80 and 105 km 
as well as ‘minima’ near 95 km. These are due to variations in the upper 
atmospheric wind speeds. 
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Figure 2.6.5 Measurement uncertainty &{v k ) for all occultations using U 3 CO3 
atmospheric lines. A precision of better than 5 ms -1 is achieved between 80 
and 110 km. The increase outside of this region of high precision is due to the 
weakness of lines above 110 km and line saturation and interferences below 
80 km. 
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II.7 Extending Relative Speed Measurements and 
the Determination of v- msi 


We now proceed to extend these measurements to lower altitudes us- 
ing absorption lines in which no instrumental counterparts are available for 
determining the rest frame positions. The profiles in Figures 2.6.3-4 were 
determined by directly comparing the instrumental COj lines to their at- 
mospheric counterparts and thus directly measuring the Doppler shifts. In 
general, instrumental lines are not available for direct comparison at lower 
tangent heights. Fortunately there is a way of getting around this problem. 
We know that the calibration of the instrument is stable, we have determined 
relative speeds using the calibrated instrumental lines, and there are more 
lines entering into the spectra as the tangent height decreases. We must use 
these new atmospheric lines that are entering into the spectra, however we 
need to know for each of these lines. 

Consider a set of new atmospheric lines which have entered the spectra. 
Their positions are Doppler shifted by (1.1). It should now be understood 
that we are referring our wavenumbers to the calibration system of the in- 
strument and are no longer referring to ‘absolute standards’ which were used 
to calibrate the v 3 COj region. As shown in the previous section (2.26) 
absolute standards are no longer needed to measure absolute relative wind 
speeds. Let us suppose that the relative speeds v k for a set of spectra k have 
been determined to a reasonable degree of precision (5 m/s or better). Then 
we can invert (2.26) to obtain an estimate of the rest frame positions of each 
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of the new lines which have appeared in the spectra for which v k had been 


estimated 


3 1 + v k /c \ c ) 3 


( 2 . 30 ) 


from the observed atmospheric line position for each new line j in each 
spectrum k. We can estimate the rest frame position from the average 


* = <$)* 


( 2 . 31 ) 


as well as its mean error <r(fy). Having determined these rest frame positions 
the new lines may be used to extend the wind below the current minimum 
tangent height. 


One problem with this technique as it is presently stated is that the line 
positions % determined axe only as precise as the relative speed profiles v k 
from which they were obtained. On the other hand, if the rest frame positions 
have a degrading precision then the relative speeds themselves will become 
less precise. An iterative scheme proved to be useful in dealing with this 
problem. Specifically, an initial set of relative speeds was determined 
from an initial set of rest frame positions over am altitude range for all 
of the occultations. From these initial relative speeds an updated list 


of line positions was determined 

i *.(0 _ S 

J 1 -f 



Then an updated relative speed profile was determined 



( 2 . 32 ) 

( 2 . 33 ) 

( 2 . 34 ) 

( 2 . 35 ) 
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This procedure was then repeated until a desired level of convergence in the 
mean uncertainties of the relative speeds &(v k ) and the line positions d’(cj) 
were both obtained. I chose value of 0.1 ms -1 and 10 -8 cm -1 respectively. 

The relative speed profiles shown in Figures 2. 6. 3-4 were recomputed 
along with the rest frame positions of the atmospheric lines by using the 
above procedure for the lines specified in Table 2.1.8. The precisions in lines 
positions o’(Cj) were found to converge to less than 10~ 6 cm -1 in usually 
three or four iterations with a simultaneous convergence of <r(v k ) less than 
0.1 ms -1 . Figure 2.7.1 shows the precisions of the positions determined from 
measuring instrumental lines (bold line) and those positions computed from 
the atmospheric lines using the iteration procedure (points). The precisions 
determined using the atmospheric lines are more uniformly distributed with 
a value near 0.28 x 10 -4 cm -1 than the positions determined from the instru- 
mental lines. The apparent degradation of the precision is actually due to 
the smaller sample size of atmospheric spectra used (about 60) compared to 
the number of high altitude spectra used (about 125). However, the results 
do show that by using the velocity profiles to determine line positions, we can 
obtain satisfactory results, particularly in the case where only a few atmo- 
spheric lines are originally present for initial line position estimates. Figure 
2.7.2 shows the precisions obtained using instrumental lines (bold lines) and 
atmospheric lines (points). The slight improvement below 100 km is not 
significant. 

At wavenumbers between 2372 and 2380 cm -1 eight v 3 COj lines become 
reasonably strong (A p > 15%) at about 100 km and are useful for extending 
the relative speed profiles. Their positions were estimated from (2.27- 
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8) using the relative speed profiles in Figures 2. 6.3-4 above 80 km. The 
line positions were then used to extend the relative speed profiles to lower 
tangent heights. The iteration procedure was used for the lines between 
about 70 and 100 km. Figure 2.7.3 shows the improvements in the lines 
precisions between the first estimate (solid line) and the converged result 
(points). An improvement of about 5 x 10 -6 cm -1 , although small, is evident. 
The precisions in the relative speeds are shown in Figure 2.7.4 and there does 
not appear to be a significant difference before and after the iteration. Since 
the range of altitudes where these COj lines appear overlap the relative 
speeds in Figures 2.6.3-4 is large (30 km) the sample size of lines for the new 
rest positions is large enough to ensure fairly precise initial estimates cj°^ of 
about 3xl0 -6 cm -1 . This is not the case in general, however. 

The above results show that the initial relative speed profiles can be used 
to determine the unknown rest frame positions of atmospheric lines. In turn, 
these rest frame positions can then be used to determine relative speeds at 
extended tangent height ranges. Since the line positions and relative speeds 
are determined simultaneously, the precisions of the initial relative speed 
profiles are maintained. This indicates that by using the calibration stability 
of the instrument, relative speeds with precisions and accuracies similar to 
those determined by using the positions measured from the instrumental lines 
can be determined from other atmospheric lines. We used this procedure to 
extend relative speed profiles to 25 km by using NjO, CO a and H a O lines 
§11.9. 
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Figure 2.7.1 Precisions of t / 3 C0 2 rest frame positions determined from at- 
mospheric lines (v < 2372 cm -1 ) before (line) and after (points) iteration 
procedure. The improvement in line positions is slight. 



Figure 2.7.2 Precisions of relative speeds determined from 1/3 COj atmo- 
spheric lines (v < 2372 cm -1 ) before (line) and after (points) iteration pro- 
cedure. There is only a slight improvement in precisions from using the 
measured instrumental line positions. 
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Figure 2.7.3 Precisions of 1/3 CO 3 rest frame positions determined from at- 
mospheric lines ( v > 2372 cm -1 ) before (line) and after (points) iteration 
procedure. There are improvements of up to 1 x 1 ‘ 0 -5 cm -1 in the line 
positions. 



Figure 2.7.4 Precisions of relative speeds determined from 1/3 COj atmo- 
spheric lines ( v > 2372 cm -1 ) before (line) and after (points) iteration pro- 
cedure. The improvement in the precisions is very small. Note that the 
precision of 5 ms -1 has been maintained to 70 km usingthese lines. The 
precision degreades above 90 km due to the weaknesses of the lines. 
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II. 8 Weighted Averages of Relative Speed Profiles 


We have obtained relative speed profiles from different sets of lines and we 
wish to combine the results to obtain the best profile. The speed profiles can 
be combined by calculating a weighted average. This takes into account the 
varying precisions of the relative speeds with height. Let us assume we have 
several profiles l for each spectrum k given by t>* with a mean uncertainty 
<r(u*). Then the weighted average is 


v k = ^2 tofu* 


(2.36) 


where the weights are given by [Bevington, 1969] 


w i 


> !/*>?) 

El’ 


(2.37) 


The weighted uncertainty is then determined by [Bevington, 1969] 

* 3 (v k ) ~ ^ ^ 2 ' 38 ^ 

Profiles combined in this method axe called “extended” or “merged” rela- 
tive speed profiles and are then used to determine initial estimates of line 
positions for further extension to lower tangent heights. Figures 2.8.1 and 
2.8.2 show the extended relative speeds v k (for sunsets) and their weighted 
precisions a(v k ) determined using this method. It can be seen in Figure 
2.8.2 that a 5 ms -1 precision or better is maintained to about 70 km. In the 
region of overlap between the two sets of 1/3 C0 3 lines between 80 and 100 
km described in the previous two sections we have approximately 30 total 
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lines. The average peak absorptance is about 50% over this range. Using 
a noise level of f n c(N) w = 5 x 10 -8 cm -1 from Table 2.4.1 the predicted 
measurement uncertainty per line from (1.180) is 


5 x 10 'em" 1 4 , 

uncertainty = = 10 cm . 

.50 


(2.39) 


For a central frequency of about 2350 cm 1 for 30 CO* lines the predicted 
uncertainty in the relative speeds is 

10" 4 - 


a(v k ) = 


2350 v/30 


~ 2.5 ms 


-l 


(2.40) 


This is consistent with the observed precisions shown in §11.7. 
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Figure 2.8.1 First extended relative speed profiles v k for the four sunset oc- 
cultations. The speed profiles using 1/3 CO2 absorption lines which did not 
have instrumental counterparts have been combined with the sunset speed 
profiles in Figure 2.6.4. Notation is the same as in Figure 2.6.4. 



Figure 2.8.2 Precisions of first extended relative speed profiles <r(v k ) for all 
occultations. The precision of 5 ms -1 is maintained from 70 to 110 km. 
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II .9 Extensions to Lower Tangent Heights by Using 

H2O and N2O 


In the spirit of the last two sections the profiles were extended by using 
three more bands. Between 40 and 80 km the v 2 HjO atmospheric counter- 
parts to the H a O instrumental lines given in Table 2.1.9 were used for an 
extension. Between 30 and 65 km lines of the 2 1/3 HjO band near 3150 cm -1 
was used and finally between 20 and 40 km lines of the 2v 1 N a O near 2550 
cm -1 was used to complete the profiles. 

The iterative technique was chosen to estimate the positions of the v 2 
water lines instead of relying on the instrumental positions as was done for 
the 1/3 C0 2 lines for a reason which will made clear shortly. In fact, we 
would like to compare the positions retrieved for both sets of lines and this is 
done in the following section. The initial positions were determined from 
the spectra between 50 and 80 km and the relative speed profiles in Figure 
2.8.1. The precisions ^(fj 0 ^) are shown by the bold line in Figure 2.9.1. The 
iterative scheme was then carried out and the precisions of the line positions 
are given by the points in Figure 2.9.1. There is an improvement in the 
mean uncertainties to better than 8 x 10 -5 cm' 1 . Figure 2.9.2 shows the 
initial (connected points) mean uncertainties for the estimated relative speeds 
between 40 and 80 km using these lines and the improved mean uncertainties 
(points). Between 45 and 60 km there is an improvement of up to 3 or 4 
ms -1 . The resultant profiles were then merged with those given in Figure 
2.9.1 and these are shown in Figures 2.9.3. The weighted mean uncertainties 
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are shown in Figure 2.9.4. The mean uncertainty <r(v*) is 5 m/6 or better 
between 50 and 110 km. The relative speeds appear to “flatten out” below 
75 km and do not appear to show the rapid variations with height which are 
noticeable above this level. 

The extension method at this point is fairly straightforward. The next set 
of lines which were found to be useful is the 2i/ 3 H 3 0 vibration band which 
is approximately at twice the frequency as its j/ 3 counterpart. However, the 
sample of nine lines chosen in the overlap regions between 30 and 65 km 
is much smaller than the i/ 3 sample size of more than 20. This does not 
cause significant problems since the higher frequency of these lines helps the 
statistics (2.25). Figure 2.9.5 shows the initial estimates and final estimates 
in the line position precisions (line and points respectively). In every case 
there is a marked improvement. Figure 2.9.6 shows the initial and final mean 
uncertainties in the relative speeds and below 40 km there is an improvement 
of several ms -1 for most occupations. Finally the resultant relative speed 
profiles were merged with Figure 2.9.3 and are shown in Figure 2.9.7. The 
mean uncertainties are shown in Figure 2.9.8. Precisions are now about 5 
m/s in most cases from about 35 to 110 km. 

The final set of lines which were found to be useful is the 2t/i N 3 0 vibra- 
tion band. There were 61 lines available in the overlap regions between 30 and 
40 km. Even though lines were relatively weak in these regions ( A p ~ 10%) 
the large sample size improved the statistics. Figure 2.9.9 shows the initial 
estimates and final estimates in the line position precisions (line and points 
respectively). In every case there is again a significant improvement to less 
than 10 -4 cm -1 . Figure 2.9.10 shows the initial and final mean uncertainties 
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in the relative speeds and below 30 km there is an improvement of several 
ms -1 for most occultations. Finally the resultant relative speed profiles were 
merged with Figure 2.9.8 and are shown in Figure 2.9.11-12 for the sunsets 
and sunrise respectively. The mean uncertainties are shown in Figure 2.9.13. 
Precisions are now between 2 and 5 ms -1 in most cases from about 25 to 110 
km. 

From these results we see that relative speeds precisions of better than 5 
ms -1 can be obtained by measuring rest frame positions from the atmospheric 
lines. The calibration stability of the instrument is adequate to maintain 
these precisions throughout the middle atmosphere. The accuracy of the 
profiles is determined by using the COj instrumental lines to calibrate the 
wavenumber scale. The iteration method of extending the relative speeds to 
lower tangent heights and using new atmospheric bands as they appear in 
the spectra does not degrade the precisions of obtaining line positions and 
speeds. 
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Figure 2.9.1 Precisions of 1/3 H 2 0 rest frame positions determined from at- 
mospheric lines before (line) and after (points) iteration procedure. There is 
a noticable improvement in the line positions of up to 10 -4 cm" 1 . 



spheric lines before (lines) and after (points) iteration procedure. The im- 
provement in precisions is small but up to 1 ms -1 in some cases. The relative 
speed precision between 5 ms -1 and 10 ms -1 is maintained between 45 and 
65 km. This degradation in the precision is due to the small sample of at- 
mospheric 1/3 H3O lines available for initial position estimates. 
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Figure 2.9.3 Second extended relative speed profiles v k for the sunset occul- 
tations. This is obtained by combining the speed profiles determined with 
the atmospheric uj HjO lines with the profiles determined using the u$ CO a 
atmospheric lines. Rapid wind variations do not appear to occur below 70 
km. The separation of the profiles due to changes in the orbital speed is 



Figure 2.9.4 Precisions of second extended relative speed profiles de- 

scribed above for all occupations. A precision of 5 ms -1 is maintained above 
50 km. 




131 



Figure 2.9.5 Precisions of 2 Vi HjO rest frame positions determined from 
atmospheric lines before (line) and after (points) iteration procedure. Preci- 
sions are improved from initial estimates by up to 0.5 xlO -4 cm -1 . 



Figure 2.9.6 Precisions of relative speeds determined from 2t/ a HjO atmo- 
spheric lines before (lines) and after (points) iteration procedure. Precision 
of about 5 ms -1 are maintained betyween 35 and 50 km. Increases above 
that are due to the weaknesses of the lines. Relative speed precisions do not 
improve signifgicantly before and after the iteration procedure. 
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Figure 2.9.7 Third extended relative speed profiles v* for the sunset occul- 
tations. This is obtained by combining the speed profiles determined with 
the atmospheric 2 i/j HjO lines with the profiles determined using the 1/3 CO3 



Figure 2.9.8 Precisions of third extended relative speed profiles <r(v*) for all 
occultations which corresponds to the speed profiles above. A precision of 5 
ms -1 is mainted between about 35 and 110 km. 
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Figure 2.9.9 Precisions of 2i>i NjO rest frame positions determined from 
atmospheric lines before (line) and after (points) iteration procedure. Im- 
provements in precision of more than 1 x 10 -4 cm -1 between the initial and 
final estimates is evident, particularly near the band origin near 2560 cm -1 
and the high J lines. 
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Figure 2.9.10 Precisions of relative speeds determined from 2 i/j NjO atmo- 
spheric lines before (lines) and after (points) iteration procedure. A precision 
of 5 ms -1 is maintained between about 5 and 38 km. Below 25 km line sat- 
uration due to pressure broadening degrades the precisions. 
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Figure 2.9.11 Fourth extended relative speed profiles v k for the sunset oc- 
cultations obtained from combining speed profiles determined using the 2v\ 
NjO absorption lines to the profiles using the 1/3 C0 2 , i / 3 H)jO and 2*/ a H a O 
atmosphwewric absorption lines. The gentle variation in relative speed with 
height below 60 km is evident as well as rapid variations with height above 
70 km. The similarities in the variations with height in the profiles is evi- 
dent. The separations between occultations is due to changes in the orbital 



Figure 2.9.12 Same as in Figure 2.9.11 except for the sunrise occultation. 
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Figure 2.9.13 Precisions of the fourth extended relative speed profiles <r(t>*) 
for all occultations described in Figures 2.9.11-2. A precision of 5 ms -1 or 
better is maintained between 25 and 110 km for all occultations. 



136 


11.10 Line Positions Obtained and Calibration 

Standards 

COMPARING RETRIEVED AND STANDARD LINE POSITIONS 

It has been shown that the line positions determined from the wind re- 
trieval method are precise to O(10 -5 cm -1 ) and accurate to 10 -4 cm -1 . These 
were compared to a set of standard for the C0 2 , H 2 0 and N 2 0 line positions 
compiled in the ATMOS line listing by Brown (1986). The differences be- 
tween retrieved values ^ and the standards recommended by Brown ?Brown 
or in the ATMOS line list Catmos are shown in Figure 2.10.1. For the v 2 H 2 0 
lines there does appear to be a difference between the positions determined 
from the atmospheric lines and the ATMOS standards positions of about 
10~ 4 cm -1 for most of the lines. The other bands appear to be distributed 
more or less about zero. Some of the N 2 0 lines have significant shifts of up 
to 5 x 10 -4 cm -1 . Since the precisions of O(10 -5 cm -1 ) we obtain are of an 
order of magnitude better than this the ATMOS instrument can be useful 
for improving rotation-vibration line wavenumbers and B values §1.6. For 
example, the 2t> 2 H 2 0 lines were only initially known to about 10 -3 cm -1 
accuracy [Brown, 1986] and the precisions we obtain of 5 x 10 -4 cm -1 is 
almost two orders of magnitude better than this. The low precision of the 
ATMOS line list for these lines is reflected in the large overall variation in 
the differences for these lines near 3200 cm -1 . 
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We can estimate the offset a»b between separate line lists ‘a’ and ( b’ 

«7 = (1 + «.b)cj* (2.41) 

where ‘a’ and ( b’ are two different sets of line positions for the same set of 
absorption lines. Table 2.10.1 gives values of a»b between various linelists. 
The line positions determined from the interferometer are compared to ac- 
cepted standards. For the C0 2 and 1/3 HjO lines comparisons are also made 
between the instrumental lines and the positions determined horn the atmo- 
spheric lines. 

We see from the data that the CO3 line positions agree within the uncer- 
tainty of about 3 x 10 -5 cm -1 for the standards and positions determined 
using ATMOS. This is to be expected since we forced the instrumental line 
positions to agree with the standards during the wavenumber calibration. 
Similarly, positions determined using the atmospheric lines agree with the 
positions of the instrumental lines since these positions were derived from 
Doppler shifts, which in turn were derived initially from the instrumental 
positions. The measured positions of the 2 v 3 H3O lines were also found to 
agree with the ATMOS linelisting standards. However, since they are only 
accurate to ~ 10" 3 cm -1 , this indicates only that we are obtaining at least 
that order of accuracy. The average positions of the 2i/j NjO lines were found 
to agree with the ATMOS linelisting and a comparison between the positions 
determined using the atmospheric lines and the standards is given in Figure 
2 . 10 . 2 . 


THE H 2 0 LINE POSITION PROBLEM 
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In the case of the 1/3 H 3 O lines curious discrepancies are apparent. First 
of all, we see that, from comparing the instrumental lines to the standards 
provided by Brown, the difference may be significant but is of the order of the 
measurement erros. Secondly, and most importantly, there is a significant 
discrepancy between the positions determined from the atmospheric H 3 O 
lines and the instrumental H 3 O lines. From the data in Table 2.10.1 we see 
that this shift corresponds to a Doppler shift of almost +15 ms -1 between 
the standards where 

Cinst < (atm* (2.42) 

Due to this discrepancy and the marginal agreement with Brown’s standards 
there is thirdly a significant offset between the atmospheric standards and 
Brown’s standards. This discrepancy poses and interesting problem. If we 
calibrated the instrument with respect to the water lines instead of the car- 
bon dioxide lines then the positions of the water lines determined from the 
atmospheric lines would be in agreement with the instrumental lines, since 
these positions would have been used to determine the Doppler shifts. This 
would cause a shift of —15 ms -1 in the recovered speed profiles. This, when 
determining the positions of the CO 3 lines from the atmospheric lines would 
cause them to be shifted —15 ms -1 from their instrumental counterparts 

Out ^ Catm* (2.43) 

Which calibration is the right one to use? Since the offset is significant 
it implies one of the two components of the gas, namely the C0 3 or H 3 0, 
is moving with respect to the instrument. In the first case when we have 
calibrated relative to the CO 3 standards we find that the H 3 O instrumental 
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frequencies are lower than those determined from the H 3 0 atmospheric lines 
- or synonomously with the C0 3 calibration via Doppler shifts. If we accept 
the accuracy of the positions determined from the £T 3 0 atmospheric lines 
then the negative shift in the H 3 0 instrumental lines indicates that they 
are moving away from the instrument. On the other hand, if we calibrate 
relative to the H 3 0 standards the C0 3 instrumental lines will have a positive 
shift relative to the C0 3 atmospheric line positions. This means that the 
carbon dioxide lines are moving toward the instrument. The first implies a 
source of H 3 0 diffusing into space at a bulk velocity of 15 ms -1 . The second 
implies C0 3 sources converging on the instrument. There is no reasonable 
explanation for this latter kind of source. Therefore we should accept the 
former. This outgassing, or leakage of a trace amount of H 3 0 vapor from the 
shuttle itself, is the diffusion of a trace cloud of water vapor into space. The 
source of this outgassing may be due to water vapor desorbing from the outer 
shell of the shuttle itself. Shaw (1988) pointed out that this may be due to 
the high level of humidity present at the launch site. This also suggests that 
H 3 0 absorption lines may not be useful for absolute line position calibration. 
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Figure 2.10.1 Differences between standard line positions (Brown’s or AT- 
MOS) and retrieved line positions for 1/3 C0 2 , i/ 2 H 2 0, 2 i / 2 H 2 0 and 2v\ N 2 0 
from atmospheric lines. 



Figure 2.10.2 Differences between standard line positions (Brown’s or 
ATMOS) and retrieved line positions for the 2t/i N 2 0 from atmospheric lines. 
This is the same data shown in Figure 2.10.1 but magnified for clarity. Notice 
that the differences in the position appear to systematically vary between the 
band center and the high J lines. These improvements in positions may be 
significant. 
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Figure 2.10.3 Peak absorptances of instrumental 1/3 CO3 lines as observed in 
all five occultations - (2) SR02; (6) SS06; (9) SS09; (11) SSll; (13) SS13. 
The peak absorptances are initially decreasing with time and appear level off 
during SS09. 
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Figure 2.10.4 Peak absorptances of instrumental i/ 2 H 2 0 lines as observed 
in all five occultations using the same labeling in Figure 2.10.3. The peak 
absorptances appear to slightly increase during time and this may be related 
to outgassing of water vapor from orbiter. 




Table 2.10.1 Offsets between various line lists ‘a’ and ‘b*. A - Positions de- 
termined from atmospheric lines. I - Positions determined from instrumental 
lines. LRB - Positions recommended by Brown (1986) for instrumental cali- 
bration. ATMOS - Positions archived in the ATMOS line listing. 


Band 

Linelist ‘a’ 

Linelist ‘b’ 

<*.b (xlO -8 ) 

^(a.b) (xlO -8 ) 

1/3 C0 2 

A 

LRB 

-0.08 

0.49 


I 

LRB 

+0.04 

0.33 


A 

I 

-0.12 

0.50 

2 i/ a HjO 

A 

ATMOS 

-2.83 

3.95 

2 1/\ NjO 

A 

ATMOS 

-0.12 

0.65 

Vj h 2 o 

I 

LRB 

-0.90 

0.71 4 - 


A 

I 

+4.63 

0.90 4 - 


A 

LRB 

+3.73 

0.61 4 - 
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II. 11 Testing the Measurement Theory 


In §1.11 a theory was developed to predict the uncertainties in the mea- 
surement of line positions attainable using the ATMOS spectra. It was shown 
that the measurement uncertainties behave as 


<r(<) = / 


<r(N)W 


(2.44) 


where / is a form factor which depends on the geometrical shapes of the 
lines and the sampling process, w is the width of the line at half maximum, 
A p is the peak absorptance and trij'f) is the background noise level. The 
atmospheric lines were used to test this theory. 

After the relative speeds v k were determined with the methods discussed 
previously, each atmospheric line position c* was shifted to a rest frame 
position c* using 


<? = (l - 7 ) <$■ (2-45) 

Then the difference between this value and the rest frame position estimated 
by (2.29) was calculated for each line 


AiJ = - ij- (2.46) 

These values were then plotted versus A k j. If the measurement theory is a 
realistic one then we would expect the spread of the distribution of Ac* to 
have an envelope width proportional to A” 1 . Figures 2.11.1-4 show this kind 
of behavior for the four bands used in the analysis - v 3 COj, */ 2 H 2 0, 2t/j 
H 2 0, and N 2 0 respectively. Indeed, the uncertainty of Ac* does decrease 
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with increasing A p . Since /, <r(JV ) and W are approximately constant for 
each band we use (2.43) 


A p <t{u) = fa{N) W. . (2.47) 

If we plot the distribution • A versus we would expect to see a 
distribution whose width is independent of the peak absorpt&nce. Figures 
2.11.5-8 show these distributions for the four atmospheric bands used in the 
analysis. Table 2.11.1 list the distribution deviation <r(Ai/ A p ) for each of the 
bands. We see that the COj lines have the smallest line noise §11.4 followed 
by the NjO band and then the H a O bands. This is in agreement with the 
ordering of the total background signal discussed in §11.4. We can also see 
upon comparing the i / 3 COj and HjO noise levels computed from the 
atmospheric lines to those computed from the instrumental lines that the 
atmospheric values are higher. This may be caused by an increase in the 
line noise due to the ratio process. The ratio of the noise levels determined 
from the atmospheric 1/3 CO a lines to the i/ a H a O lines is about 1.6 from 
Table 2.11.1. A similar ratio determined from the values in Tables 2.3. 1-2 is 
about 1.7. These two values are approximately the same which again suggests 
that the measurement precision theory of § 1.11 can be used to describe the 
measurement uncertainties of line positions from the spectra. 
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Figure 2.11.1 A? for the 1/3 C0 2 atmospheric lines as a function of the mea- 
sured peak absorptance Ap. The spread of the data is larger for small peak 
absorptances and becomes sma’ler as the peak absorptance increases. 



Figure 2.11.2 A? for the i/ 3 HjO atmospheric lines as a function of the mea- 
sured peak absorptance A p . The spread of the data is larger for small peak 
absorptances and becomes smaller as the peak absorptance increases. 
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Figure 2.11.3 A<f for the 2i/ 3 HjO atmospheric lines as a function of the 
measured peak absorptance A p . The spread of the data is larger for small 
peak absorptances and becomes smaller as the peak absorptance increases. 






Figure 2.11.4 A? for the 2 v-i N 2 0 atmos’pheric lines as a function of the 
measured peak absorptance A v . The spread of the data is larger for small 
peak absorptances and becomes smaller as the peak absorptance increases. 
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Figure 2.11.5 A p A? for the i/ 3 CO3 atmospheric lines as a function of the 
measured peak absorptance A p . The spread of the data appears to remain 
constant as a function of the peak absorptance in agreement with the mea- 
surement precision theory. 



Figure 2.11.6 ApAc for the u 3 HjO atmospheric lines. The spread of the data 
appears to remain constant as a function of the peak absorptance. 
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Figure 2.11.7 A p Ac for the 2vi H 2 0 atmospheric lines. The spread of the 
data appears to remain constant as a function of the peak absorptance. 



Figure 2.11.8 A p A? for the 2v\ N a O atmospheric lines. The spread of the 
data appears to remain constant as a function of the peak absorptance. 
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Table 2.11.1 Estimates of line noise determined from atmospheric lines. 


Band 

mean (xl 0 - 8 cm _1 ) 

V 3 COj 

6.1 

Vi H ]0 

9.8 

2 i / 2 HjO 

9.5 

2 i/iN a 0 

7.5 
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11.12 Estimating the Doppler Shift due to Orbital 
Motion and the Earth’s Rotation 

THE EPHEMERAL DATA 

The relative speed profiles presented in Figures 2.9.11-12 contain three 
independent Doppler shifts, the projection of the orbital velocity of the shut- 
tle along the line of sight, the rotational speed of the earth along the line of 
sight, and the atmospheric winds along the line of sight. We wish to compute 
the effects due to the first two and then to subtract these from the relative 
speeds profiles to obtain residual wind speed profiles along the line of sight. 
The positions of the orbiter with respect to the earth were provided by God- 
dard Space Flight Center and were tabulated in terms of geocentric distance 
r, longitude <f > , and latitude 0. Positions were evaluated at times t at six 
second intervals for each occultation. The orbital data are given in terms of 
the coordinate system of the earth rotating with the angular speed 

il = - (2.48) 

where T, is the sidereal day of 23h 56m. In order to examine the data in 
terms of a stationary frame it was necessary to translate the longitudes ap- 
propriately. This stationary frame can be arbitrarily chosen for the purposes 
of computing the relative motions of the orbiter with respect to the tangent 
point. I chose an initial time t° = 120d00 : 00 : 00 GMT. The sampled times 
ti for each point i were converted to an elapsed time 


Ti = U — t°. 


(2.49) 
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The longitudes fa were also translated to the reference frame at t° by 

fa = fa + SlTi (2.50) 

so that a longitude fa given at the elapsed time r,- corresponds to the longitude 
fa at the time f°. The positions are given in terms of this geocentric stationary 
coordinate system. It is useful to transform this spherical coordinate system 
into a Cartesian system via 


Xi = Ti COS Oi cos <fa t 

(2.51) 

yi = r \ cos sin fa, 

(2.52) 

Zi = ^ sin Oi. 

(2.53) 


These positions were computed for all of the occupations and axe presented in 
Figures 2.12.1-2 for x-y and x-z planar projections respectively. The figures 
show that the part of the orbit tabulated in the Goddard ephemeris closely 
follow a circular orbit, however there are notable changes in the orbit. Figures 
2.12.3-4 are close up diagrams of the x-y planar projections for the sunset 
and sunrise occupations respectively. There is evidently a regular drift in 
the orbit. This drift is responsible for the shifts in the relative speeds for the 
sunset occupations shown in Figure 2.9.11. The gap between the profiles for 
SS06 and SS09 in Figure 2.9.11, for example, is related the change in position 
of the orbit between these two occupations which is evident in Figure 2.12.4. 

The ephemeris data also give positions of the tangent points in the same 
terms as for the orbital data. Again the longitudes fa “* are translated by 




(2.54) 
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and the Cartesian coordinates are given by 


.‘*n« = r‘“« COS0‘“« cos^* ng , 

(2.55) 

C“* = r‘“* cos 8*— sin#"*, 

(2.56) 

= r‘“* sin 0‘“*. 

(2.57) 


The x-y planar projections of the tangent points as well as orbital data for 
the four sunset occultations is shown in Figure 2.12.5. The orbital drift 
is evident. It is useful to note that the ephemeris data are collected at 
time intervals At ~ 6s. The occultation takes about 5 min which gives us 
n =~ 300s /6s =50 data points to use for the analysis below. 

LEAST SQUARES METHOD 

To get the Doppler shift from the positions of the orbiter and tangent 
points in terms of the Cartesian coordinate system we need to compute the 
rates of change of these positions. We wish to determine v = dr /dr where 
r = r(r) is a function to be determined from the data r,-, r,-. It was found that 
errors in the positions of the orbiter and tangent point caused uncertainties 
in the velocity of over 20 ms -1 between points using simple finite differences 
for the time derivative of the positions. A least squares method was used to 
determine the smooth trajectories of the orbiter and the tangent point as a 
function of elapsed time r. A polynomial of order m was fit to the observed 
data of the form 


*» = (o + ClTi + + ••• + (mT ™ , 

Vi = Vo + VlT{ + TJjT? + ... + rjmT ™ , 

= Co + Cl T » + Ca T i + ••• + CmT - /" • 


(2.58) 

(2.59) 

(2.60) 
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These are three independent sets of equations. If we consider the set of 
equations for the X{ and that there are n independent data points, we can 
write 


x = T£ + e 


(2.61) 


where 


x = 


*1 

X 2 


, T = 


1 Ti T, 


1 Tj T , 2 


1 r n r* 

and the column of model errors is given by 


~fm 



^5 


6 

, * = 

it 

• 


• 

U 


€ = 


Cl 

Cj 


(2.62) 


(2.63) 




In order to get the best set of parameters which best fit the measured 
positions Xi in a least squares sense we need to find the least squares estimator 
£* for the coefficients £ that satisfies 

~[(x-I(f(x-Tf)]=0. (2.64) 

Here, (x-Trf(x-Tr ) is the residual sum of squares. By differentiating 
(2.63) it can be shown that 


(T r T)£* = T r x. 


(2.65) 


So solving for £*, 


( 2 . 66 ) 


r = (T r T) -1 T r x. 

The matrix T r T is an (m + 1) x (m + 1) symmetric matrix. 

The assumptions for the least squares method are that the errors arise 
from distributions which have a mean of zero 

E(e) = 0 (2.67) 

and have a common unknown variance variance <r a , but zero covariance 

E( ee T ) = cr a I. (2.68) 

Under these conditions £* is an unbiased estimator for £ [Myers, 1986] since 
E(x) = T£ and it can be shown that 

E( 0 = t- (2.69) 

The estimator £* also has the minimum variance [Koerts, Abrahamse, 1969]. 

The parameters xt * have been estimated by (2.65) and we wish to esti- 
mate the errors e*. By using (2.60) 

e* = x - T C = Nx (2.70) 

where 

N = I - T(T r T)- 1 T r . (2.71) 

The n x n matrix N has the following properties [Graybill, 1976] 


N = N r , 


(2.72) 
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N = NN = N r N = NN r , (2.74) 

trace N = n — (m + l) = n — m — 1. (2.75) 

The last property can be proven since if A and B are square matrices of 

the same order then trace(A -f B) = trace A + trace B and trace (AB) = 

trace (BA). Then 

trace N = trace [i - T(T r T)" 1 T T ] 

= trace I — trace jT(T r T) _1 T r ] 

= trace I - trace [(T t T)- 1 (T t T)] 

= n - (m + 1) (2.76) 

since T r T is an (m + 1) x (m + 1) matrix. 

We can now find an estimate s 7 for the unknown variance a 7 of the errors. 
One way is to estimate it by the mean square of the residuals 

s 7 = n -1 €* r e* = n -1 (Nx) T Nx = n -1 x r N r Nx = n -1 x r Nx. (2.77) 

This estimate is biased [Myers, 1986] since E{s 7 ) ^ a 7 

E{s 7 ) = n _1 £[x r Nx] = n'^^Ne] (2.78) 

using (2.60) and (2.72). At this point we make use of a theorem [Graybill, 
1976] on the expected value of a quadratic form characterized by a matrix 
Q. Given a random vector e with a mean E(e) = /i and E(ee T ) = <r 7 1 then 

E(e T Qe) = a 7 trace Q + /i r QA*. (2.79) 

In (2.77) we identify fx —* 0 and Q-+Nso that 

E(s 7 ) = o -3 ?! -1 trace N = (n — m — l)n -1 tr 3 . (2.80) 


155 


We see then that even though a 2 is a biased estimate the quantity a* 3 = 
s 3 n/(n — m — 1) is an unbiased estimate. Therefore 

<r' 3 = {n-m- l)" 1 [x - T£*f [x - T£*] (2.81) 

is an unbiased estimate of the variance. Now the variance-covariance matrix 
X for £* can be evaluated by computing X = 25{(£ — £*)(£ — £*) T } H 
has been shown [Myers, 1986] that 

X = E{(( - CM ~ O r } = (2.82) 


Having estimated the coefficients we wish determine the velocity vec- 
tors from the least squares fit. We estimate the position x(r) by the function 


**M = '£C i r i - 

,=o 


If we define the vector 


T = 


then the variance is given by [Myers, 1986] 


(2.83) 


(2.84) 


<r 3 (x*) = r T Xr. 


(2.85) 


The velocity can be estimated from (2.79) by 

<(r) = ^ = ±i f -r>-. 

dT j—0 


( 2 . 86 ) 
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If we define the vector 


T = 


dr 

dr 


(2.87) 


then the variance is estimated by 

ff s (v*) = r r Xr. (2.88) 

The least squares method is similarly applied to the y and z variables as well 
for both the orbiter and tangent point. 

CORRECTION FOR THE EARTH’S ROTATION 

Finally we account for the rotational motion of the earth. At the equator 
the tangential velocity of the earth rotation is 467 ms -1 and causes significant 
Doppler shifts. The velocity vector of the rotation at the tangent point r t “‘ 
is given by 

v‘“« = O X r tan * (2.89) 

where Cl is the angular velocity of the earth’s rotation. When the vector 
functions r*‘“*(r) and r ,orb (r) have been respectively computed and their 
respective velocities v**“* and v* orb (r) determined the relative speed along 
the line of sight can then be estimated. The total contribution to the relative 
speed profiles caused by orbital motion and the earth’s rotation is given by 

Av* = v* orb - v*‘“«. (2.90) 


The unit vector between the orbiter and the tangent is calculated by 

n = (r* orb - r * tBn «) / r * orb -‘“t (2.91) 

where r* orb ~‘“* is the distance between the orbiter and the tangent point. 
Finally the velocity along the line of sight is given by 


Av* = Av* • n. 


(2.92) 
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The quantity Aw* must be subtracted from the relative speed profiles deter- 
mined in §11.9 to obtain residual wind speeds along the line of sight. The 
variance of Aw* can be estimated by expanding (2.91) 

Aw* = Aw* n* -f Aw* n v + Aw* »*. (2.93) 

Since the error in n is negligible and the Aw* are uncorrelated then 

fr 2 (Aw*) = * 3 (vl) n* + fr 2 (w*) nj + 9 2 (v*) h\ (2.94) 

where the variances in the components are determined by (2.87). 

These calculations were carried out for the ephemeris data described in 
Figures 2.12.1-2 for data corresponding to tangent heights from —40 to 300 
km. In each occultation the elapsed time r was determined relative to the 
time when the tangent height was nearest 30 km. This allowed us to obtain 
the best fit for the orbital/rotational speeds with tangent heights between 30 
and 110 km since this altitude range represents the region where the ‘origin’ 
of the polynomial fit is located and the fit is better in this region than at the 
ends near 300 km and -40 km. In Figure 2.12.6-7 the relative speed profiles 
§11.9 together with their orbital/rotational motion contributions are plotted 
for the sunset and sunrise occultations respectively. The slopes of the or- 
bital/rotational motion curves follow the trends of the relative speed profiles 
determined from the atmospheric lines. Figure 2.12.8 shows the uncertainty - 
©■(Aw*) from the least squares analysis using (2.90) for all five occultations. 
The uncertainties in the orbital/rotational speeds were about 3 ms -1 near 
60 km tangent height. SS11 and SS13 have precisions by as much as 5 ms -1 
worse than the other three occultations above 120 km since there are not as 
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many ephemeris data available for these two occultations. However, these 
precisions are still better than 5 ms -1 below 100 km. 

PRECISION OF RESIDUAL SPEED PROFILES 

The residual speed along the line of sight, which is the wind motion, 
is computed by subtracting the orbital-rotational motion from the relative 
speed profiles. The residual speeds for the four sunset and the sunrise occul- 
tations are plotted in Figures 2.12.9-10 respectively. We immediately see the 
striking similarity between the sunset profiles. The measurement precision 
is determined from the precision displayed in Figures 2.12.8 and 2.9.13 

{^(residual) = <7 3 (relative) + er 2 (orbit — rotation). (2.95) 

The precisions for the residual speed profiles are shown in Figure 2.12.11 
This final result shows a total measurement uncertainty of about 5 m/s be- 
tween 25 and 110 km for the five occultations. Since the amount of absorber 
of CO] above 110 km rapidly falls off and the orbital-rotational uncertain- 
ties becomes large with increasing tangent height the uncertainties in the 
residual wind speeds becomes large. Below 20 km the uncertainties becomes 
large due to the saturations of NjO lines as well as orbital-rotation uncertain- 
ties. Between these regions we have obtained the desired total measurement 
precision of 5 ms -1 . 

THE LINE OF SIGHT ORIENTATION 

The orientation of the line of sight relative to the earth’s surface was also 
determined. This calculation was performed by measuring the angle between 

A 

the vector n and the east-west unit vector (f> using standard trigonometry. 
In Figures 2.12.12-13 this angle is plotted as a function of tangent height 
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for the sunset and sunrise occultations respectively. From these data we can 
see the line of sight is inclined to the zonal direction by about 17-18° and 
21-23° for the sunsets and sunrise respectively. Therefore we can conclude 
that residual speeds measured are very nearly the zonal component of the 
wind speeds. In Figure 2.12.9 a positive residual speed is eastward for the 
sunsets while a negative residual speed in Figure 2.12.10 is eastward for 
the sunset. We see that in both figures an eastward prevailing wind of 25 
and 50 m/sec between tangent heights of 20 and 70 km are present in the 
sunsets and sunrise respectively. Above 80 km rapid variations of the zonal 
wind with height axe present. The similarity of zonal wind speed with height 
between the various sunset wind speed profiles in Figure 2.12.9 is particularly 
interesting and various possible causes of the observed wind speed features 
is discussed in §§III. 

The Doppler shifts introduced by the orbital motion of the instrument 
as well as the rotation of the earth have been estimated to 5 ms -1 or better 
between tangent heights of 20 and 120 km from the ephemeris data. The 
ephemeris data can be used to determine absolute zonal wind speeds from the 
relative speed profiles determined in §11.9 without a significant degradation 
of the total precision from 5 ms -1 over the tangent height range. 
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Figure 2.12.1 x — y planar projection of orbital positions computed using the 
ephemeris data showing both sunset (SS) and sunrise (SR) positions for all 
of the occultations. The orientation of the coordinate system was arbitrarily 
chosen, however, the circular shape of the orbit can be seen. The identifica- 
tion of points is not important in this illustration. There are variations in 
the orbital path. 



Figure 2.12.2 Same as in Figure 2.12.2 except for the x — z planar projection. 
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Figure 2.12.3 Close up x — y planar projection of orbital positions for the 
sunrise occultations. There is a noticable changing of the orbit between 
occultation (1) SR01 through (7) SR07. 
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Figure 2.12.4 Close up x — y planar projection of orbital positions for the 
sunset occultations. There is a noticable changing of the orbit between occul- 
tations (1) SS01 through (13) SS13. This changing of the orbital geometry 
causes the shifts between the relative speed profiles of up to 200 ms -1 in 
Figure 2.9.11. 
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Figure 2.12.5 Close up x-y planar projection of orbital and tangent point 
positions for the four sunset occultations - (6) SS06 - (9) SS09 - (11) SSll - 
(13) SS13. The progression of points as time increases during the occultation 
is indicated. This changing of the line of sight geometry causes the shifts 
between the relative speed profiles of up to 200 ms" 1 in Figure 2.9.11. 
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Figure 2.12.6 Orbital-rotational motion plotted with relative speed profiles 
for the sunset occultations - (6) SS06 - (9) SS09 - (11) SS11 - (13) SS13. The 
points on the approximately ‘straight’ curve represent relative speeds due 
to the orbital-rotational motion. The separations of the orbital-rotational 
motion are in agreement with the separation of speed profiles obtained by 
measuring the Doppler shifts of atmospheric lines. The rapid variations above 



Figure 2.12.7 Same as above except for the sunrise occultation - (2) SR02. 
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Figure 2.12.8 Uncertainty in least squares fit for orbital-rotational motion for 
all five occultations. The precision of 5 ms -1 or better is maintained between 
20 and 100 km for all occultations. 



Figure 2.12.9 Residual speed profiles for the sunset occultations - (6) SS06 - 
(9) SS09 - (11) SS11 - (13) SS13. 
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Figure 2.12.10 Residual speed profile for the sunrise occultation - (2) SR02. 



tations. A precision of about 5-6 ms 1 is maintained between 25 and 110 
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Figure 2.13.2 Concentration c(h) for HjO from ATMOS standard at mo- 
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Figure 2.12.12 Angle between the line of sight and the zonal direction for the 
sunset occultations. The angle is approximately 17°.5 between the line of 
sight and the zonal direction indicating that the wind speeds are essentially 



Figure 2.12.13 Angle between the line of sight and the zonal direction for the 
sunrise occultation. The angle is approximately 22° between the line of sight 
and the zonal direction indicating that the wind speeds are essentially the 
zonal component. 
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11.13 Density Weighted Residual Winds 


The wind speeds determined in the previous section were computed un- 
der the assumption that the observations were made solely on gases at the 
tangent points. However, the line of sight penetrates all levels above the 
tangent height. We have seen that there are wind variations with height. If 
these variations are rapid enough, which may be the case for the zonal wind 
profiles determined in the previous section, they may cause a ‘smearing’ of 
the results and degrade the vertical resolution of the profiles. We wish to 
examine how this affects the vertical resolution of the results. 

To begin with we can examine the number density of each gaseous com- 
ponent used to obtain wind speeds, namely CO a , H a O, and N a O. Figures 
2.13.1-3 show the concentrations as a function of height c(k) as given in the 
ATMOS data base. We wish to determine the total fraction of absorber 
present as a function of distance about the tangent point. Figure 2.13.4 
shows the typical geometry of a tangent ray. Refraction in the atmosphere is 
not taken into consideration since this model is only assumed to be a first ap- 
proximation to begin with. The parameter x is the distance from the tangent 
point along the line of sight. Given the tangent height t and the distance x 
the altitude above the surface h t (x) is given by geometry 

h t (x) = \/(r 0 -I- 1) 7 + x 3 - r 0 . (2.96) 

Let the total path length through the atmosphere for a given tangent height 
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t be represented by 

Lt = \J (t*o + Jhnax) 5 ~ (**0 + $)* (2.97) 

where h — is “top” to the atmosphere and can usually be picked to be above 
the region where the wind measurements were carried out. This is a reason- 
able assumption since the amount of infrared absorber decreases rapidly as a 
function of height. We can then determine the number of molecules per unit 
area N t along the line of sight 

N t ~ f dx c(h t (x)). (2.98) 

J—Lt 

Having determined this we want to determine the total fraction ft(x) of 
absorber within a distance z to tangent point along the line of sight 

ft(x) = y J 'dx c{h t {x)). (2.99) 

The integrated fraction was determined for the concentration profiles 
given in Figures 2.13.1-3. A maximum height of 150 km was assumed for 
h mhX . Figures 2.13.5-8 show the total fraction (2.98) as a function of x for 
various values of the tangent height t. It is found that 90% of each species is 
within 400 km of the tangent point, the remainder being distributed about 
the ray. For N 2 O the 90% value threshhold is about 250 km. 

It is of even more importance to estimate the range of altitudes which 
contain certain fractions of the total absorber along the line of sight. Figures 
2.13.5-8 were transformed using (2.95) to examine how the integrated fraction 
is related the h t (x) — t. For a given fraction of the total absorber along the 
line of sight a value of a “layer thickness” ht(x) — t can be determined. The 
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layer thickness for the three components are plotted in Figures 2.13.8-10 for 
various tangent heights t. We see that for COj over the range of tangent 
heights where the atmospheric 1/3 lines are useful that 90% of the absorber 
is contained in the first 10 km or less of the above the tangent point. Since 
the spectra are sampled at approximately 4 km intervals, then 70 and 80% 
of the absorber is contained between a tangent height and the next highest 
sampling height. Between a tangent point and two sampling spectra more 
than 90% of the absorber is present. We therefore expect that the vertical 
resolution of the wind speed at a given tangent height i6 being degraded 
primarily by the next highest level. 

Figure 2.13.11 is a diagram which defines a “layers” model. Basically, we 
can characterize the atmosphere for spectrum k at the tangent height t k by 
the gas lying within ±R/2 of the tangent height where R is the sampling 
resolution of approximately 4 km. We see from the diagram that the layer 
thickness for the tangent spectrum k is R/2 or about 2 km. The layer thick- 
ness for the next level k + 1 is the full resolution R or about 4 km. From 
Figures 2.13.8-10 we see that approximately 60% of each absorber is con- 
tained within the first 2 km. This corresponds to the first dotted line above 
t k in Figure 2.13.11. In the next layer between 2 and 6 km an additional 20 
to 30 % of the absorber is present. In the next layer contains about 10% of 
the absorber, and so on. So we can assume, to a first approximation, that the 
winds measured from Doppler shifts in spectrum k are a weighted average 
of the “true” winds, the weights being the percentage of absorber each layer 
[Shaw, 1985]. 

The above argument suggests the following approximate weighting scheme. 
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Let each measure wind speed at the layer k be given by v k and let the “true” 
speed be given by v' k . Then by the above arguments we have that 

v k ~ 0.6v' fc + 0.2v' fc+l + 0.1t> ,fc+a + ••• (2.100) 

which is a set of equations for all the spectra k and this can be inverted. Fig- 
ures 2.13.12-13 show the wind profiles predicted by this method. Comparing 
these figures to Figures 2.9.9-10 show that in the vicinity of large wind vari- 
ations that the weighted profiles have an enhanced vertical resolution. The 
peaks are shifted slightly in altitude and the magnitudes are larger by up to 
50 m/s near 105 km. 

The layers model can be used to improve the vertical resolution since the 
ATMOS instrument has a limited sampling resolution of about 4 km between 
spectra. The interferences of lines in different layers in the spectra is caused 
by the rate by which the density changes with height and the curvature of 
the earth and is not related to the vertical sampling interval of about 4 km of 
the instrument. If the vertical sampling interval can be decreased in future 
missions it may be possible, however, to determine the contribution of each 
layer to the observed wind profiles. 

Positive speeds in Figures 2.13.12-3 correspond to an eastward motion 
in the northern latitudes and a westward motion in the southern latitudes 
respectively. Figure 2.13.14 shows a zonal average from the four sunset occul- 
tations and Figure 2.13.15 shows the standard deviations about these mean 
values. Between 25 and 65 km there was an eastward wind of about 25 
ms -1 which remained steady within the precision of 5 ms -1 over the period 
of observation and range of longitudes as indicated by the small standard 

'1 
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deviations in Figure 2.13.15. The consistency in the shapes of these profiles 
also confirms the stability of the instrument calibration. The wind motion 
reverses at 65 km and above 80 km. All of the sunset occupations show max- 
ima near 85 and 105 km. The variations in these maximum speeds appear 
to be real. However, if the wind sheets are thin, the variations may be due 
in part to a lack of vertical resolution. These variations in maximum speeds 
may also be due to tides [Chapman and Lindzen, 1970] or breaking gravity 
waves [Fritts, 1984] (see §§III for the definitions of these terms). The minima 
near 90 km correspond to an eastward wind of about 25 ms* 1 . 

The sunrise occupation in Figure 2.13.13 shows a single eastward maxi- 
mum of about 100 ms -1 near 105 km. From 80 to 95 km there is a westward 
motion of about 50 ms -1 which is opposite to the direction observed in the 
northern hemisphere. A wind shear of about 20 ms -1 km -1 occurs between 
84 and 78 km. From 30 to 75 km the eastward motion of about 50 ms -1 is 
nearly twice that observed in the northern hemisphere. 

Because the observations were made over a 24 hour period we can con- 
clude that there were vertically stratified zonal winds at sunset between 25 
and 120 km altitudes near 30°N over a wide range of longitudes. The single 
occupation available for the southern hemisphere also shows winds with mag- 
nitudes, directions, and shears which have similarities to the northern hemi- 
sphere profiles. The measurements were made near the vernal equinox and 
since the asymmetries of the mesospheric circulation between hemispheres 
may be dependent on the seasons it would be useful to have more measure- 
ments at other latitudes and seasons [Labitzke, 1985]. 



Figure 2.13.4 Tangent point geometry of the line of sight. 


if 


x (Km *to' 


Figure 2.13.5 Concentration fraction as a function of x for COj computed 
using data in Figure 2.13.1 for several tangent heights; (a) 60 km, (b) 70 km, 
(c) 80 km, (d) 90 km, (e) 100 km. 
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Figure 2.13.6 Concentration fraction as a function of x for HjO at various 
tangent heights; (a) 30 km, (b) 40 km, (c) 50 km, (d) 60 km, (e) 70 km, (f) 
80 km. 
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Figure 2.13.7 Concentration fraction as a function of x for N 2 O at various 
tangent heights; (a) 20 km, (b) 25 km, (c) 30 km. 




175 



Figure 2.13.8 Concentration fraction as a function of layer thickness for COj 
at various tangent heights; (a) 60 km, (b) 70 km, (c) 80 km, (d) 90 km, (e) 
100 km. 



Figure 2.13.9 Concentration fraction as a function of layer thickness for HjO 
at various tangent heights; (a) 30 km, (b) 40 km, (c) 50 km, (d) 60 km, (e) 
70 km, (f) 80 km. 
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Figure 2.13.12 Density weighted zonal wind speed profiles for the sunset 
occultations. Positive speeds refer to an eastward flow; (6) SS06. (9) SS09, 
(11) SS11, (13) SS13. 



Figure 2.13.13 Density weighted zonal wind speed profiles for the sunrise 
occultation SR02. Negative speeds refer to an eastward flow. 
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Figure 2.13.14 Zonal average of the wind speeds for the four sunset occul- 
tations near 30° N. A positive speed is eastward [van Cleef, Shaw, Farmer, 
1987]. 



Figure 2.13.15 Standard deviations of the zonal average winds shown in Fig- 
ure 2.13.14 [van Cleef, Shaw, Farmer, 1987]. 
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11.14 Conclusion 


Precisions of about 3 x 10~ B cm -1 for the assumed rest frame positions 
of the t/3 CO3 instrumental lines were sufficient to allow for stable frequency 
calibration and wind speeds to be accurate to better than 1 ms -1 . The cali- 
bration scale was found to be stable to better than 2xl0~ 5 cm -1 and there 
were no indications of random calibration offsets between spectra. Thus 
Doppler shifts can be measured to accuracies of better than 5 ms -1 . Analy- 
sis of the orbital ephemeris data indicated that relative speeds between the 
orbiter and the earth with accuracies of 5 ms -1 or better could be obtained. 
We therefore successfully obtained zonal winds to precisions of 5 ms -1 be- 
tween tangent heights of 25 and 120 km. 

The relative positions of the lines used to obtain wind speeds were ob- 
tained to precisions of 5 x 10 -6 cm -1 or better. There appeared to be sig- 
nificant improvements in the precisions of these line positions as compared 
to other sources [Brown, 1986]. The line positions were compared to other 
estimates of rest frame positions and were found to agree with them within 
the known accuracies of about 10 -4 cm -1 (10 -3 cm -1 for the 2vi HjO lines). 
The positions therefore agreed with the uncertainties of the standards. The 
accuracies of the rest frame positions, however, had no effect on the accuracy 
by which Doppler shifts could be measured since instrumental lines of CO3 
were available for absolute calibration. 

The zonal wind speeds obtained by this method showed significant varia- 
tions with height of up to the order of 100 ms -1 . The similarity in the sunset 
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wind profiles, taken at different longitudes but at similar latitudes and local 
times, indicate that the winds may be directly related to the atmospheric dy- 
namics at sunset. Since there was only one sunrise occultation was available 
in the Filter 3 spectra similar conclusions cannot be made at this time. 
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CHAPTER III 


LARGE SCALE DYNAMICS OF THE 
CIRCULATION OF THE EARTH’S 
ATMOSPHERE 


In this chapter we examine several types of activity which are character- 
istic of the large scale circulation of the middle and upper atmosphere. This 
discussion is motivated to understand the ‘oscillating’ nature of the ATMOS 
wind profiles obtained from the spectral data. The different classifications of 
motion discussed in this chapter can be attributed to some of the observed 
features of the zonal winds shown in Figures 2.13-15. It was found that at- 
mospheric tides, planetary scale oscillations produced by the daily or diurnal 
ultraviolet absorption and subsequent heating of the atmosphere by bound- 
ary layer water vapor and stratospheric ozone, are significantly large enough 
to produce the kind of wind speeds measured during the mission, particu- 
larly above stratospheric heights. In the stratosphere the wind motion is 
dominated by the thermal wind which can produce zonal wind motions of 
up to several tens of ms -1 and these appear in the ATMOS wind profiles at 
approximately the correct magnitude predicted [Salby, 1981] of several tens 
of ms -1 . 
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This chapter presents the linear theory of global scale atmospheric dis- 
turbances. The linearized theory predicts that the amplitudes of these dis- 
turbances, in terms of velocity, pressure aijd temperature perturbations from 
the mean state, will grow with height proportional to p where p\, is the 
mean background density. The exponential increase with height is due to en- 
ergy transporting through the exponentially decreasing background density 
imposed by this linear theory [Kato, 1980]. The linear theory presented in 
this chapter includes radiative damping or infrared cooling [Andrews, Holton, 
Leovy, 1987] primarily due to the radiative properties of the 15 pm j/ 2 C0 2 
vibration band [Dickinson, 1984] but neglects dissipative mechanisms such as 
the kinematic viscosity which varies as p [Andrews, Holton, Leovy, 1987; 
Pedlosky, 1987]. 

In the previous chapter we showed how zonal wind speeds between 25 and 
120 km were obtained with a precision of 5 ms' 1 from the ATMOS spectra. 
The results indicated that there were variations of about 100 ms -1 in the wind 
speed over changes of only a few kilometers in the tangent height. The profiles 
in Figures 2.13.12-13 show that the amplitudes of the wind variations become 
larger with increasing height. The variations appear to begin at about 70 km 
with an amplitude of about 20 ms -1 and increase to almost 150 ms -1 at about 
110 km. The sunset measurements occurred at a variety of longitudes at the 
same latitude and the wind structure appeared to be preserved between the 
occultations. This regularity of the wind patterns suggests that they are 
caused by global scale dynamics. Wind measurements by other researchers 
[Reed, 1969; Manson, Meek, 1986; Wallace, Tadd, 1974; Muraoka, Sugiyama, 
kawahira, 1988 - for example] indicate that the behavior observed by ATMOS 
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is typical of mesospheric and lower thermospheric dynamics. 

The atmosphere heats and cools on a regular daily basis and this forces the 
upper atmospheric disturbances are periodic in nature. The classical or lin- 
ear theory of atmospheric tides first described by Siebert (1961) and later by 
Chapman and Lindzen (1970) and Kato (1980) predicts atmospheric distur- 
bances of pressure, density, temperature and wind due to atmospheric heat- 
ing which propagate to thermospheric heights. The exponentially decreasing 
density of the atmosphere with height allows the amplitudes of pressure, 
density, temperature and velocity perturbations waves to grow with height 
[Fritts, 1984]. The classical theory of tides [Siebert, 1961; Chapman and 
Lindzen, 1970] predicts unlimited growth in wave amplitudes with height. 
This is due to the linearity of the theory. However, instabilities from the 
nonlinear terms in the fluid dynamical equations equations cause the waves 
to break which can lead to turbulence. Adiabatic instabilities due to large 
local potential temperature perturbations also may contribute to the break- 
down of the tidal waves [Lindzen ,1981]. It is desirable to see if we can 
propose a theory which can ‘explain’ the nature of the wind speeds obtained 
by ATMOS. 

In order to test the linear theory a model of atmospheric heating rates 
proposed by Groves (1982) was used as a periodic heating term and the 
response of the atmosphere is examined. The linear theory predicts that 
the leading diurnal tidal mode may be responsible for the regular structure 
observed in the ATMOS results. However, since the linear theory does not 
‘break’, this response is overestimated. 
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III.l The Atmospheric Temperature Distribution 

The earth’s atmosphere has an average altitude dependent temperature 
distribution similar to the one given in Figure 3.1.1 [from Andrews, Holton, 
Leovy, 1987]. Near the surface of the earth convection and adiabatic expan- 
sion of the atmosphere cause the temperature of the atmosphere to decrease 
with height [Kato, 1980]. The temperature minimum near 12 km is caused 
primarily by infrared emission of HjO vapor. The region of the temperature 
decrease is known as the troposphere and the region where the minimum 
occurs is known as the tropopause. Above the tropopause the atmospheric 
temperature rises to a maximum near 50 km. The region between about 10 
and 50 km where the temperature rises is known as the stratosphere and this 
is due mainly to the absorption of solar ultraviolet radiation by O3 [Andrews, 
Holton, Leovy, 1987]. The maximum near 50 km is called the stratopause 
and occurs because the O3 concentration decreases rapidly above this level 
and reduces the solar energy absorption. Above the stratosphere the tem- 
perature decreases again primarily due to infrared radiative cooling of CO] 
[Kato, 1980]. The temperature continues to decrease to a minimum near 85 
km where the temperature again again begins to increase due to radiative so- 
lar heating by dissociation of molecular oxygen and by ionization of O [Kato, 
1980]. The region between 50 and 85 km where the temperature decreases 
is known as the mesosphere and the region of the temperature minimum at 
about 85 km is known as the mesopause. The region above 85 km where the 
temperature increases is known as the thermosphere. Above the mesopause 
the amount of ionization of atomic species increases with height and the 
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atmosphere can conduct electricity and interact with magnetic fields. This 
region is known as the ionosphere and its electromagnetic interactions play 
a dominant role in the atmosphere above 150 km [Volland, Myers, 1977]. 



Figure 3.1.1 Midlatitude temperature profile based on U.S. Standard Atmo- 
sphere (1976) [from Andrews, Holton, Leovy, 1987]. 
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III. 2 Fluid Equations of Motion 


The atmosphere’s circulation is governed by the Navier-Stokes equation 
of fluid dynamics which can be expressed, as seen in a rotating coordinate 
frame, by [Pedlosky, 1987] 


Du 

Dt 


-f 2fl x u J = -VP + p + T 


(3.1) 


where u is the velocity as observed in the rotating frame, the mass density 
of the fluid is defined by p, the pressure is given by P, the constant angular 
velocity of the earth’s rotation is ft = 7.272 X 10~ 8 s -1 , the potential of 
conservative forces is $, and other non-conservative forces such as friction 
are represented by the term T. The factor 2fi x u, the Coriolis acceleration, 
comes about from the coordinate transformation to the rotating frame of the 
earth. The total time derivative D/Dt is given by 


D 8 

— = (- u • V 

Dt dt 


(3.2) 


and represents the total time rate of change of any property following a fluid 
element. The conservation of mass is governed by the continuity equation 
and this is commonly expressed as [Pedlosky, 1987]. 




(3.3) 
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For the atmosphere, the ideal gas law is used for the equation of state. 
This is given by the relation 

P = 4jRT (3.4) 

M 

where R is the gas constant, T is the kinetic temperature and M is the average 
molecular weight. The energy balance of the atmosphere is described by the 
thermodynamic equation [Pedlosky, 1987] 

m _ m_± h 

Dt Cp T 

where the potential temperature d is defined as 

/ pO\ R / C P 

* = T (^p J . (3.6) 

Here P° is a constant reference pressure at temperature T° and 7 i is the 
heating rate per unit mass. Since Cp — R + Cy where Cy and Cp are 
the specific heats at constant volume and pressure respectively, then k = 
R/Cp = 1 — 7 -1 where 7 = Cp/Cy = 1.4 is the ratio of specific heats 
[Pedlosky, 1987]. The five relations given in (3.1, 3, 5) and the equation of 
state (3.4) are used to describe atmospheric dynamics at all time and length 
scales. 

The equations of motions (3.1, 3, 5) can be transformed into spherical 
coordinates representing the height z above the earth with radius a, latitude 
9 and longitude <j>. We denote the eastward or zonal component of the velocity 
by u, the northward or meridonial component of the velocity by v, and the 
upward or vertical component of velocity by w. The zonal, meridonial and 
vertical momentum equations are [Pedlosky, 1987] 
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Du wu vu 

— — i tan 6 — 2(1 sin 6 v + 20, cos 6 w 

Dt a a 

= 1 — + 

pa cos 6 d<f> p ’ 

Du tun u J 1 9P J~t 

~ H 1 tan 0 + 2fl sin 0 u = — H , 

Dt a a pa da p 

Dw u 2 +v 2 13P P, 

Dt a p oz p 


(3.7) 

(3.8) 

(3.9) 


respectively where the gravitational acceleration is given in (3.1) by V$ = 
— gz ~ —9.8 ms -1 z. In this coordinate system the continuity equation (3.3) 
is 


Dp 


Dt 


+ P 


f dw 2w 1 ( d . . du \ ) , . 

\ 02 a a cos 6 \d6' ' d<f>) J v ' 

These equations (3.7-10) along with the thermodynamical relationship (3.5) 

(3.11) 


m _ M_±nj 

Dt~ CpT U 


and the spherical coordinate representation of the time derivative operator 
(3.2) 

7^ /} /) «t 41 

(3.12) 


D d d v d u d 

= J- yj 1 1 

Dt dt dz add a cos 0 d<t> 


are used to describe the dynamics in the spherical geometry of the earth. 


189 


III. 3 Hydrostatic State and the Thermal Wind 

The simplest static solution to the dynamical equations (3.7-11) is the 
time independent hydrostatic case where the atmospheric velocity vanishes 
everywhere u = 0. This state can only occur if there are no horizontal tem- 
perature gradients. Assuming that this state represents the background is 
what is known as the hydrostatic approximation. If we assume that non- 
conservative forces vanish then from (3.7-11) the conditions for this solution 
are 


dPa 

9<f> 

= o, 

(3.13) 

dPi 

90 

= o, 

(3.14) 

dl\ 

dz 

= ~Ph9, 

(3.15) 

n 

= 0. 

(3.16) 


The hydrostatic pressure field (denoted by the ‘h’ subscript) is independent 
of longitude (j> and latitude 0. The height dependence of the hydrostatic with 
height z varies according to (3.15). By this relationship and the equation of 
state (3.4) the hydrostatic density and temperature (similar to Figure 3.1.1) 
must also only depend on the height z and the heating "H must also vanish 
for the hydrostatic state to be realized. 

The hydrostatic case is a highly idealized approximation to the true at- 
mospheric state. Figure 3. 3. 1-2 show the zonal mean temperatures for alti- 
tudes of up to approximately 160 km respectively during equinox and solstice 
conditions [from Salby, 1981]. The average over latitudes is similar to the 
temperature profile in Figure 3.1.1 and there Me variations between seasons 
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as well. From these data we see there are significant deviations of up to 
10% in the temperature at a given height between different latitudes for each 
season and that the zonally averaged temperature background T\, is not in- 
dependent of latitude 9. This temperature variation violates the conditions 
for the hydrostatic state which implies that the velocity field will not vanish 
and the horizontal temperature gradients induce a zonal velocity field. This 
is called the thermal wind and by assuming this to be the background state 
is known as the thermal wind approximation. 

The simplest case for the thermal wind is a steady state [6/dt — 0) 
zonally symmetric solution (d/d<j> = 0) from the temperature data. The 
zonal mean wind can be derived from the latitudinal and vertical tempera- 
ture gradients and assumes that the vertical and meridonial velocities vanish 
[Andrew, Holton, Leovy, 1987] v — w — 0. In the absence of dissipation, the 
dynamical equations become (3.7-11) 


dpb 

6<f> 

= o, 

(3.17) 

tan 0 + 2U sin 0 Ub 

1 ap *> 
Pba 60 ' 

(3.18) 

ul 

— — 23T2 cos 9 «b 

a 

1 0Fb 

Q 9t 

pb oz 

(3.19) 

du\, 

d<f> 

= o, 

(3.20) 

n 

= 0. 

(3.21) 


These equations show that the fields are zonally symmetric and exist in the 
absence of heating. A differential equation which can describe the zonal 
velocity in terms of the background temperature T\> can be derived from 
these equations. By combining (3.18-9) and using the equation of state it can 
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be shown [Kato, 1980] that to a first (linearized) approximation the zonal 
mean winds Ub can be derived from the zonal mean temperature Tb by the 
thermal wind equation 


« (ln r b ) Uk = __X_^ (lnTb ). 


(3.22) 


The differential equation which describes the relationship between the pres- 
sure and density can be determined by combining (3.18-9) 


dPb cot 0 dPb 


dz 


+ 


39 


= ~Pb9- 


(3.23) 


The mean zonal wind determined from the mean temperature profiles 
in Figures 3.3. 1-2 are given in 3.3.3-4 for equinox and solstice conditions 
respectively [after Salby, 1981]. Measurements of zonal mean temperatures 
are routinely conducted during the various month’s of the year [see Barnett, 
Carney, 1985]. The zonal mean winds change as a function of height and 
latitude with the seasons. An eastward jet of about 50 ms -1 develops in the 
southern hemisphere during the equinox near 50° S around 9 scale heights 
or 60 km and continues to increase to about 80 ms -1 during the southern 
winter. The eastward maximum zonal wind of about 70 ms -1 between 40 
and 80 km in sunrise zonal wind profile SR02 in Figure 2.13.13 determined 
from the ATMOS experiment can be explained by this jet near 50° S since 
the mission was flown about one month after the vernal equinox. The sunset 
wind speeds SS06, SS09, SS11 and SS13 in Figure 2.13.12 show an eastward 
flow of about 20 ms -1 near 50 km. This is about 10 ms -1 less than the value 
predicted for the equinox conditions of about 30 ms -1 . However, between 
the equinox and summer solstice the zonal wind speed will decrease with 
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time and changes direction which is evident by comparing the zonal winds 
in Figures 3.3.3-4 at 30° N. 

The thermal wind due to latitudinal gradients of the zonal mean tem- 
perature can explain the eastward wind below 70 km which was observed 
observed in the ATMOS wind data. The consistency the zonal winds fea- 
tures of the sunset wind profiles in Figure 2.13.12 can be explained by the 
zonal symmetry of the thermal winds. Thermal wind features below 70 km 
can be expected in future ATMOS missions. Seasonal as well as latitudinal 
variability in the zonal winds should be expected [Barnett, Corney, 1985]. 
The ATMOS wind profiles confirm the steadiness of the zonal wind over a 
wide range of longitudes. The deviations of the averages of the four sunset 
profiles in Figure 2.13.15 are between 5 and 10 ms -1 below 60 km. This ob- 
servation may be useful for calibrating spectra in which instrumental lines are 
not available for calibration §11.1 since bith observation and theory indicate 
that the zonal mean winds are very steady and predictable. By assuming 
that the thermal wind is steady, it may be useful in future missions to ad- 
just the wavenumber calibration of the spectra to obtain the best mutual 
correlations with wind speeds below 60 km. 

The temperature profiles Tb in Figures 3.3. 1-2 and the corresponding 
zonal winds Ub in Figures 3.3.3-4 show that the atmospheric state is highly 
non-uniform with latitude and height. Above 13 scale heights or about 85 km 
the thermal wind vanishes since the zonal mean temperature in thermosphere 
effectively becomes independent of latitude [Salby, 1981] which can be seen 
in Figures 3.3. 1-2. However, the ATMOS measurements indicate that there 
exist winds with zonal speeds of over 100 ms -1 and large wind shears. This 
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cannot be explained by the thermal wind. We can show that time dependent 
global waves play an important role in the circulation of the middle and up- 
per atmosphere. From the perspective of modeling atmospheric circulation, 
the mean zonal-wind and temperature profiles can serve as a background 
upon which time dependent fluctuations occur [Lindzen, Hong, 1974; Wal- 
tersheid, Venkateswaran, 1979]. Rossby waves, which are homogeneous free 
atmospheric oscillations having periods on the order of days, can interact 
with the zonal-mean circulation and modify the zonal-mean winds predicted 
by the thermal wind equation (3.22) [Dickinson, 1969]. Atmospheric tides 
[Chapman and Lindzen, 1970], forced by the periodic heating and cooling of 
atmosphere by absorption of ultraviolet solar radiation by ozone and water 
vapor, have periods of a day and fractions thereof and interact with Rossby 
waves and the zonal mean flow to produce net accelerations as well [Teit- 
elbaum, Vial, 1981; Zurek, 1986]. The mean zonal winds conversely affect 
the structure of the planetary waves and tides [see Lindzen and Hong, 1974; 
Waltersheid, Venkateswaran, 1976] and these types of atmospheric activity 
may be highly interdependent. 

This section has shown that the eastward winds of about 30 ms -1 and 
70 ms -1 in the northern and southern hemisphere ATMOS wind profiles 
below about 70 km can be described by the regular thermal wind associated 
with latitudinal temperature gradients of the lower and middle atmosphere. 
However, these represent are slowly varying with a time scale of the order of 
the seasons. Other time dependent effects also occur the wave like character 
of the ATMOS wind profiles §11.13 may be caused by other dynamics. These 
are discussed in the next few sections. 
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(o) TEMPERATURE <*) 



Figures 3.3.1 Zonal mean temperatures for equinox conditions [from Salby, 
1981]. 



LATITUDE 


Figures 3.3.2 Zonal mean temperatures for winter soltice conditions [from 
Salby, 1981]. Summer soltice conditions are approximately the same with by 
interchanging N and S. 
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Figures 3.3.3 Zonal mean winds for equinox conditions. Solid contours are 
positive (eastward) speeds while dashed contours are negative (westward) 
speeds [from Salby, 1981]. 



Figures 3.3.4 Zonal mean winds for winter soltice conditions. Solid contours 
are positive (eastward) speeds while dashed contours are negative (westward) 
speeds [from Salby, 1981]. Summer soltice conditions are approximately the 
same with by interchanging N and S. 
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III. 4 Linear Waves on a Hydrostatic Spherical 

Atmosphere 

We have seen that winds occur in the upper atmosphere. Indeed, the 
atmosphere can support a variety of free and forced wave motions. In the next 
few sections we briefly review the linear wave theory of global scale dynamics. 
The waves which have meteorological significance range from acoustic gravity 
waves [Kato, 1980] which are free oscillations with periods of several hours or 
less due to stratification of the background temperature, atmospheric tides 
[Chapman and Lindzen, 1970] with periods of a day or less and planetary 
or Rossby waves with periods of several days or more [Andrews, Holton, 
Leovy, 1987]. Of these three types of waves the latter two, atmospheric 
tides and Rossby waves, are significantly affected by the variation of Coriolis 
acceleration with latitude 6. I will not examine acoustic gravity waves in 
this report due their limited horizontal length scale [Gossard, Hooke, 1975; 
Kato, 1980]. However, they can be important in affecting the zonal mean 
circulations due to wave breaking [Fritts, 1986] and subsequent transfer of 
momentum due to turbulent eddy diffusion [Andrews, Holton, Leovy, 1987 ]. 
We will examine, therefore, waves which are of the greatest importance for 
the large scale circulation of the middle atmosphere since it is not possible 
to examine all waves in detail. 

We shall only examine waves in an atmospheric background state which 
is vertically stratified without horizontal temperature gradients. This cor- 
responds to the hydrostatic approximation. It was mentioned that the ac- 
tual atmosphere is not in a hydrostatic state and possesses a mean thermal 
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wind. The presence of background winds can affects the way atmospheric 
waves propagate through the atmosphere [Lindzen, Hong, 1974; Waltersheid, 
Venkateswaran, 1979]. Although background winds tend to complicate lin- 
ear wave calculations, qualitatively the results are similar to those obtained 
using the more tractable hydrostatic background approximation and can be 
found in the literature [Lindzen, Hong, 1974; Waltersheid, Venkateswaran, 
1979]. 

The system of linear equations describing fluid motions on a sphere was 
originated by Laplace (1799, 1825) and is useful for understanding atmo- 
spheric tides and global modes (Rossby waves). The purpose of the following 
discussion is to examine some of the most important solutions. 

The nonlinear equations which describe global scale dynamics are de- 
scribed in Appendix A. The final equations (A. 61-7) are presented in dimen- 
sionless form and represent the behavior of the system in terms of the length 
and time scales for global scale disturbances. The scaling of the variables 
shows which terms in the dynamical equations are of primary and secondary 
importance. These derivations, although extensive, are the underlying basis 
for the linear wave theory of disturbances on a vertically stratified rotating 
spherical atmosphere. The description of linear wave theory which follows is 
based on the scaling derivations. 

Linear atmospheric wave theory is derived from the equations of motion 
(A.61-7) by ignoring the nonlinear terms to a first approximation. The lin- 
ear balance equations can then be written. The continuity equation (A.62) 
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becomes 


£ 4 + a 


dw ( dA \ 
dx ~ [ A ~ ~di) 


dv' 1 du' 
w + — + z r— = 0. 


dfi 1 — fi 2 dtp 
The horizontal momentum equations (A.63-4) are 

ap 


du ' 
1 dt 
dv' 


‘-di-^ = ~W 


uv , I f-, 2\ &P 

aT + ' ,u = -l 1 -'*) a?> 


the vertical momentum equation (A. 65) is 




dP 

dx’ 


the thermodynamical relation (A.66) is 


(3.24) 


(3.25) 

(3.26) 


(3.27) 


30 

'.F- gi +>Pw = M, 


(3.28) 


and the relationship between the potential temperature, pressure and density 
perturbations (A.67) is 

0 = -AP - p. (3.29) 

7 

The derivation of these equations is described in Appendix A. Basically, in 
dimensionless format, the zonal, meridonial and vertical velocity components 
are represented by u', v', and w respectively while the pressure, density and 
potential temperature perturbations about the hydrostatic state are repre- 
sented by P, p and 0 respectively. The vertical coordinate x defined in 
(A.54) represents the number of atmospheric e-folding scale heights above 
the surface, and a latitude coordinate has been introduced fi = sin 8. The 
hydrostatic temperature distribution is contained in the function A (A.22) 
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and the dimensionless local bouyancy frequency M (A. 49; see Appendix B). 
Time dependent external heating is represented by the dimensionless func- 
tion (A.52). This linearized approximation is used to describe atmospheric 
dynamics provided that the (dimensionless) variables have magnitudes much 
smaller than unity. 

The standard procedure [Siebert, 1961; Chapman, Lindzen, 1970; Salby, 
1979; Kato, 1980; Andrews, Holton, Leovy, 1987] for determining the solu- 
tions to the set of equations (3.24-9) is to pose the form 

(u', v ', P, w, p, H) = 3?{[ti'(0,a:), v'(0,x), P(0,x ), w'(6,x), p(9,x), 'H(B,x)] 

x exp[i(<rt + a^)]} . (3.30) 

The zonal wave number s must be an integral value due to the periodicity of 
the azimuthal coordinate <f>. The time t is the dimensionless representation 
described in the Appendix A, the rotational period of the earth being 2x in 
this representation. The (dimensionless) frequency number «r is continuous 
in general. 

We define the number f = ae = <r/2 for convenience. The periodic 
horizontal momentum equations (3.25-6) are then 

ifv! — p. v 1 = —isP, (3.31) 

ifv' + p u' = - (l - /i 2 ) (3.32) 

By using these equations the horizontal velocity field (u',v') can be deter- 
mined in terms of the pressure P 



(3.33) 
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= jrzj> f ( J -»')%- ‘>‘ p ■ ( 3 - 34 > 

It is important to note that for / < 1 which corresponds to atmospheric 
tidal distrubances [Chapman, Lindzen, 1970] there are two latitudes where 
fi = ±f and this causes the denominators in (3.33) and (3.34) to vanish. The 
latitudes, which depend on the frequency of the wavesthrough /, where this 
occurs are called the critical latitudes. Brilluon (1932) first showed that the 
numerators in (3.33-4) also vanish at these latitudes and they do not impose 
any theoretical difficulties. For the diurnal tide §111.6 the critical latitude is 
at 30° N and S. 


Following Siebert (1961) we define the dimensionless divergence 

ir r-f . & v> I 

X = V • u — A—— + — — |- -u . 

ox d/i 1 - fl a 


(3.35) 


Since u' and v' depend on the pressure P the horizontal part of the 
divergence (3.35) is 


dv' is , 

X h = — + - u = if CP 

Ofl 1 — fl* 


(3.36) 


where the operator C is defined as 


by substituting (3.33-4) into (3.36). The relationship (3.35) can be written 
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ifFp 


-(a- 


dA\ 

dx J 


I w + X = 0, 


and 


if F$ -f /flw = "H 


(3.40) 


(3.41) 


respectively. From (3.29) we can solve for the density p and substitute into 
(3.39-40) to obtain the vertical momentum and continuity equations 


i> = -tPP + A~, 


-ifF\P - if Ft 
7 




w + X = 0. 


(3.42) 


(3.43) 


The continuity equation (3.43) and thermodynamic equation (3.41) can be 
added to yield 

ifFP = w - 1 {X - H ) . (3.44) 

This equation can be differentiated with respect to z (A.54) to give 


dw dP 


d 1 dA 
dx A dx 


(X-H). 


(3.45) 


The vertical momentum equation (3.42) can be rearranged 

BP 

A^— = + /PP. 

dx 


(3.46) 


We now solve for •& in (3.41) and substitute this and (3.44) into (3.46) to 
show 

BP 


This result can then be substituted into (3.45) to give 


(3.47) 


. dw 


dx 

dx 


—(£-«)] 


(3.48) 
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Differentiating this result with respect to z gives 


t d dw k \d 3 x ax fa 3 n a?Al 

A dx A dx ~ lA [a® 2 K dx v&r 2 dx)\' 

If we now differentiate (3.38) with respect to z 


k dX k d k dw 

A— A— A— — 

dx dx dx 


•fr\ dP 


(3.49) 

(3.50) 


We now substitute (3.47) and (3.49) into this equation to obtain the following 
differential equation for the dimensionless divergence X 

a*x ax jam m\ _ i (, 

dx 1 Sx dx ) FA‘ C Y^ X ( ix)\ ' 1 

This equation can be further simplified by introducing an auxiliary function 
y, which we will also term the vertical structure function , as 


X = ye*' 3 + H (3.52) 


Equation (3.51) then transforms to 

= W- 7 ’]. (3.53) 

This equation can be solved by using the method of separation of variables 
as described by Siebert (1961) and Chapman and Lindzen (1970). Consider 
a set of functions {©£*(/*)} for n = 0,1,... which is assumed to be complete 
over the interval — 1 < fi < 1 for each chosen pair cr, s. Then the vertical 
structure function y and heating W can be expanded on this complete set by 

y = £y»(*) 0 nV), ( 3 - 54 ) 

n 

w = E«"(*)0"W- 


(3.55) 



203 


By calling the separation constant —A'* for each n component the following 
pair of equations can be derived 


£0- + A- 0- = 0, 


+ * n Kp -- v = 1 
+ [FA a 4 Vn 7 FA n ' 


(3.57) 


The eigenvalue problem for the operator £ in (3.56) is known as Laplace’s 
Tidal Equation [Longuett- Higgins, 1968; Chapman and Lindzen, 1970] since 
it was originally derived by Laplace (1799, 1825) to describe the free oscil- 
lations of an ocean. The eigenfunctions 0**, commonly called Hough func- 
tions (after Hough, 1898), form a complete orthogonal set for each set <r, 
s. A method for computing the eigenfunctions and eigenvalues described 
in Appendix C is based on a method similar to one described by Longuett- 
Higgins (1968). The orthogonality of the Hough functions is also proved. 
The equation (3.57) is called the vertical structure equation [Siebert, 1961]. 
The solution for y n is determined by the eigenvalue A£* and the heating term 
‘ H together with appropriate boundary conditions. 

The eigenvalues of (3.56) are related to the frequency of the mode through 
<r and, for each azimuthal mode number s, a relationship between the eigen- 


values A and <r can be determined. This relationship for the values of 
— 1 < s < +2, as computed by Longuett-Higgins (1968), is shown in Figures 
3.4. 1-4. The modes where a > 0 (s < 0) are called westward ( eastward ) 
migrating modes since the longitude <f> of constant phase decrease (increase) 
as time t increases. Negative values of the eigenvalue A exist only for modes 
with <r < 1 [Longuett-Higgins, 1968; Andrews, Leovy, Holton, 1987]. 
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Once the vertical structure equation (3.57) has been solved for each eigen- 
mode n we can determine the physical fields. In terms of the Hough function 
expansion, equation (3.38) can be written 


X n -A 


dw n 

dx 


-ifK'Pn 


(3.58) 


and from equation (3.48) 


A 


dwn 

dx 


= 7 


dX n 

dx 


-nX n 



(3.59) 


By using (3.52) these can be combined to give the pressure in terms of the 
vertical structure function y n 


Pn = 


ifK 1 



(3.60) 


With this result and equation (3.44) the vertical velocity becomes 



Combining (3.60) with (3.59) the density is 


(3.61) 


Pn = 


•7*S 


-Ae*' 2 


\ <Py n _ (. _ 1 dA\ dy n (l 1 dA\ j 

X dx 3 \ A dx J dx \4 2A dx J ^ n J 

The velocity fields u' and v' are determined from equations (3.33-4) 


(3.62) 


< = Pn(x)U'(n), (3.63) 

< = iP n {x) (3.64) 
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where 


K 

yl<r* 




(3.65) 


(3.66) 


The vertical components of the physical fields (3.60-4) along with the 
Hough eigenmodes ©„ and the horizontal velocity modes U' n and describe 
the time and space solutions of the linear wave equations for fluid motion on 
a spherical atmosphere. 

The linear theory serves as tool by which periodic atmospheric distur- 
bances may be studied. In the sections that follow we will show that both 
homogeneous solutions and nonhomogenous solutions of the vertical struc- 
ture equation (3.57) exist and that planetary waves and atmospheric tides 
are representative of these kinds of solutions respectively. It will be shown 
that the homogeneous solutions, the planetary waves, are not likely to be 
responsible for the upper atmospheric winds observed using the ATMOS in- 
strument. However, the nonhomogeneous solutions give rise to perturbations 
from the mean state large enough to account for the global extent and the 
magnitudes of the winds observed above 70 km by ATMOS. 
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Figure 3.4.1 Eigenfrequencies of eastward traveling modes when a = — 1 
[adapted from Longuet- Higgins, 1968]. The frequency number a (the number 
of oscillations per day) is plotted against the eigenvalue A for various modes. 




Figure 3.4.2 Eigenfrequencies <r/2 of stationary modes when s = 0 [adapted 
from Longuet-Higgins, 1968]. The frequency number <r (the number of oscil- 
lations per day) is plotted against the eigenvalue A for various modes. 


V>lc* 



[adapted from Longuet-Higgins, 1968]. The frequency number a (the number 
of oscillations per day) is plotted against the eigenvalue A for various modes. 
The horizontal line at <r = 2 intercepts the dispersion curves which gives 
the diurnal eigenvalues A„. The Rossby 5-day wave © and 16-day wave © 
are marked on the intersection of the dispersion curves with the atmosphere 
eigenvalue Aq = 8.88 giving the dimensionless frequencies <r of oscillation. 



[adapted from Longuet-Higgins, 1968]. The frequency number <r (the number 
of oscillations per day) is plotted against the eigenvalue A for various modes. 
The horizontal line at a = 2 intercepts the dispersion curves which gives the 
semi-diurnal eigenvalues A n . The leading Hough mode is marked 
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III. 5 Free Planetary Waves - The Homogeneous 
Solutions to the Linear Wave Equations 

The linear wave theory described in the previous section is used to de- 
scribe both free and forced oscillations of the atmosphere. One set of solutions 
to the linear wave equations (3.24-9) is the homogeneous solution in which 
the atmosphere undergoes free oscillations. These kinds of waves, also known 
as free traveling planetary waves or Rossby waves[ Andrews, Leovy, Holton, 
1987], can exist throughout the lower and middle atmosphere in the absence 
of external forcing and have been well-documented. They appear in surface 
pressure data [Madden, 1978], upper tropospheric radio-sonde data [Mad- 
den, Labitzke, 1979; Venne, Stanford, 1979] and in stratospheric satellite 
data [Rodgers, 1976; Hartmann, 1976]. 

The most prominent of all of the atmospheric free oscillations is the west- 
ward traveling 5-day period wave with zonal wavenumber a = 1. Other 
traveling modes have also been observed, among them are the a = 1, 16-day 
wave [Madden, 1978; Madden, Labitzke, 1981; Madden, 1983] and the a = 3, 
2-day wave [Salby, 1981; Rodgers, Prata, 1981]. Figure 3.5.1 shows the lati- 
tude variations of the temperature perturbations from the mean state for the 
5-day component determined by Rodgers (1976). There are amplitude max- 
ima of about 0.5 K near 50° N/S. Figure 3.5.2 shows a five day component 
of the pressure variation as observed in the stratosphere over the surface of 
the earth [Hirota, Hirooka, 1984]. The patterns show the a = 1 character 
of the waves with westward traveling surfaces of constant phase. The wave 
structure is symmetric about the equator and has amplitude maxima near 
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50° N/S latitude in agreement with Rodgers (1976). 

The existence of a wave with a period of 5 days and zonal wavenumber 
a = 1 in the earth’s atmosphere is in agreement with the theory of free 
planetary modes [Salby, 1979] and is a homogeneous solution of the linear 
wave theory §111.4. To determine the nature of the solutions, we solve (3.57) 
in the absence of heating W = 0 [Salby, 1979] 



(3.67) 


after defining the auxiliary function 


AP 

FA 3 ' 


(3.68) 


The homogeneous equation (3.67) is considered to be an eigenvalue problem 
for the parameter A [Salby, 1979; Andrews, Leovy, Holton, 1987] subject to 
appropriate boundary conditions. The value(s) of A depends on the verti- 
cal temperature structure through the function u>. The frequency number <r 
which corresponds to the eigenvalue A for each given a are then computed 
from the dispersion relations portrayed in Figures 3.4.1-7. This method pro- 
vides the frequency for the free oscillations. The latitudinal Hough functions 
Q£* are then determined using the method in Appendix C. In Figure 3.5.1 the 
Hough function for the 5-day wave is plotted along with the observed tem- 
perature fluctuation amplitude [Rodgers, 1976] and the similarity between 
the peaks near 50° N /S and the minimum at the equator is striking. 

A single value of eigenvalue A for free oscillations of the atmosphere of 
about 8.88 was determined shown to exist for free atmospheric oscillations 
from the ATMOS standard temperature profile in Figure 3.5.3 by the method 
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described in Appendix D. The corresponding frequencies of oscillation can be 
determined from the graphs of the dispersion relations (see Figures 3.4.1-7). 
The value of the Lamb mode eigenvalue A = 8.88 is plotted on the a = 1 
dispersion relation in Figure 3.4.5. The point marked ® corresponds to the 
eigenfrequency of <r = 0.099 which corresponds to a period of about 5.1 days. 
Since this wave is solution to the homogeneous problem it can propogate 
without the need for external forcing. The 5-day wave is believed to be 
excited by random disturbances in the atmosphere and fluctuations of mean 
winds [Andrews, Leovy, Holton, 1987] and can persist for weeks at a time 
[Salby, 1979]. The 5-day mode as well as the other homogeneous planetary 
waves are ‘simply the result of the hydrostatic distribution of mass in the 
atmosphere, the earth’s rotation, and spherical geometry’ [Salby, 1981]. 

We can study how free traveling planetary waves can influence observed 
wind speeds. The vertical structure equation (3.67) is homogeneous for free 
oscillations and its solutions are unique to a multiplicative constant. This 
multiplicative constant can be determined by adjusting the amplitudes of the 
surface pressure oscillations predicted by the theory to fit observed ampli- 
tudes. An analysis by Burpee (1976) indicates that the component of the 
surface pressure oscillation identified with the 5 day wave is slightly less than 
1 mb near 10-20° N latitude and some of the surface pressure data in that 
study is shown in Figure 3.5.4. The five day oscillating component is evi- 
dent in the figure. Using these observations, the amplitude of the computed 
pressure fluctuation from the linear wave theory §111.4 was adjusted to the 
experimental values to have an amplitude of 1 mb at 15° N. The fractional 
deviation of the pressure from the hydrostatic value (given by P' / = eFAP 
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in Appendix A) as function of the number of scale heights x and fi (= sin#, 
$ = latitude ) is shown in Figure 3.5.5. The amplitude P'/P^ remains below 
10 % until about * = 17 which corresponds to about 120 km altitude. The 
corresponding zonal wind speed u is shown in Figure 3.5.6. The zonal winds 
are less than 1 ms -1 below * = 17 or 120 km altitude at all latitudes for the 
5-day wave. 

We therefore do not expect planetary waves to significantly influence the 
zonal wind speeds at mesospheric heights and below and conclude that this 
kind of atmospheric disturbance is not likely to be responsible for the large 
variations in winds measured during the ATMOS mission. We must therefore 
look for other effects. In the next section we show that the nonhomogenous 
solutions to the linear wave theory §111.4 which are driven by solar heating 
of the atmosphere can indeed influence the circulation of the atmosphere in 
a significant way and is likely to be responsible for the wave-like features 
observed consistently with the ATMOS instrument. 
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Figure 3.5.1 Solid curve: temperature amplitude of the 5-day wave as a func- 
tion of latitude at about 42 km altitude for November 1973, as measured by 
the Selective Chopper Radiometer on the Nimbus 5 satellite [after Rodgers, 
1976]. Broken curve: latitudinal structure of the Hough function correspond- 
ing to the 5-day Rossby mode [from Andrews, Holton, Leovy, 1987]. 
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Figure 3.5.2 The 5-day wave at 1 mb or about 45 km as observed by the 
Stratospheric Sounding Unit on the TIROS-N satellite, for 6 successive days 
in August- September 1980 [Hirota, Hirooka, 1984]. The wave-number 1 
Fourier component of the geopotential height anomaly (which can be directly 
related to pressure - see Houghton, 1987). The westward-travelling pattern 
except south of 50° S. Shaded areas denote negative anomalies; contour levels 
are 20 m [from Andrews, Holton, Leovy, 1987]. 




Figure 3.5.4 Time series of surface pressure averaged for groups of stations 
in 5° longitude bands from 10°-20° N in the Caribbean [from Burpee, 1976]. 
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Figure 3.5.5 Fractional deviation of the pressure from the hydrostatic back- 
ground pressure as a function of number of scale heights x. Notice that 
the deviation remains below 10 % until x ~ 14 which corresponds to about 
120 km. This type of oscillation does generate significant winds until great 
heights are reached (see Figure 3.5.6). Notice the two peaks at 50° N/S which 
agrees with the data presented in Figures 3.5.1-2. 



Figure 3.5.6 Computed zonal wind speeds as a function of number of scale 
heights x for the 5-day Rossby wave. Notice that the deviation remains below 
5 ms -1 until * ~ 18 which corresponds to about 140 km. 
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III. 6 Atmospheric Thermal Tides 

In this section we show that forced tidal oscillations, due to the peri- 
odic heating and cooling of the atmosphere, can give rise to large pertur- 
bations of the pressure, density, temperature and velocity about the mean 
background state. Atmospheric thermal tides axe forced global-scale daily 
oscillations. They are driven by the diurnally varying heating of the earth’s 
atmosphere due to ultraviolet radiation absorption by water and ozone [An- 
drews, Leovy, Holton, 1987] between the surface and up to about 75 km 
[Chapman, Lindzen, 1970; Groves, 1982a- b]. 

The energy absorbed by the atmosphere is transformed into mechanical 
energy which is transported to great altitudes [Chapman and Lindzen, 1987]. 
Energy is also radiated away by infrared cooling (also known as radiative 
damping) primarily by the 15 fii n band of COj, [Dickinson, 1984] which 
tends to restore the atmosphere to its undisturbed background temperature 
distribution [Chapman and Lindzen, 1970; Kato, 1980]. The first mechanism 
is usually incorporated into the boundary conditions of the vertical structure 
equation (3.57) described in Appendix D. Radiative damping is contained in 
the heating function 7i and its inclusion in the linear wave theory §111.4 is 
described in Appendix J. 

The diurnal heating of the atmosphere with respect to local time can be 
represented by a distribution similar to Figure 3.6.1. The heating vanishes 
on the ‘dark’ side of the planet while it has a maximum near local noon 
and decreases to zero at local sunrise and sunset. A Fourier analysis of 
this distribution yields a steady component, a diurnal component (with a 



217 


24-hour period), a semidiurnal component (with a 12-hour period), and so 
on. In order to estimate the atmospheric response a knowledge of the Fourier 
components of the heating is necessary. The phase of the heating ‘follows’ the 
sun [Seibert, 1961; Chapman and Lindzen, 1970; Andrews, Leovy, Holton, 
1987] and it is easily verified that the frequency number a and the zonal 
number s in (3.30) are equal in this case. It can be verified that the phase 
of the heating and the linear response travels westward and hence the term 
migrating tides is used to describe this class of linear waves. The component 
<r = 8 = 1 is called the diurnal tide and the component a = s = 2 is called the 
semi-diurnal tide. The dispersion relation for the eigenfunctions of Laplace’s 
tidal equation are plotted for the diurnal and semi-diurnal modes in Figures 
3.6.5-6 respectively. 

For each tidal mode described by <r, a complete set of Hough eigenfunc- 
tions O n and eigenvalues A n can be determined by the method in Appendix 
C. The corresponding vertical structure equation (3.57) can be solved and 
the local time, latitude and height dependence of the tidal modes can be 
determined §111.4 provided the Hough components of the heating term Ti are 
known. 

In Figure 3.6.2 is shown theoretical northerly wind speeds [Lindzen, 1967] 
for a hydrostatic background model and rocket measurements [Reed, 1969] 
for the diurnal component. The amplitudes of these two sets of data are 
similar, however, the hours of maximum amplitude vary with altitude. The 
discrepancies between the phases with altitude are probably due to the ne- 
glect of mean winds, errors in the heating and perhaps atmospheric cooling 
in the theoretical model as pointed out by Lindzen and Hong’s calculations 
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(1974) for the semi-diurnal components. In the lower atmosphere (below 60 
km) the linear tidal theory can be used to describe wind speeds since the 
amplitudes of the waves are small (horizontal speeds < 10-20 ms -1 ) and the 
linearization of (A.61-7) is a good approximation. 

To further understand the tidal phenomena, we can calculate the solution 
to the linear theory using a model heating function and temperature profile. 
The ATMOS standard temperature profile in Figure 3.5.3 was used to model 
the hydrostatic background. The heating model given by Groves (1982a- 
b) was chosen since it includes results for diurnal and semi-diurnal heating 
components of both water vapor in the boundary layer and stratospheric 
ozone. The forcing model data describe the vertical structure of the heating 
Tin for each Hough component (see 3.55). The model also describes the 
heating components as function of the season. Since the ATMOS mission 
was flown near the vernal equinox the heating profiles corresponding to this 
time of the year was used. Figures 3.6.3 and 3.6.7 show the latitude and 
height dependences using (3.55) for the diurnal and semi-diurnal components 
respectively as given by the model [Groves, 1982-a-b]. 

The solutions to the linear theory §111.4 was applied to the diurnal and 
semi-diurnal modes using the heating rates by Groves (1982a-b). The calcu- 
lated temerature and horizontal velocity components in Figures 3.6.4-6 and 
3.6.8-10 for the diurnal and semidiurnal components respectively. These re- 
sults show the exponential increase of the amplitudes of the physical fields 
with scale height x as described by (3.60-6). From these figures we see 
that above x = 10 or about 70 km the upper atmosphere is extremely ac- 
tive mechanically. Dissipative mechanisms such as molecular viscosity in the 
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lower thermosphere and infrared cooling by COj [Dickinson, 1984] in the 
stratosphere and mesosphere along with non-linear wave breaking near the 
mesopause [Fritts, 1984] and subsequent generation of turbulence play a role 
in damping the amplitudes of these waves. Radiative damping using the 
Newtonian cooling approximation [Chapman and Lindzen, 1970; Andrews, 
Holton, Leovy, 1984] can be added to the linear wave theory §111.4 with 
little added complexity to the integration procedure (see Appendix J) and 
was used in my calculations. Lindzen and Hong (1974) have modeled mean 
winds as well as molecular viscosity and radiative damping effects in their 
calculations of semi-diurnal tidal perturbations and have found that velocity 
amplitudes reach a maximum value of about 100 ms -1 about 150 km. 

Although the inclusion of radiative damping into the linear wave theory 
dampens the amplitudes by about 20 % [Chapman and Lindzen, 1970], it is 
not sufficient to prevent the growth of the waves with height. Other mecha- 
nisms become important. Large amplitudes of the potential temperature ■& 
near 70 km [Andrews, Holton, Leovy, 1987] can cause local adiabatic insta- 
bilities since the static stability N 2 becomes negative and the local buoyancy 
frequency N becomes imaginary (see Appendix B). The atmosphere can then 
no longer support the wave and local turbulence and subsequent cascade of 
energy to smaller length and time scales occurs [Lindzen, 1981]. At about the 
same altitude where adiabatic instabilities set in, non-linearities also become 
significant [Andrews, Holton, Leovy, 1987]. When the horizontal velocity 
become of the order greater than 0.1 (in the dimensionless representation of 
Appendix A) or about 50 ms -1 then advective instabilities due to the u • Vu 
terms in the fluid equation §111.2 at wind shears become important [Wein- 
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stock, 1976]. The linearized approximation §111.4 which was used to compute 
the atmospheric response to the diurnal heating of the atmosphere is then 
no longer adequate for describing the dynamics. 

From the results we conclude the linear theory calculations of the tidal 
response to diurnal atmospheric heating suggest that the zonal winds above 
70 km observed during the ATMOS mission are caused primarily by the tidal 
waves propagating into the thermosphere. The waves in this region break 
due to non-linearities and static instabilities generated by tidal perturbations 
[Lindzen, 1981; Andrews, Holton, Leovy], Non-linearities in the theory be- 
come important and make the problem complex and hence difficult to treat 
theoretically. I have currently reached this point in my study of atmospheric 
tides. My study of the ATMOS wind profiles is by no means completed and a 
more general method for incorporating the non-linear and adiabatic breaking 
of the tides needs to be developed. 
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Figure 3.6.1 Variation of solar heating at a point in midlatitudes (heavy 
curve). The horizontal line represents the diurnal or zonal average of the 
heating. This is assumed to be balanced by a zonal-mean infrared cooling so 
that the net heating vanishes. The term H is represented by the departure of 
the heavy curve from the dashed curve [from Andrews, Holton, Leovy, 1987]. 



Figure 3.6.2 Observed northerly diurnal wind speeds in the stratosphere de- 
termined from rockets launched at Ascension Island (8°S) and theoretical 
calculations by Lindzen (1967) based on classical tidal theory [Reed, 1969]. 
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Figure 3.6.3 Diurnal heating rates per unit mass H as a function of scale 
heights x and latitude 9 used as the driving term in the tidal model for 
equinox conditions [adapted from Groves, 1982a-b]. There is heating due of 



Figure 3.6.4 Amplitude of the diurnal temperature perturbation determined 
using the linear wave theory as a function of scale heights x and latitude 
9. The exponentially increasing amplitude of the temperature perturbation 
with scale height x is evident. 
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Figure 3.6.5 Amplitude of the diurnal zonal velocity perturbation determined 
using classical the linear wave theory as a function of scale heights x and lat- 
itude 6. This theory predicts that zonal winds can grow to great amplitudes 
> 100 ms' 1 or more at mesospheric heights and propogate into the upper 
atmosphere. The great amplitudes determined here do not occur in the at- 
mosphere, however, since nonlinear and adiabatic instabilities break the wave 



Figure 3.6.6 Amplitude of the diurnal meridonial velocity perturbation de- 
termined using classical the linear wave theory as a function of scale heights 
x and latitude 6. 


224 



Figure 3.6.7 Semi-diurnal heating rates per unit mass % as a function of 
scale heights x and latitude 6 used as the driving term in the tidal model for 
equinox conditions [adapted from Groves, 1982a-b]. There is heating due of 
water vapor near the surface x = 0 and ozone in the stratosphere 4 < x < 12. 



Figure 3.6.8 Amplitude of the Semi-diurnal temperature perturbation de- 
termined using the linear wave theory as a function of scale heights x and 
latitude 8. The exponentially increasing amplitude of the temperature per- 
turbation with scale height x is evident. 
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Figure 3.6.9 Amplitude of the semi-diurnal zonal velocity perturbation deter- 
mined using classical the linear wave theory as a function of scale heights x 
and latitude 6. The semi-diurnal winds become large > 100 ms -1 at greater 
altitudes than the diurnal winds, however, the exponential growth of the 
wave will eventually force it to break and prohibits its amplitude growth 
with scale height x. 



Figure 3.6.10 Amplitude of the semi-diurnal meridonial velocity perturbation 
determined using classical the linear wave theory as a function of scale heights 
x and latitude 0. 
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SUMMARY 

I have presented many of the aspects of spectroscopy which were impor- 
tant for obtaining wind speeds in the upper earth’s atmosphere. The primary 
objective was to determine whether or not such measurements were indeed 
possible considering the technical limitations of the ATMOS interferometer. 
These have been discussed in some detail, primarily as an instructive guide 
for future research with ATMOS or other similar instruments. 

During the course of the research I have determined that the instrument 
is adequate for obtaining accurate wind speeds to precisions of 5 m/s. The 
technique described, however, is somewhat limited in its ability to obtain 
as much wind data as possible during the course of what may be a typical 
ATMOS mission due to other priorities. The lack of instrumental lines avail- 
able for calibration in all but Filter 3 means that only a fraction of the total 
number of occultations are available to obtain wind data. 

The calibration analysis indicated that the instrument was stable to bet- 
ter than 2 x 10 -s cm -1 from spectrum to spectrum and does not reveal 
significant fluctuations. This stability indicated that the instrument was 
providing consistent frequencies. Comparison between positions determined 
from atmospheric lines via the instrumental calibration indicated that the 
ATMOS data can be used to improveme accepted frequency standards, par- 
ticularly those line positions where the accuracy and/or precision are not well 
known. The atmosphere itself provides a natural laboratory for spectroscopic 
measurements of trace gases, many species which are difficult to observe on 
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earth laboratories. 

The instrumental lines which were present in the spectra indicated that 
the amount of C0 2 present tended to decrease in time. This information was 
consistent with the theory of outgassing by the instrument. The equivalent 
width of the H 2 0 lines appeared to remain relatively constantly throughout 
the mission. This information in addition to a Doppler shift of about +15 
ms -1 away from the instrument suggests that the water vapor is outgassing 
from a source primarily external to the instrument while the C0 2 is primarily 
located in the instrument, since its concentration decreased with time and 
no Doppler shift was measured. 

The winds obtained from the ATMOS spectra showed a similarity among 
the sunset occultations revealing consistent phase shifts and amplitudes. The 
source of this activity was studied from a theoretical point of view and the 
concept of the thermal wind and atmospheric tides was discussed. The study 
showed that models of the radiative transfer between the sun’s radiation, 
ozone and water vapor drive upper atmospheric circulation in a manner which 
can be considered to be consistent with the results of the ATMOS experiment. 
The regularity of winds measured during the sunset occultation indicates that 
mesospheric winds might indeed be induced by the diurnal and semidiurnal 
tides. The theory presented was found to be inadequate in describing the 
detailed structure of the winds since the amplitudes of the waves predicted by 
the theory grow exponentially with height due to the exponential decrease of 
the hydrostatic background density. The presence of an adiabatic instability 
in the theoretical model was shown to occur at about twelve scale heights or 
about 85 km above the earth’s surface. 
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The breakdown of the classical theory suggests that a more elaborate 
study into the nonlinear evolution of the tides at heights comparable to those 
studied by ATMOS is necessary. The response of the atmosphere to the 
heating of the ozone layer on a global scale suggests that upper atmospheric 
winds may contain much information regarding the latitudinal distribution of 
this species. A nonlinear theory of tides is probably need to better understand 
upper atmospheric wind patterns. The the vertical structure of atmospheric 
winds at each point in local time may contain information regarding the 
global distribution of the heating sources. 


229 


Appendix A. Scaling the Nonlinear Equations of 

Motion 


The dynamical equations in §111.2 can be used for many meteorological 
problems. If we wish to use them to describe motion over an entire sphere 
the dynamical equations (3.7-12) are used. The horizontal length scale is 
assumed be of the order of the radius of the earth a while the vertical scale 
is of the order of 100 km. The time scale will be of the order of the rotation 
period. Due to the latitude dependence of the Coriolis force oc 2D sin 0, it 
can vary significantly over the horizontal scale. 

Due to the complexity of the fluid equations in spherical coordinates 
it useful to introduce several changes of variables. The latitude 9 can be 
transformed into the new coordinate 


/x = sin 9 


(A.1) 


where — 1 < /i < 1. Then the derivatives 




(A.2) 


can be determined. It is also useful to redefine the horizontal velocity func- 
tions [Margules, 1893] as 


u' = 

yjl - fl 2 ti, 

(A.3) 

v' = 

Jl- (l 2 V. 

(A.4) 
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By substituting these definitions into the total time derivative (3.12) 

d a , a l / , a u' a \ , 4 % 

Dt at +w dz + o v dn + i ~p?a<j>) ’ ^ A ' 5 ^ 

The continuity equation (3.10) becomes 

Dp f dw 2w if dv' 1 du'\ \ 

Di + ' , (sr + T + i(^: + rr^8^)} = 0 ( A - 6 > 

and the horizontal momentum equations (3.7-8) are 


Du ' wu' . / 1 dP / J-*. 


(A.7) 


Dt a a 1 — p 2 ' 1 pa dp. v p 

(A.8) 

The vertical momentum equation (3.9) is 


- - («" + f n ) - 2(J«' = - g + 

a l-/i J ' ' pdz p 


(A.9) 


We now wish to examine dynamics on the scales appropriate for global 
scale disturbances. Atmospheric oscillations can be classified by the time 
scale and horizontal distance scale over which the motions vary. In the case 
of planetary and tidal waves, the characteristic horizontal distance is of the 
order of the radius of the earth a. Since the horizontal extent of these kind 


of atmospheric disturbances is as large as it is, the Coriolis acceleration will 
vary significantly over the range of latitudes in which these kinds of global 
disturbances propagate. 

The problem is mathematically manageable by following a scaling tech- 
nique similar to Pedlosky (1987). It is useful to introduce a change in no- 
tation. Dynamical fields will be represented in both dimensioned and non- 
dimensioned form, a “ * ” superfix representing the dimensional form. For 
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Here, u\ v\ w, z, and t are dimensionless variables. The condition that the 
dimensionless velocities are much less than 0(1) will be used in many of the 
arguments. The r* and t* derivatives become 


d U 8 -n 8 

dt * ~ adt~ U dt ’ 
d _ d _ 1 d 

(A.12) 

dr* dz* D dz ' 

(A.13) 

The total time derivative D/Dt* scales as 


° =* D D 
Dt* Dt \a J Dt 

(A.14) 

where 

D d d , d v! d 

EE 1- W )- V 1 

Dt dt dz dfi 1 — fi 3 d<f> 

(A.15) 

in terms of the dimensionless variables. 


The dynamics are superimposed on a hydrostatic background state. Let 
the pressure Ph(z) and density Ph(z) represent this “standard” hydrostatic 
state and further assume that they depend on z* alone. Since we are inter- 
ested in the magnitudes of the variables the neglect of horizontal temperature 

gradients is unimportant in the scaling transformations. 

The hydrostatic 

condition requires that 

dP h 
dF = 

(A. 16) 

Let the background be characterized by a temperature 2h(z) and a constant 
mean molecular weight M. The background pressure and density are related 

by the ideal gas law (3.4) as 


n , _X Ph(z)RTh(z) 

“ M ■ 

(A.1T) 
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Defining the scale height as 


" RT £ ] 

(A. 18) 

then 


Pk = Ph.gH‘ 

(A. 19) 

The hydrostatic condition (A. 16) can then be written 


dPh Ph 

dz * H ' 

(A. 20) 


Integrating this equation gives the background pressure field in terms of the 
surface value Ph(0) and the scale height H(z*) 

Pi(z') = Ph( 0)e~ S' *C/H(C) ( A . 21 ) 


It is useful to introduce a dimensionless inverse scale height by 

A(z) = (A.22) 

The hydrostatic pressure is then 

P h (z) = Pk{ 0)e-/o'^ A( O (A.23) 

in dimensionless units. By using the ideal gas law the hydrostatic density 
becomes 

Ph(z) = ^—A(z)e~^o d( A M. (A.24) 

The pressure P* and density p* must now be scaled appropriately and the 
quasi-geostrophic scaling given by Pedlosky (1987) is useful. In this scaling 
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we define the total pressure and density as the sum of the hydrostatic value 
and a departure from this value given by “prime” notation 


P* = P h (z) + P'(z,0,<M)> (A.25) 

P * = Ph(z) + p’(z,6,<t>,t). (A.26) 

The condition of quasi-geostrophic balance requires that the horizontal pres- 
sure gradients should be of the same order of magnitude of the Coriolis force. 
This can be written as 


p* 2ftu* = 0 ( 2ntfph), 


(A.27) 


and this should be of the order P'/a. Then, P' = 0{2SlUapk). So the natural 
scaling for P* is 

P* = Pk + Ph • 2 ClUaP (A.28) 


where P is dimensionless. To scale the density the vertical pressure gradient 
associated with P' is 



This suggests that p' be scaled as 



Therefore we write 



By defining the Rossby number 

U 

2 fla’ 


(A.29) 


(A.30) 


e = 


(A.31) 
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which is 1/2 since U = f la using this particular scaling, and by defining the 
parameter 

„ _ (an «)* 

~ gD ’ 

the pressure and density can be written by using the hydrostatic condition 
(A. 19) and (A.22) as 



P * = i^(l + eFAP ) , (A.33) 

P* = Ml + el'V). (A.34) 


With the scaling parameters used above F « 0.9. 

By using these scaling laws the continuity equation (A.6) can be deter- 
mined in terms of these dimensionless variables and the hydrostatic condition 
to be 

g + j-LjQ} =0. (A.35) 
The horizontal momentum equations (A. 7-8) become 


+ Atou'l — pv' + A/i0. — p? w = 


1 8P 


1 + eFp d<j> 




(A.36) 


e {^T + w + rr^(“° + ' , °)} + '‘ u ' = 

The vertical momentum equation (A.9) is 


_ ( 1 -^ 2 )^ 

1 + eiV ' 9/i 


(A.37) 



A p_i» (AP) _, 
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+ * 5 .. 

2SIUPI 


(A. 38) 



Figure A.l Spherical coordinate system. 


Equations (A. 35-8) represent the continuity and momentum equations 
respectively. However, the system of equations cannot be solved without 
the thermodynamic relationship (3.11). This equation must be scaled in an 
appropriate manner consistent with (A. 35-8). From (3.6) we have 


jjmjo) ( p* y 
Rp • \Pk(o )) ' 


(A. 39) 


Then 


Ini r = -In 

7 L p n(0) 


(1 + ePAP) 


— In 


Rph 

MPni 0 ) 


(1 + eFp) 


= lndhH — ln(l + ePAP) — ln(l + eFp) 
7 
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~ lnd^H — eFAP — eFp 
7 

— lndj, + In ^1 + —eFAP — eFp'j . 

By using the approximation ln(l 4- y) ~ y we have 

r = t? h ^1 + ^eFAP - . 

The fractional change of the potential temperature is 

sr „ (i 4 „ \ 

If we define the dimensionless potential temperature perturbation 

d = -AP - p 

7 

then the local potential temperature (A.41) becomes 

r = d h (l + eFd). 


(A.40) 


(A.41) 


(A.42) 


(A.43) 


(A.44) 


The buoyancy frequency, described in Appendix B, of the hydrostatic 
state can also be defined by 


Nl = f 

dh dz * 

where the hydrostatic potential temperature is defined by using (3.6) 

i/t 


(A.45) 


gjvg) ( p» y /7 
h_ P*R {pm) * 


(A.46) 


By combining the hydrostatic condition (A. 19), (A.22) and the relationship 

1 dph _ f _ £dA\ 

Ph dz* D \ Adz) 


(A.47) 
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with (A. 46) we obtain 


JW2 = i f«A - I . 

D \ A dz ) 


It is also useful to define dimensionless hydrostatic buoyancy term 

A dz 


so that 


J*2 = §v. 


(A.48) 


(A.49) 


(A.50) 


This relationship can be substituted into (3.11) and using the scaling laws 
in the previous two sections it can be shown that 


eF^ + C 1 + J^w = (l + eF0)^^-. 
Dt gJJil 

If we define the dimensionless heating rate per unit mass 

... «A7f* 


(A.51) 


(A.52) 


and use the relationship (A.43) the thermodynamic relationship can be ex- 
pressed in the form 


eF-— + (1 + eF6)J tf 2 w = (1 + cFi?)7f. 


(A.53) 


The fluid equations can be expressed in a still more fundamental form by 


a transformation of the vertical coordinate 


= r dC A (C). 

Jo 


(A. 54) 


Then 


8 _ dx d _ d 
dz dz dx dx 


(A.55) 
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This vertical coordinate allows us to express the hydrostatic background 
pressure (A.23) as 

Ph(x) = J\(0) e— . (A.56) 

The hydrostatic density (A.24) becomes 


Pht*) = ^ A ( z ) e '*• 


gD 


(A.57) 


The hydrostatic condition is then 

dP h 

dx 


= -Ph, 


(A.47) becomes 


A d/»h 


Ph 


dph. _ f . dA \ 

dx \ dx )' 


and the dimensionless hydrostatic buoyancy frequency (A.49) is 

dA 


AP = k\- 


dx' 


(A.58) 


(A.59) 


(A.60) 


The dimensionless total time derivative (A. 15) is 

d da . a v! a 
Dt~ at + Aw dx + v dfi + i - m 3 a<j>‘ 


(A.61) 


We now have showed how the five dynamical equations behave under 
these scaling transformations. Since the scaling parameter A = D/a ~ 1/64 
is very small, terms involving it can be ignored to a first approximation. 
By ignoring terms of order A as well as frictional terms, the five non-linear 
equations (A. 35-8, 51) become 

1 du r 


«,f£ +(1+ .,,,{ A £_( A -£) 


dv' 

”’ + ^ + r 


fi 3 a<f> 


|=0, (A. 


(A.62) 
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Du' 
' Dt 


— flV 


ap 


1 + eFp d<f> ’ 

{ ik + + ” 0 )} + 1 = -TTTTpV ~ 


\p - ~(AP) - P = o, 


,D i? 


eF— + (1 + eF0) APw = (1 + eF8) H 


respectively. These five equations, along with (A.43) 


= -A P-p 
7 


(A. 63) 
(A.64) 

(A.65) 

(A.66) 

(A.67) 


are the nonlinear equation which describe global scale disturbances 


Appendix B. The Atmospheric Buoyancy 

Frequency 
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The hydrostatic background temperature can sustain periodic free buoy- 
ancy oscillations, also known as ‘gravity waves’ [Houghton, 1986; Pedlosky, 
1987]. Consider a fluid element A at the height z A and another B at zb 
where zb = z A + dz. Suppose the fluid element A is displaced upwards to 
the level zb occupied by B. Let the displacement of fluid element A over the 
distance dz be slow enough that the pressure of the fluid element A continu- 
ously adjusts to its surrounding but rapid enough so that thermal dissipation 
and external heating of the fluid element is negligible. The fluid element is 
then undergoing an adiabatic displacement. Mathematically it can be stated 
that the potential temperature of the fluid element is preserved during the 
displacement. We can say that 

dd = 0 (B.l) 


following the fluid element. From (3.23) this means that 


~dP = -dp. 

7 P P 


(B.2) 


The change in density of the element A can then be determined 
. (dp\ . 1 pdP 

Now the density of the fluid element B is 

® *• 


Pb = Pa + 


(B.3) 


(B.4) 
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The difference in density of the displaced fluid elements A and the background 
fluid elements B is 


(pa + ip a ) = f) iz - < B - 5 > 

The restoring force per unit mass caused by this density fluctuation is 

g, * \ ( 1 dP ldp\ , 

-(PA+dpA-PB) = * 

= IT,**' < B ' 6 > 

Since the restoring force is proportional to the displacement an oscillation 
with frequency N is implied 


N ' = It < B ' 7 > 

N is called the Brunt-Vaisala frequency, also known as the buoyancy or grav- 
ity wave frequency, and represents the natural frequency of vertical oscillation 
in a density stratified fluid. If N 3 > 0 the atmosphere oscillates vertically 
about an equilibrium state and the atmosphere is adiabatically stable. How- 
ever, if N 2 < 0 the displacements of the fluid elements continue to grow and 
the atmosphere is unstable. 
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Z-^Z (B) p fB) 

Jx. 

z ( 3 ) pCA> 

Figure 3.4.1 Fluid element A at z* is slowly raised to z* + dz* . 
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Appendix C. Method for Determining the 
Eigenfunctions @£ a (/x) 

To determine the eigenfunctions of Laplace’s tidal equation a decompo- 
sition described by Longuett-Higgins (1968) and Moura (1976) was found to 
be particularly useful for computational purposes. The streamline approach 
described by [Longuett-Higgins 1968; Moura, 1976] was found to be useful 
in determining the atmospheric tidal eigenfunctions and the implementation 
of this method is reviewed in this appendix. 

For convenience, we define the differential operator 

b^(i-V)A. (c.i) 

Then (3.124-5) are written as the pair of equations 

ifu' — fiv' = —iaP, (C.2) 

ifv' + fiu 1 = —DP. (C.3) 

Since 

p = '£p^ 

«' = £ Pn K, (C.4) 

n 

v' = Y.iP.K, 

n 

and by omitting the <r, a superscripts since they are understood at this point 
of the development, then it can easily be shown that for each mode n 


fU'-nV' = —a©, 


(C.5) 
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fV'-pU' = VO. 


(C.6) 


By taking 


we get the relation 


£< c - 5 >+ t^( c - 6) 




By taking 


we also find 


^(A. 6 ) + r ^(A.5) 


f(WL + -l—u\ -p(— + —?— i v\ -U' = V 2 0 (C.8) 

1 \ dp 1 - m 2 7 \ <*M 1- p 2 ) 


after defining the horizontal Laplacian 

d 3 \ d 


V 3 = — (l - i* 2 ) — . 

dp ' ' dp 1 — p 2 


(C.9) 


At this point it is useful to introduce the velocity potential $ and stream 
function 9 [Longuett-Higgins, 1968] by the transformations 


U' = -£>¥ + *$, 
V' = -V$ + s*. 


(C.10) 

(C.ll) 


Substituting these transformations into (C.7) and (C.8) we obtain the pair 
of equations 


(/V 2 + «) « - (mV 2 + V) $ = 0, 

(/V 2 + s)$- (mV 2 + V) 9 = -V 2 0. 


(C.12) 

(C.13) 



246 


At this point we have two equations for the three unknown functions 
and 0. By using the expansions (C.4) and Laplace’s tidal equation (3.149) 
in (3.129) it can be shown that 


% + T^ u ‘ = -V® < c - 14 > 

where it is understood that A is the eigenvalue. However, by using the trans- 
formations (C.10-1) 

^ = -V a $. (C.15) 

dfi 1 — n 3 

Therefore 

V a $ = A/0 (C.16) 


and 

V a © = ^V 4 $. (C.17) 

When this is substituted into (C.13) we find 

$ - [fiV 2 + V)* = 0. (C.18) 

This equation together with (C.12) constitutes the eigenvalue problem and 
is equivalent to solving the eigenvalue problem (3.149). 

In order to solve the system of differential equations (C.12) and (C.18) 
we expand the functions $ and 9 over an orthogonal set. The “standard” 
representation of the associated Legendre functions satisfy the eigenvalue 
problem 

V a P^ = — m(m + 1)P« . (C.19) 


(_L v . + /v> + .) 


The “standard” representation is normalized as 

2 (to + «)! 


/ 1 in PM = 


2 to + 1 (to — s)! 


(C.20) 
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and satisfies the recursion relations 


VP m 


m + S p. , m^a + 1 
2m + l ro " 1 + 2m +1 m+1 ’ 

(m + l)(m + j) m(m - a + 1) 

2m + 1 m_1 2m + 1 m+1 


(C.21) 

(C.22) 


It can be shown from these recursion relations that 


(pV’ + v) P' 


(m - l)(m + l)(m + <) 

2m + 1 m_1 

m(m + 2)(m - s + 1) 

2m + 1 m+1 * 


(C.23) 


We use the expansions 


* = (C.24) 

m>« 

» = E PA- (C.25) 

m>» 

With this representation and the above recursion relations then after a series 
of algebraic manipulations that (C.12) and (C.18) become the system of 
equations 


where 


+ Pm+l^m+1 + — 0> 


La + Lm 


+ Pm+l^m+1 + — 0 


(C.26) 

(C.27) 


T — 

P f* 

(C.28) 

"m — 

m(m -f 1) m a (m + l) a ’ 


/(m + l)(m + j) 

(C.29) 

Pm “ 

(m — l)m a (2m + 1)’ 

9m = 

/m(m — 5 + 1) 

(C.30) 

(m + l) a (m + 2)(2m + 1)’ 
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Solving for \? m in (C.26) and substituting into (C.27) we get the equation 


^ 4?m Rm lm+J^m+2 Qm— 2^m— 2 — 0 


where the auxiliary symbols 

Rm 

Qm — 

Pm = 


= -(A»- 


Pm+l9m Pm9m— 


■'m+1 


'll — 1 1 

-i r 


9mQm+l 

L m+ i ’ 

PmPm—l 

Lm-1 


are used. Then (C.31) becomes two sets of equations 



( R' 

Ps+2 

0 

0 




' #. ' 


Q. 

Rf+2 

^.+4 

0 

. . . 



$«+2 










(-a) 1 - 

0 

Qi+2 

R$+ 4 

Pt + 6 




*M-4 


0 

0 

Q*+ 4 

Rt + e 




$.+6 

• 

V 

• 



'■ / 

- 


k : > 

( 

■R«+i 

P i+3 

0 

0 


1 

1 

( *.+. 


o- 


Qt+1 Rm + 3 P *+ S 0 
0 Q,+ 3 iZ«+s -F»+r 

0 0 -R»+7 


•y 


$•+3 

$«+5 

$*+7 


V 


(C.31) 


(C.32) 

(C.33) 

(C.34) 


= 0, (C.35) 


= 0 (C.36) 


for the symmetric and anti-symmetric solutions respectively. If a ± super- 
script denotes these two classes of solutions. Then (C.35-6) can be written 
as 

[e ± X-X :t ]$ :fc = 0. (C.37) 
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This is an eigenvalue problem where A* = — l/e ± and A ± is the matrix in 
brackets. The eigenvectors $ ± are the expansion coefficients in terms of the 
standard representation of the associated Legendre functions. However, in 
the form (C.35-6) the matrix A ± is not symmetric. Therefore the eigenvectors 
will not be orthogonal. This is due to non- normalized form of the standard 
representation P^. From (C.16) and by using an expansion 


0 = £ 


m> J 


and by using (C.19) we have 

0 ± _ 

^ ±f * m* 

This can be written more compactly in matrix form as 


0± = --r-Ad ± $ ± 


A± /' 


where M± are the diagonal matrices 

l M. 0 

M + = 


0 M m +i 0 

0 0 M .+ 4 




Mr = 


0 0 
0 M .+ 3 0 

0 0 M t+6 


(C.38) 


(C.39) 


(C.40) 


(C.41) 


(C.42) 
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and the matrix elements axe given by 


M m = m(m + 1). 

Then (C.37) can be transformed into 

[e* J - M ± A ± M ± ~ 1 ] Q ± = 0. 


(C.43) 


In this form the matrix is still not symmetric because the standard represen- 
tation is not normalized. Define the ‘'normalized representation” by 


P' = 

1 m 


2m + 1 (to - a)! 


2 (to + a)! 

From (C.20) it is trivial to show that 

J^dfl P^Pm' — fimm 1 ' 

The eigenfunctions 0 can also be expanded upon these 

© = E ©mP* 

m>* 


(C.44) 


(C.45) 


(C.46) 


where the overbar refers to the normalized representation. Define the nor 
malizing factor 

N m = 


N 


2m + 1 (to — a)! 


(C.47) 


2 (to + a)! 

so that (C.44) becomes 

Pm = NmP, ; (C.48) 

and substitute into (C.46). Using (C.38) it becomes clear that the compo- 
nents in the standard and normalized representations are related by 


©£ = JV.6*. 


(0.49) 
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In matrix form 


where the matrices // ± axe 


= 


0 ± =M’ ± Q ± 


N. 0 0 

0 N.+ 2 0 

0 0 N ,+ 4 


(C.50) 


(C.51) 


Af~ = 


N . + i 0 0 

0 N . + 3 0 

0 0 N.+ 5 


(C.52) 


Since diagonal matrices commute then in the normalized representation (C.43) 
becomes 

e ± X - (aT ±_1 M ± ) A ± (■A/' ±_1 M*) _1 ] 0* = 0. (C.53) 

The matrix (,A/’ ±-1 Ad ± ) is symmetric. It can be shown 

that 

' R, T, 0 0 

T, Ri+2 T,+ 2 0 

0 T,+ 2 Rt+ 4 X»+4 * • • 

0 0 T»+4 R t + 6 


(,A/’ + 1 A4 + ) A + (A/' + 1 M + ) 1 = 


• / 


(C.54) 
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Rt+i 

T t+ r 

0 

0 

T . + 1 

R*+3 

T .+ 3 

0 

0 

T,+2 

•fr«+5 

T ,+ 5 

0 

0 

r. +5 

Ri+7 


where 


rp ^rn+J 0 

- L m — Vr 


2m + l(m + « + l)(m + « + 2) 


(C.55) 

(C.56) 


2m + 5 (m — a + l)(m — 5 + 2) 

Since these matrices are real-symmetric then the eigenvectors 0* are or- 
thogonal with real eigenvalues e ± = — l/A* [Mathews, Walker, 1970]. The 
eigenvectors are then normalized. The phases of the eigenvectors are chosen 
so that 0 + (O) > 0 and dQ~/dfi{ 0) > 0 for the symmetric and anti- symmetric 
eigenfunctions respectively. 

Now that the eigenvectors 0* in the normalized representation have 
been determined we wish to compute t7 #± and V f± . We use expansions over 
the standard associated Legendre functions 


U* 

= V' u 1 * P‘ 

/ ^ u m *TTI> 

(C.57) 


m>i 



= E V m P ™- 

m> t 

(C.58) 


From (C.10-1) we see that U' and V' are determined from $ and \I?. Now 

(C-59) 

and since 

= -AVeJ/JH.. (c.6o) 
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therefore 

*£ = (c.ei) 

The stream function $ can be determined from (C.26) 

1 

^ m = j (Pm+l^m+1 "h • (C.62) 


After substituting the standard representations of U\ V\ $, and $ into 
(C.10-1) it can be shown that 


tti± _ , $ ± _ (m + 2)(m + a + 1) ± (m - l)(m - a) 

m TO+1 2m + 3 m ~i 2m -1 


(C.63) 


V* = .$±_$± (”* + 2)(m + s-H) ± (m - l)(m - a) 

m s *m *m+i 2m + 3 m " 1 2m - 1 


(C.64) 


The derivations in this section were used to obtain the Legendre func- 
tion expansions for the eigenfunctions 0 and eigenvalues A of Laplace’s tidal 
equation. This technique can be used for arbitrary values of / and a. The 
differential operator C is represented by a matrix and the eigenvectors and 
eigenvalues are obtained by using algebraic methods. Once the expansion 
coefficients 0 m are determined the velocity eigenfunctions U* and V' are 
determined by using recursion relations. 

This method presented in this appendix was used for the latitudinal struc- 
ture of the linear wave theory §111.4 and susequent calculations §111.5-6. 
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Appendix D. Method for Solving the Homogeneous 
Vertical Structure Equation 


In this appendix we present a method which is used to determine the 
eigenvalues and eigenfunctions of the homogeneous vertical structure equa- 
tion (3.67) for the auxiliary temperature dependent function u defined in 
(3.68). Equation (3.67) can be written in the form 

0 +*'(*;*)» = 0 (D.l) 

where 

= (D.2) 

is a complicated expression and plays the role of the refractive index [Chap- 
man and Lindzen, 1970; Salby, 1979]. In regions where k 2 > 0 we have 
oscillatory behavior of y while in regions where k 2 < 0 the behavior of y is 
evanescent and reflections can occur [Salby, 1979]. The homogeneous verti- 
cal structure equation (D.l) is solved with two boundary conditions. At the 
surface x = 0 we impose the condition that the vertical velocity w vanishes 
[Chapman and Lindzen, 1970; Salby, 1979]. Using (3.44) this is expressed 

{£+ &-§)»} 

The upper boundary condition must satisfy a radiation/finite energy con- 
dition [Chapman and Lindzen, 1970; Salby, 1979]. The radiation condition 
in the WKB approximation is discussed in Appendix E and we will use the 


= 0. (D.3) 

c=0 
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result for the upper boundary condition. From the WKB approximation, at 
some upper level x top the solution y must satisfy the boundary condition 

= { ik - iJ ir)*\ =°- < D - 4 > 

V / l*top 

This condition along with (D.3) are used to determine the eigenvalues A of 
the system for the temperature profile described by uf. 

To solve the eigenvalue problem we can use a finite difference approxima- 
tion. We will divide the region 0 < x < x top into N + 1 levels with spacing A. 
The lower boundary is given by x 0 = 0 and the upper boundary xjf — x iop . 
In (D.l) the second derivative of y is approximated by 

<?y _ Vm+i - 2 y m + y m -i /r , 

dP 2 ? ’ (D - 5) 

TO 

The differential equation (D.l) can be expressed in terms of the difference 
equation for 1 < m < N — 1 

ym-i + (—2 + A 2 fc£j y m + y m +i = 0. (D-f>) 

Using the finite difference method the lower boundary condition (D.3) can 
be expressed in the form 

,~ l + A " s)] y ° + y ' = °- (D,7) 

The upper boundary equation is derived from (D.4) and is 

VN - 1 + -1 + A (ik - J VN = 0. (D.8) 

The system of equations (D.6-8) can be cast into the form 



My = 0 


(D.9) 
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where y is the column vector of y values 


yo 

yi 

VN 


(D.10) 


and the matrix M is defined by the coefficients of the y m in (D.6-8). The 
equation (D.9) has a non-trivial solution if 


det[M] = 0. (D.ll) 

The eigenvalues of the system are those values of A in which (D.ll) is true. 

The eigenvalues A of the system can be determined by a search method 
in which the parameter A is varied and the magnitude of the determinant 
is evaluated. This technique was carried out for the ATMOS standard tem- 
perature profile described in §111.5. Trial values of A were varied from -10 
through 40 and | det [M] j vs. A with a spacing AA of 0.2 is shown in Figure 
D.l. Two ‘spikes’ were found in the data at A ~ 8.88 and 13.18. These 
minima are magnified in Figures D.2-3. Increasing the resolution of AA to 
0.05 near A ~ 9 increases the sharpness of the first peak. This indicates 
that an eigenvalue has been detected. Near the eigenvalue A ~ 13 increasing 
the resolution of AA was not found to increase the sharpness of the peak. 
Indicate that there is no true eigenvalue near A ~ 13 and I will disregard it 
in this study. 

We see that a solution to the homogeneous vertical structure equation 
exists with a unique eigenvalue A ~ 8.9 and vertical structure function y(x). 
Corresponding to this eigenvalue is a multitude of eigenfunctions to Laplace’s 
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tidal equation (3.56) with various frequencies and azimuthal wave numbers 
for free oscillation. 
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Figure D.l Graph of values of det[M] a* the parameter A was varied between 
10 and 40. There are two distinct ‘spikes’ in the data at 8.88 and 13.18 
respectively. The first value corresponds to the free oscillating mode or Lamb 
mode. The second spike is ‘approximately’ an eigenvalue and it depends on 
the temperature profile [see Salby, 1979]. 
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A 


Figure D.2 Region around A ~ 9 at higher resolution. The magnitude of the 
minimum increases by a factor of 40 compared with Figure D.l as the search 
resolution AA is improved. This suggests the presence of the eigenvalue. 



Figure D.3 Similar to Figure D.2 in the region around A ~ 13 at higher 
resolution. In difference to Figure D.2, the magnitude of the minimum does 
not increase any further with increased search resolution AA. This suggests 
that the dip near 13.18 is an artifact of the calculation and probably does 
not correspond to a true eigenvalue of the atmosphere. 
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Appendix E. The WKB Approximation 


One particularly useful approximation to the solution of equations of the 
form 

0 + *’(*)» = 0 (E.l) 

is to use the WKB approximation [Richmond, 1971; Mathews, Walker, 1970]. 
This method is valid when k(x) varies sufficiently slowly with x so that the 
change in k in one wavelength is smaller than k itself. The form of (E.l) 
suggests solutions of the form 


y oc 

(E.2) 

Then (E.l) becomes 

-OT + *^' + * 5 = o. 

(E.3) 

If we assume <f>" is small then to a first approximation 

<f> ~ ±k 

(E.4) 

so that 

<j) ~ ± J dxk. 

(E.5) 

The “smallness” of <j>" requires that 

(E.6) 

A second approximation is determined by iteration. 

<f>" ~ ±k'. 

(E.7) 
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This is substituted into (E.3) to yield 

(< j>'f ~ k 2 ± ik\ (E.8) 

The condition (E.8) gives 

a = ±* + ( e . 9 ) 

so that 

<f> ~ ± J dx k + ^\nk. (E.10) 

Then the two choices of sign yield two approximate solutions 

3 / ( ± ) = ^e ±i /' i ** (E.ll) 

V* 

the general solution being 

y = 4= [c+e*/ a<ixfc + c_e-‘ /’*'**] . (E.12) 

For the case where k 2 (x) > 0 the solutions are oscillatory. In the case where 
k 2 < 0 the solutions will have an exponential behavior. 


Appendix F. The Upper Boundary Conditions 


The vertical structure equation (3. 57) in the homogeneous region where 
H = 0 is described by (3.67). At the upper boundary at x iop we can assume 
the solutions to be represented by the WKB approximation described in 
Appendix E. In this appendix we derive the relationship between dy/dx and 
y at the upper boundary appropriate to atmospheric gravity waves. 
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First consider the case where k 3 > 0 in the WKB approximation (E.l). 
In this case there is an oscillatory boundary condition. The analysis is fur- 
ther simplified if we impose a ‘radiation condition’ that the energy flow is 
vertically upwards out of the heating region x < xs [Wilkes, 1949; Siebert, 
1961; Chapman and Lindzen, 1970]. The radiation condition is specified by 
choosing c_ to vanish. The reasons for this are discussed by Kato (1980) and 
he shows that for gravity waves, the energy propagation is opposite in direc- 
tion of the phase velocity. If c_ = 0 we have downwards phase progression 
which gives outward energy flow. The radiation condition assumes that there 
exist no energy sources outside of the integration region and that the energy 
propagates toward z = +oo from the lower atmosphere. As pointed out by 
Chapman and Lindzen (1970) this approximation neglects reflection of waves 
at greater heights which can be a source of disturbance. If this condition is 
used then our solutions above the top layer N are of the form 



y 

C + if' dxk 
= — 7=e 

y/k 

(F.l) 

and 






dy 

dx 

= X 

[, * *'l 
[* + 2 * ] »■ 

(F.2) 

For simplicity we define 







K = 

, ik' 
k + - - 
2 k 

(F.3) 

so that 


dy 

dx 





= iKy 

(F.4) 


represents the radiation condition. 
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In the case where k 3 < 0 we can write that k = tJfcj so that the WKB 
solution to (E.l) is given by 

1 — r dmk i f* dm fctl . . 

y = [c+e + c~e Jm * . (F.5) 

The two solutions represent exponentially decaying and increasing ampli- 
tudes with height in addition to that predicted by (3.52). The exponential 
growth is assumed to be unrealistic in the context of this atmospheric model 
which demands c_ = 0 unless a source of further attenuation exists at higher 
altitudes. If the solution at the upper boundary is specified by 


then 


Again we define 


we gives 


(F.6) 

(F.7) 

(F.8) 

(F.9) 


Appendix G. Integrating the Vertical Structure 

Equation 


In this appendix we describe how the vertical structure equation (3.57) 
can be solved. A method of solving second order ordinary differential equa- 
tions with boundary conditions at different points is used which follows the 
integration scheme proposed by Chapman and Lindzen (1970). Two bound- 
ary conditions are specified at the surface of the earth and at altitudes where 
the heating term is zero. At the surface we assume that the vertical velocity 
of the atmosphere vanishes. At the upper ‘free’ boundary we assume that en- 
ergy propagates outwards with no reflections. This is known as the radiation 
condition. 

The method described in the previous two appendices describes how to 
obtain eigenvalues A£* and eigenfunctions ©£' of the horizontal structure 
equation (3.56), i.e.- Laplace’s tidal equation. However, this is only half 
of the problem for understanding the three dimensional dynamics since the 
vertical structure must also be known for each eigen-mode S^* . The eigen- 
value A*' provides this coupling. The driving terms in (3.57) are the Hough 
components of the heating distribution 7C*'(x). 

A fairly technique for integrating second order ordinary differential equa- 
tions introduced by Richtmyer (1957) and described by Chapman and Lindzen 
(1970) was used. Consider an equation of the form 
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We wish to determine a numerical solution for Y(z) given the other functions 
F(x), G(z), jET(z), and J(z). Suppose we wish to integrate this solution over 
the range xo < x <xn and this range is divided into N equal subintervals of 
width A. Then the derivatives of Y at each point x m can be approximated 
by 


(S) - 


Yn-fl Y m —l 
2A ’ 

Yn+1 ~ + Ym - 1 

A a 


(G.2) 

(G-3) 


Then the differential equation (G.l) can be cast into the finite difference form 


where 


-AmYn+l + B m Y m + C m Y m _i = D„ 


A-m — Fm + 

Bm = -2F m + A 3 H m , 
Cm = F m ~ ~AG m , 
D m = A 3 J m . 


If we introduce a change of variables 


Y m — Yn+1 


then 

Y m — 1 = ®m-llm "i" fim—l • 

Thus the difference equation (G.4) can be put into the form 

y A m y D m 0m—lCm 

“ “ C„ m+1 + + 


(G.4) 


(G.5) 

(G.6) 

(G.7) 

(G.8) 


(G.9) 

(G.10) 

(G.ll) 
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(G.12) 

(G.13) 


By comparing this to (G.9) we see that 

= - ' Am (G.12) 

a D m — Pm-i C m 

Pm = ~B T T 7T-- ( G - 13 ) 

^ra T flm-lvm 

If ao and /?o are known then all of the other and (3 m can be deter- 
mined. The simplest lower boundary condition is usually chosen to represent 
a smooth spherical earth where the vertical velocity w vanishes at the sur- 
face x — 0 [Chapman and Lindzen, 1970]. From equation (3.61) this can be 
expressed as 


In terms of finite differences 


and this can be rewritten 

2/o = [l 


.‘(in)'..."- 

(G.14) 


(G.15) 


(G.16) 

we can identify 



(G.17) 

0. 

(G.18) 


From (G.12-3) all of the other ow and (3 m are determined. 

After calculating all of the a™ and /3 m the function Y can be determined 
from (G.10) provided that an upper value is given, since 

Yff - 1 = ckn-iYn + Pn-u 
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Yn- 2 = Olff-iYff-l + Ptf-2, 

Yi = a 1 Y i +p l . (G.19) 

The result (G.19) implies that we need to know the value for the upper 
boundary Ypf to complete the integration. 

For the case in which the vertical structure equation is homogeneous 
(§111.5, Appendix D) the upper boundary value is arbitrary since there is a 
free constant of integration which can have any value. The upper boundary 
dynamics - i.e. radiation condition - is imposed on the solution during the 
determination of the eigenvalue A. In the case of the inhomogeneous equa- 
tion (3.57) the upper boundary dynamics are imposed during the integration 
procedure. We use the WKB approximation (F.4) or (F.9) depending on 
whether k 2 > 0 or k 2 < 0 at the upper boundary respectively. 

For the case where k 2 > 0 we have a radiation boundary condition and 
apply (F.4) at the level N. Then (F.4) is approximately 

y”Zl?-' ~iK N y N . (G.20) 

Since, in this integration scheme, 

yif-i = ajv-i yN + Pn- i> (G.21) 

it can be shown that 

Pn - i 

y N ~ - — 

1 — ajv_ i - t A Ajv 

This results then allows a full integration of the vertical structure equation 
below the top layer. 
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For the case where k 2 < 0 we have an evanescent boundary and apply 
(F.9) to the upper boundary level. Similar to the derivation of (G.15) we 
apply this condition at the level N so that 

VN-VN-l ~_ KNyN ' (G.23) 

Again using (G.16) we find that 

^ — . 

* 1 - aw-, + A K„ 

This final boundary condition allows the vertical structure equation to be 
integrated. 

From this appendix we can integrate the second order vertical structure 
equation (3.57) consistent with two boundary conditions specifying the be- 
haviour of the atmosphere at great heights and at ground level. The upper 
boundary condition assumes that the heating of the atmosphere from below 
is eventually radiated outward by the tidal gravity waves (or, in the case of 
the homogeneous equation, any constant value will do). Although it does not 
assume sources of energy above the upper level, the inclusion of such sources 
can be built into theory by an appropriate adjustment of the upper boundary 
condition. The lower boundary is assumed to correspond to a flat surface 
which has no vertical component of motion and hence wind. Although this 
does not take into account features such as continents and possible tidal de- 
formations of the oceans, they can be assumed to be negligable to at least 
©(A) in the scale of the dynamics and of secondary importance. 




269 


Appendix H. Energy Transport in Linear Wave 

Theory 


The linear wave theory presented in can be used to describe the physical 
structure of the velocity, pressure and density perturbations to the hydro- 
static state generated due to homogeneous oscillations and inhomogeneous 
waves, particularly in the case of atmospheric thermal tides in which the 
periodic ozone and water water vapor heating of the atmosphere can drive 
upper atmospheric circulation. In this section we outline how these waves 
also transport energy through the atmosphere. In the absence of nonlinear 
turbulence and radiative dissipation the linear waves may transport the en- 
ergy absorbed by the atmosphere outwards. In this section we examine the 
energy transport of these waves. 

For notational simplicity we will temporarily omit the “ * ” superscript to 
denote dimensioned quantities. The energy flux per unit mass S is given by 
[Kato, 1980] 

S = Pu. (H.l) 

The divergence of the flux density is 

V • S = V • Pu 

= VP • u + P V • u. (H.2) 

By using the momentum equation (3.1) and continuity equation (3.2) we 
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have 


V.S— 

P [Dt \ 2 ) p 1 Dt 8 

From the thermodynamic equation (3.11) it can be shown that 

P DP 1 D (Z\_ n 

p 2 Dt 7-1 Dt\p) 

Then (H.3) can be shown to be 

D ( luP 1 P\ 1 

m + L +± V .S«W + g.u. 

Dt\ 2 7-1 PJ P 

We see that the energy per unit mass E is given by 

1 P 


E = 


+ 


7“1 P 


and (H.5) can be expressed 

DE 
P Dt 


-fV-S— pg-u = p W. 


(H.3) 

(H.4) 

(H.5) 

(H. 6 ) 

(H.7) 


The energy equation (H.7) expresses how the energy of the wave is ex- 
changed. Some of the heating is used to increase the energy density, some is 
transported, and some goes into doing work against gravity. Since the tides 
are periodic in time (and local time) we can average the energy equation over 
a sphere to obtain the average energy exchange as a function of altitude. Let 
us define the spherical average, which is also an average over one period, by 

1 t 1 t 2w 

(’•’) = — — / dp, I d(f> • • • . (H. 8 ) 

4x J - i Jo 

The spherical average of (H.7) is then 

(p^f) + < vs )- 0 >g •“) = </>«). 


(H.9) 
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The fields will be assumed to vary as §111.6 

/ = /" e~ (nt+ *>. (H.10) 

a=l 

It is also understood that we use the real parts of the dynamical fields in 
(H.7) when taking these averages. If we denote the complex conjugate of a 
function / by then 

R/ = i(/ + /'). (H.11) 

Only products of non-conjugate and conjugate pairs will survive the averag- 
ing process. 

If we write the pressure and density in the form §111.3 


P = ft(l + P'), (H.12) 

P = Ph(l + ?' ) (H.13) 

then the energy flux divergence is 


V-S = V • Pu = PV • u 4- u • VP 

= P h (l + P')V • u + u • Vft(l + P') 

= ft [(1 + P') ( V -u-|)+u. VP'.] (H.14) 

The spherical average of (H.14) is then (omitting “ <r ‘ r ” superscripts) 

(V • S) = ft (v • u - + u • VP'^ . (H.15) 


Since the fields in (H.15) are the real parts of the solutions, in terms of the 


complex functions we have using (H.10-1) 

<v-s> = ^E/>[^ + ^ + i 


dv 

dn 


1 ia 
r 1 — fi 3 
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( dw 2w 1 dv' l iff A 1 , 

+P ' [i; + T + -r^ + -rT^ U )- + ” P 1 


r 1 — /i a 
dP'* 


rl-/t s 


r fi/i r 9/i 


+tu- 


+ w' 


t dP' 


a . ■ - » , ■ • < H16 > 

We now wish to express this in terms of the dimensionless variables. From 

§111.10 

u' = V£P„1 C, 

n 

n 

to = A£7^to n 0 n , (H.17) 

n 

P' = ePA^^nOn, 

n 

/>' = cF'Ee.e*- 

n 

By using these transformations (H.16), the orthogonality of the Hough func- 
tions and (A. 14) 

(V-S) = 4EErhMJ+»!ft)]+o(A). (h.u) 

In a similar manner the power lost to gravity can be evaluated 
(pgw) = PhgeF (p'w) 

= aEE «-(/*«. + /-.»!)• (H.i9) 

8 r 0 a n 

The term on the right hand side of (H.7) represents the energy expended in 
heating density perturbations 


(P«) = P*tF{p'H) 

ft(0) ePlf l 

8 r 0 k 




(H.20) 
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The total rate of change of energy following the fluid particles can be evalu- 
ated from (H.7) 

(/>^) = ((’W ) -tar» ) -<V- s >. (H.21) 

The flux density can be computed from (S • n) where n unit vector pointing 
outward from the spherical surface at level x. 

(S - n) = + ®iA). (H.22) 

The energy sources and sinks given by (H. 18-20) can be used to compute 
the spherically averaged energy exchange of the linear wave solutions. The 
energy flow of the linear waves are given by (H.22). We see that due to the 
orthogonality of the Hough functions that the total energy terms are the 
sum of the various component terms. This allows us to ex ami ne each mode 
separately. 


Appendix J. Including Radiative Damping into the 

Linear Wave Theory 


The primary absorption processes of solar radiation in the atmosphere 
are the ultraviolet absorption of water near the surface and ozone in the 
stratosphere [Andrews, Holton, Leovy, 1987]. The 15-/xm band of COj is the 
primary source of infrared emission in the atmosphere and it provides the 
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medium by which the atmosphere radiates [Dickinson, 1984]. The hydro- 
static temperature distribution can be assumed to be produced by a balance 
between radiative and non-radiative processes — see Andrews, Holton, Leovy 
(1987) for details. We will assume that any temperature perturbation will 
lead to an increase or decrease in infrared cooling which tends to restore the 
atmosphere to its original temperature [Chapman and Lindzen, 1970]. 

Including this type of dynamics into the linear wave theory §111.4 requires 
only a slight modification. We can substitute (3.6) into (3.5) to obtain 


DT TDP M .. 
— - = k — —— + —K 
Dt P Dt C P 


(j.i) 


Since the infrared cooling opposes the growth in temperature, following 
Chapman and Lindzen (1970) we include a term into the right hand side 


DT 

Dt 


TDP M 

= '‘p-5i+c; li -' iT 


(J.2) 


where a is the infrared cooling coefficient and 6T is the deviation of the 
temperature from the equilibrium temperature. This is known as the Newto- 
nian cooling approximation [see Andrews, Holton, Leovy, 1987] and a is the 
Newtonian cooling coefficient. We can write (J.2) in the form 

DT T DP M - 

— — = k — — — (- — — 7i 
Dt P Dt Cp 

with the ‘modified’ heating 

7 k = H-^-a6T. (J.4) 

M 



When we include a Newtonian cooling term, the linear wave theory §111.4 is 
still used with the replacement of 7i — ► "ft. 


c '4 
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The vertical structure equation (3.57) requires that the heating be in 
the dimensionless representation of Appendix A. This is achieved with the 
transformation (A.52) 


7i = 


fcA7f* 

gDSl 


(J.5) 


where now we have introduced the dimensional **’ notation (see Appendix 
A). The temperature perturbation, described in Appendix K, is then 


ST = T h • eFT 


(3.6) 


where Tj, is the background temperature and T is the temperature pertur- 
bation function. Using (J.4) in (J.5) it can be shown using the definition of 
A and the vertical depth k described in Appendix A that the dimensionless 
heating is 

* kAH* _ .■ . 

U ~ gDQ. ~ aT 

where the dimensionless Newtonian cooling coefficient is defined 

a — (J*®) 

Using (K.6) we can then write (J.7) in terms of the pressure and density 
perturbation functions P and p 

kM 


U = 


gDSl 


— a (A P — p). 


(3.9) 


The dimensionless heating can be shown to be related to the vertical 
structure function y n for each Hough component of (J.9) using (3.60) and 
(3.62) 


74 

t/A„ 


<Py n 1 dA dyn _ (l J_ dA\ 

dx 2 A dx dx l 4 2A dx J ^ n 


. (J.10) 
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With this result, the vertical structure equation with Newtonian cooling can 
be derived. Substituting (J.10) for the heating function which drives (3.57) 
we obtain the following second order differential equation for the vertical 
structure function y n 



fVn _a_ 1 dA dy n | f A" n _ 1 
dx* ifF A dx dx FA a 4 


1 \ at 

= H rt e~ m l\ 

7 FA 


a /l 1 dA\ 
ifF I 4 + 2A I Sn 

(J.ll) 
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Appendix K. Linearized Temperature Perturbation 


The ideal gas law (A. 17) is a non-linear relationship between the pressure, 
density and the temperature. Using the notation of Appendix A for dimen- 
sional and nondimensional variables (A.33-4) describe the total pressure P* 
and density p* respectively in terms of the hydrostatic background P j,, pi, 
and perturbation functions P, p 

P * = P,(l + ePAP), (K.l) 

p* = pi, (1 -f eP p) (K.2) 

where e and F are scaling factors and A contains information regarding the 
hydrostatic temperature distribution with height. For simplicity, we can 
define a dimensionless background temperature perturbation function T by 

T* = T h (l + eFT) (K.3) 

where T* is the total temperature, Tj, is the background temperature. Using 
the ideal gas law 

n. P-RT- 

M 

and (A. 17) for the background it can be verified that the dimensionless per- 
turbation functions P, p and T are related to each other by 

AP = p + T + tFpT. (K.5) 

Linearizing and rearranging we obtain 



T = AP - p. 


(K.6) 
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